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Chapter 8 


SEPARATION OF TRAJECTORY INTO KEPLERIAN, INTERMEDIATE, AND 


GAS-DYNAMIC PHASES; BOUNDARY CONDITIONS BETWEEN PHASES 


8.1 Definition of Trajectory Phases 


The function of a guidance system in a vehicle entering the atmos- 
phere of a planet is complicated by the radically different operational 
environment it encounters during the course of entry. For example, 
the density of the atmosphere changes from an insignificantly small 
quantity at orbital altitudes to the dominating factor affecting the 
shape of the trajectory at low altitudes. 

It is convenient to consider the entry trajectory as made up of 
three separate phases: 

Gis Renuleaen Phase (Free-Fall Phase): 

That segment of the trajectory, at high altitudes, where gas- 
dynamic forces are insignificantly small (P=0). 

(2) Intermediate Phase: 

That segment of the trajectory where accelerations due to gas- 
dynamic terms are of comparable magnitude with other terms in 
the dynamical equations of motion. 


(3) Gas-Dynamic Phase: 


That segment of the trajectory where gas-dynamic accelerations 
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are the predominant terms affecting the shape of the trajectory. 

Figure 8.1 shows the various phases of the entry trajectory and 
gives definitions of boundaries between phases. 

In the Intermediate Phase, either a pure Keplerian description 
of the trajectory or the gas-dynamic simplifications* in the equations 
of motion are not warranted. Chapman ‘15? considers the entry trajectory 
as made up of two phases only** which are fitted together by matching 
boundary conditions. He defines the onset of the Gas-Dynamic Phase as 
that point at which drag has reduced the vehicle velocity by about 
0.01 of the initial velocity. It is shown in this chapter that the 
Intermediate Phase may span an operating band of sufficient width that 
three, rather than two, separate segments of the trajectory should be 
considered in a piecewise analysis. 

The investigation described in this chapter is unique in that 

three separate operating regions are suggested for guidance analysis; 
the boundaries between these operational phases are defined on the 
basis of the rate at which angular momentum and energy are being trans- 
ferred from the vehicle to the planetary atmosphere. For this purpose, 
a Conservation Parameter is defined which is related directly to the 
rate of transfer of energy and angular momentum. The Conservation 
Parameter may be determined in flight from measurements of specific 
force levels and flight path angle. Quantitative estimates are given 
for the specific force level and altitude at the phase boundaries 
between the three operating regimes for vehicles with various lift and 
“3 See Chapter 9 


** Keplerian and Gas-Dynamic 
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. 


Gasedynamic forces exceed threshold value detectable 
by Specific Force Measuring Subsystem 


/ Magnitude of the percentage change 
in distance of the entry vehicle from 
planet center is very much smaller 
(0.1 times) the percentage change in 
/ horizontal component of vehicle 
/ velocity with respect to inertial 
space. 





Magnitude of the percentage change 
in distance of the vehicle from the 
planet center is approximately 
equal to the percentage change in 
horizontal component of vehicle velocity 
with respect to inertial space 


PLANET 


Key 
K: Keplerian Phase 
I: Intermediate Phase 
G. Gas-Dynamic Phase 


Fig. 8.1 : Phases of the Entry Trajectory and Definition of 
Boundary Conditions 


24h, 


drag characteristics. 

The Keplerian Phase and the Gas-Dynamic Phase have been examined 
extensively in the current literature. The Intermediate Phase on the 
other hand, has been largely ignored in trajectory studies. From the 
guidance standpoint, this is a very important phase, particularly in 
the degenerate orbital entry profile where perturbations in the Inter- 
mediate Phase in the vicinity of perigee are the mechanism by which 
eventual entry is induced. The Intermediate Phase is also important in 
the direct entry profile; this transition from Keplerian motion to Gas- 
Dynamic flight spans an altitude band of the order of 20 miles in the 
case of entry into the Earth's atmosphere*. The distance flown while in 
the Intermediate Phase may be a significant portion of the total 


vehicular range. 


8.2 Conservation of Energy and Angular Momentum 


Angular momentum was defined in equation (4-1) as: 
P= RVIg lye a) 
The time rate of change of angular momentum is therefore: 
P = RVig + RVy¢g (8-1) 


Defining the following auxiliary parameter: 


Fg ee 


* Assuming constant flight path angle and constant drag coefficient. 
If the drag coefficient is reduced by a factor of $ during the 
course of this phaSe, such as might occur if the flow regime changes 
from free molecular to continuum, the altitude band of this phase 
increases to approximately 24 miles. 
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With equation (8-2), equation (8-1) is written: 
| CRS): 
P= - —~z Ng (8-3) 


Using the non-dimensionalizing procedure described in Chapter 6, 


equation (8-3) is written in dimensionless form as follows: 


p'=- agD ev (8-4) 


The Conservation Parameter was defined as follows: 


— (8-5) 


(Conservation Parameter) = 








The name Conservation Parameter was chosen because of the close relation 
of this quantity to conservation of energy and angular momentun. 

The magnitude of the rate of transfer of angular momentum from the 
vehicle to the atmosphere is the product of the Conservation Parameter, 
the rate of change of altitude, and the horizontal component of velocity 
with respect to inertially fixed coordinates. During the Intermediate 
Phase, the change of velocity is small, hence the VI¢ term in equation 
(8-4) is nearly constant in this phase. 


Total energy per unit mass was defined in equation (4-10) as: 
(aoe 9 ne 


ere awa 
<> ae + (4210) 








In dimensionless form, this was witten in equation (6-54) as: 
Veg 
ee. «1. 
E(tot) 2 oF Soa) 


Equation (6-55) showed that the rate of change of total energy may 


be expressed as: 
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E" (yan) + Eee) = ite Vy =r'f, + Vigtg (8-6) 


where f,. and tg are the radial and horizontal components of external 
specific forces in mean surface g's of the planet. 
Differentiating equation (6-54), substituting into equation (8-6), 


and rearranging gives: 


E! VW 
__(kin) = é f _I¢ £ S.. 
: my Camara (8-7) 


It was shown in Chapter 6 that for a spherical planet: 
: R 
Equation (8-8) is written in dimensionless form as: 
fy = (1-§) v'3g (8-9) 


Therefore, equation (8-7) may be written: 


EY e 
1+ __(kin) = ré te = (i-§ ) 6 | (8-10) 
E" (pot) E 5 


Equation (8-10) shows that in a Keplerian trajectory ( fr = 0, 
E —elO) : 


EB" (kin) = 7 E" (pot) (8-11) 


This corresponds to a continuous trade-off of kinetic and potential 
energy throughout the trajectory; total energy remains constant. For 
small flight path angles, the radial component of external specific 


force is approximately equal to ny and the horizontal component is 
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approximately equal to -np. The predominant term in the bracket on the 


ee) — 
right-hand-side of equation (8-7) is = f0 because r' is generally 


very much smaller than Vig until near the terminal phase of the mission. 
Therefore, in the vicinity of atmospheric perturbations where velocities 
are generally high and flight path angles are relatively small, equation 
(8-10) reduces to the following approximate form: 
ee pein) a as) vig r (8-12) 
(pot) 

Equation (8-12) is plotted in Fig. 8.2. The Keplerian Phase of the 
trajectory corresponds to € in the vicinity of +1 while the Gas- 
Dynamic Phase corresponds to € in the vicinity of zero. All other values 
On E correspond to flight in the Intermediate Phase. In the Keplerian 
Phase, the rate of change of kinetic energy is equal to the negative of 
the rate of change of potential energy. In the Gas-Dynamic Phase, the 
magnitude of kinetic energy transfer rate is very much greater than the 


magnitude of the potential energy transfer rate. 


8.3 Relating the Flight Trajectory and the Conservation Parameter 
Gas-Dynamic forces become significant, from the guidance standpoint, 
when they are of sufficient magnitude to be detectable by accelerometers 
carried within the vehicle. The trajectory is described by a pure 
Keplerian transfer ellipse as long as the angular momentum of the entry 
trajectory is conserved. Under these conditions, EG = 0 in equation 


(8-8); therefore: 


NId _ : (8-13) 


or 
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OV 
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av 
i (8-14) 
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vo 


Equation (8-14) shows that in order for angular momentum of the 
entry trajectory to be conserved, the percentage change in distance of 
the entry vehicle from the center of the planet is equal to the negative 
of the percentage change in horizontal component of vehicle velocity 
with respect to inertial space... 

Angular momentum is not conserved when external forces (e.g., lift, 
drag, and thrust forces) exist in the ly direction. The auxiliary 
parameter é was defined in equation (8-2). The magnitude of € is an 
expression of the ratio of the percentage change in radial distance to 
the percentage change in horizontal component of vehicle velocity with 
respect to fixed inertial coordinates. E must be equal to unity for 
conservation of angular momentum, 

During the course of entry, the magnitude of the ratio of percen- 
tage change in radial distance to the percentage change in horizontal 
velocity will decrease toward zero as gas-dynamic forces become pre- 
dominant factors in specifying the ultimate path of the vehicle. Hence, 
the magnitude of E is an indication ofthe degree to which angular 
momentum of the trajectory is conserved; as E becomes very much smaller 
than 1.0, angular momentum is transferred rapidly to the gaseous envelope 
of the atmosphere through the mechanism of lift and drag. 

From equation (8-9), the horizontal component of external specific 


force is; 


l- i 
f.=- (b8)ry (8-15) 
p €é ‘r Ig 
In powerless flight, tg results entirely from lift and drag gas~-dynamic 
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forces: 


fy = - (npcosX + nz sind ) = - nocosY (1 +(L/D)tan¥ ) (8-16) 


For small flight path angles, fy is always negative’. With the left 
side of equation (8-15) negative, the following conditions are required 


of the right hand side: 


lite Yig* Oe: é must be less than +1.0 or must be negative. 


if Vrg7 0 : € must be greater than +1.0. 


If R< O (i.e., descending flight): E must be greater 
than +1.0 he it must be negative. 

i R > 0 (i.e., ascending flight): & must lie between 
O and +1.0. 

Fig. 8.3 shows the resulting flight regimes for a vehicle entering 
the planetary atmosphere. Boundaries separating the various phases of 
the trajectory are shown on Fig. 8.3 in terms of the auxiliary parameter 
é . it is noted that a true Keplerian description of the trajectory 
applies only for E = +1.0; it was assumed in this analysis that the 
Keplerian description of the trajectory is adequate for small departures 
of € from 1.0 (i.e., 0.9<€ < 1.1). 

Fig. 8.3 can best be visualized as a cylindrical graph; i.e., 
the points at é = i{°O should be connected together. The Keplerian 


* For descending flight (i.e., Or negative), it is seen from equation 
(8-16) that f£% is positive only when |(L/D)tan¥| >1.0, a con- 
dition which could arise only for relatively large negative flight 
path angles with a vehicle having a relatively high lift-drag ratio. 
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Phase is a small region in the neighborhood of g = +1.0, The region 
within ~t @.1 of é = 1.0 is the region in which the gas-dynamic terms 
in the equations of motion are negligibly small in comparison to the 
other terms; i.e., the terms involving ny, and n, in equations (6-47) 
and (6-48) are zero for & = 1.0 and are negligible for 0.9L EX eles 
The Gas-Dynamic Phase, on the other hand, is a small region in the 


V_V 
vacanity of ae O. In the band -0.1< € £ +0.1, the 4 


term in equation (6-48) is less than 0.1 of the lift and drag terms and 

can reasonably be ignored. The Intermediate Phase spans all other values of 
é not included in the Gas-Dynamic Phase and the Keplerian Phase. In 

the Intermediate Phase, all vat in the equations of motion mst be 
retained for a reasonably accurate guidance analysis of entry. 

The fact that the parameter e has a sharp behavior during the 
course of an entry mission can be inferred from Fig. 8.3. For this 
purpose, the R and Vz bands shown in this figure are instructive. 
g lies between zero and 1.0 for climbing flight; i.e., v is positive 
in this region. For g negative and for E>4 1.0, the vehicle is in 
descending flight (Y = negative). The vehicle is always slowing down 
( Vxg <0) in ascending flight; if it is slowing down at the same time 
it is descending, then gE must be negative. The speed of the vehicle 
increases only in descending flight and only when E 71.0. At these 
points where V4 = 0, E ist°O , The behavior is sharp because & 
moves from infinity to zero during the time interval between the instant 
that V1g = 0 and R =O. In elliptical flight with small drag forces, 
for example, these two points occur very close to each other in the 


vicinity of perigee and apogee. They would occur simultaneously in 
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pure vacuum flight. In the vicinity of perigee when small drag forces 
exist, Vig goes to zero shortly before R goes to zero. In the 
vicinity of apogee if some drag exists, Vig goes to zero a very short 


time after R=0. 
To illustrate the behavior of é as a consequence of the fore- 
going conditions, it is instructive to consider a couple of simple 
examples. Fig. 8.4 shows a typical braking pass in the planetary 
atmosphere such as occurs in the vicinity of perigee during the un- 
stable orbital profile. When the vehicle is approaching the planet 
at a very high altitude, gas-dynamic forces are negligible and G is 
near 1.0. In this, the Keplerian region, the vehicle is descending 
(R <0) and velocity is increasing (Vag > 0) due to kinetic-potential 
energy trade-off. As drag forces increase, é increases to 1.1 when the 
Keplerian-Intermediate Phase boundary is crossed; velocity is still 
increasing due to energy trade-off, but the rate of increase decreases 
as drag forces become more pronounced. Ultimately in the Intermediate 
Phase, the effects of drag and the effects of energy trade-off cancel; 
at this point Vg = 0 and g passes from +°° to =-co , At some 
later time perigee is crossed (R =0)4 é has increased rapidly from 
-cO to its value of zero at perigeal passage. As the vehicle now 
ascends through the outer reaches of the atmosphere, R is positive and 


‘x¢ is negative; here g continues to increase as atmospheric drag 
forces become less significant until at about E = 0.9, the vehicle is 
again effectively in Keplerian flight. During this braking pass, a 

pulse of energy and a pulse of angular momentum are transferred to the 


planetary atmosphere. 


Fig. 8.5 shows a gliding direct entry profile. Unless skipping 
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motion is induced by the planetary atmosphere”, the radius is always 
decreasing; i.e., R< O. During the Keplerian portion of flight and 
during part of the Intermediate Phase, velocity is increasing due to 
kinetic-potential energy trade-off. e increases from its character- 
istic value of 1.0 in Keplerian flight to +°O as drag forces increase. 


When drag forces reduce Vig to zero, Sl moves from + to - at in- 


ih 
fi Ly. & increases from -coO toward zero as gas-dynamic drag builds 
up with altitude loss; when & = -0O.1, the vehicle enters the Gas- 
Dynamic phase and remains in this phase of flight thereafter. 
Equations (8-4), (8-12) and (8-15) expressed rate of change of 
angular momentum, energy, and horizontal force level as functions of 
Ee ™ In all eases, e appears in the combination (ae). The magnitude 
of this combination was defined in this thesis as the Conservation 
Parameter; this is a convenient quantity for defining the boundary 
conditions between phases. At the Keplerian-Intermediate Phase boun- 
dary, the Conservation Parameter is approximately equal to 0.1; and 
at the Intermediate-Gas-Dynamic phase boundary, the Conservation 
Parameter is approximately equal to 10. The information of Fig. 8.3 
is shown on Fig. 8.6 as a function of 2). From this figure it is 
seen that the vehicle is in ascending flight when (1-8) is positive 
and in descending flight when this quantity is negative. Vig is 
negative in all regimes except -1.0< Gf) < O; i.e., the horizontal 


velocity of the vehicle can increase only in the Keplerian Phase and 


during a portion of the Intermediate Phase. 


* See Chapter 9 for a discussion of skipping flight. 
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8.4 Altitude Description of Boundaries Between Operational Phases 
of the Entry Trajectory 


The horizontal component of specific force was written as a func- 
tion of the Conservation Parameter in equation (8-15) and as a function 
of lift and drag forces in equation (8-16). Eliminating fg from these 


two equations gives: 


(1:§) ='vig (8-17) 


Ny cosd + a sin & = E 


The following relations were given in Chapter 6 for the dimensionless 
quantities required in Equation (8-17) (assuming the exponential atmos- 


pheric model and no atmospheric winds): 





ny = (ae) CEL) py ve enkh (8-18) 
n, = Np (L/D) (8-19) 
r' =v sin¥ (8-20) 
v_.=vecosy + rflcoosy (8-21) 


1g 
Substituting equations (8-18) through (8-21) into equation (8-17) 


gives: 
2 ze a cosy | 3 (z V ie 
This is written in logarithmic form as follows: 


Kah + C, + 2.59 X 10-7 log | = Lege | 


r Vv Ces, 


Ly C 


(8-23) 
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where 





—— R 
= iL CP (SL) “(m)o 
me 8° os ( CGie Se) al 


The following table gives values of quantities on the left side of 


equation (8-23): 






Table 8.1: Conversion Quantities for use in Equation (8-23) 










a 
et] 


eae 


The left side of equation (8-23) expresses dimensionless altitude 








3.87 X 107 





directly for the Earth (for v >0.1). It,is plotted as "altitude 
function" on subsequent graphs. At very low velocities, the term 
on the left side of this equation involving planetary rotation becomes 


significant. For example: 


If: v cos¥ = 0.01 
W = 0° (equatorial trajectory) 


1/r = 1.0 


2.59 X 107? log (= + tices 2 2aalo 


8.7 miles (Earth) 


lle 


At velocities actually existing in the altitude range for which equation 
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(8-23) is useful (i.e., in specifying boundaries of the Intermediate 
Phase), the third term on the left side is negligible. 

In order to illustrate the behavior of equation (8-23), the vehicle 
with (L/D), = 2-0 listed in Table 5.1 was plotted on Fig. 8.7. Mass 
characteristics chosen for this graph was M/S = 1.0 slugs/ft“. The 
effect on the location of the resulting curves for other M/S is 
indicated in the figure. Four separate flight path angles were selected; 
these ranged in value from ¥ =+ 0.01° to +10°. The effects of 
varying the angle of attack from zero to the value corresponding to 
CO ie, is shown for each flight path angle. 

The curves of Fig. 8.7 were plotted for constant flight path 
angles. In general, the trajectory of the lifting entry vehicle will 
be characterized by a changing flight path angle as it passes through 
the Intermediate Phase, hence the true altitude band in this phase may 
be greater or less than that indicated by a constant ¥ curve. For 
example, if the vehicle has a flight path angle of -10° at an altitude 
om 15 Xx 107? on Fig. 8.7, it corresponds to a near Keplerian trajec- 
tory (abscissa near 1.0). If the magnitude of the flight path angle 
rapidly decreases to -0.01°, as in the early phase of a skip, the 
trajectory (at the same altitude) may be described as in the Gas- 
Dynamic Phase (abscissa less than 0.1). These figures show that for a 
given flight path angle, angular momentum is essentially conserved down 
to an altitude that is strongly dependent on the characteristics of 
the vehicle and on the flight path angle. Once the vehicle has descended 
to this altitude, angular momentum rapidly is transferred to the atmos- 
phere. The altitude range over which the abscissa of Fig. 8.7 moves 


from approximately 1.0 to approximately 0.1 is called herein the 
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Intermediate Phase. Below this altitude, the trajectory is described 
as in the Gas-Dynamic Phase. 

At flight path angles with magnitudes less than 10°, the upper 
and lower altitude bounds of the Intermediate Phase is a very weak 
function of lift coefficient. This weak dependence arises from the 
fact that (at ¥ small) the horizontal component of external specific 
force is made up almost entirely of gas-dynamic drag. This weak 
dependence on lift was verified by plotting curves similar to Eq. 8.7 
for many classes of vehicles with and without the lift term included. 
The lift term generally caused such slight displacement of the phase 
boundaries as to be virtu ally unperceptible. Lift terms become 
Significant only with large lift-drag ratios at flight path angles 
greater than five to ten degrees. 

With the weak dependence of the locus of the Intermediate Phase 
boundaries on lift having been established, Fig. 8.8 was plotted to 
show the phase boundaries for various drag parameters, 

The curves of Fig. 8.8 show the Keplerian-Intermediate phase 
boundary and the Intermediate-Gas-dynamic phase boundary for four 
separate flight path angles, Y = + 0.01°, 0.1°, 1°, and 10°. The 
altitude band spanned by the Intermediate Phase at a given drag 
parameter and a given flight path angle is approximately 25 miles 


for Venus, 20 miles for Earth, and 10 miles for Mars. 


2.5 Description of Boundaries Between Operational Phases of the Entry 
Trajectory in Terms of the Horizontal (14) Component of Specific Force 


The horizontal component of the specific forces measured by an 
accelerometer oriented in the 14 direction was given in equation (8-16). 
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Using equations (8-18) and (8-19) this is written: 


fy = . = b° Sp ae ve e7Kh (1 + = tany¥ ) cos ¥ (8-25) 


In order to express fy in Earth g's, the following is used: 


_ G(m)o 


Teg MeGhE in Earth g's (8-26) 


Using equation (8-26), equation (8-25) is written in logarithmic 


form as follows: 


4 
Knot + C, + Cs = 2.59 X 107? og (1 + = tanY) + log (2S F0 , 


+ log C, cos¥ = log = sh) | (8-27) 
2 
Vv 
where C. is given by equation (8-24) and: 


—" G (m)o 
c= = ion log Gone” (8-28) 


The following table lists values of quantities given on the left side 


of equation (8-27). 


Table 8.2: Conversion Quantities for use in Equation (8-27) 





It is noted that the left side of equation (8-27) expresses dimension- 
less altitude directly, in the case of the planet Earth. 
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Equation (8-27) is plotted on Fig. 8.9 for the classes of vehicles 
listed in Table 5.1 with M/S = 1.0 slug/ft® and for a zero-lift vehicle 
with Cys/M = 20. For small flight path angles and moderate lift-drag 
ratios, the curves are essentially independent of Y . The curves are 
drawn for values of Cp and C,; at (L/D) wax for the lifting vehicles. 


Equation (8-15) may be written as follows: 


fg =- =) v~ sin¥ cosy E + ticos (8-29) 


v cos ¥~ 
At the altitudes and velocities which are important at the boundaries 
of the Intermediate Phase, the quantity in square brackets in equation 


(8-29) is approximately equal to 1.0. Therefore: 


fp =. 1-§) Seimeoss t = - (1:8) ae x (8-30) 


5 


Equation (8-29), and its simplified version (8-30), show the 
interesting result that transition from Keplerian flight through the 
Intermediate Phase to Gas-Dynamic flight is independent of the surface 
loading and lift-drag characteristics of the vehicle. The magnitude 
of the horizontal component of specific force in transition depends only 
on velocity and flight path angle. Hence, if flight path angle, velocity, 
and horizontal specific force are measured or computed by the guidance 
system, then the magnitude of the Conservation Parameter can be speci- 
fied; this is an index of the degree to which angular momentum is 
being conserved and is indicative of whether the vehicle is in the 
Keplerian, Intermediate, or Gas-Dynamic Phases of the trajectory. 

Equation (8-30) is plotted in Fig. 8.10 to show the boundaries 
between the various operational regimes as a function of the horizontal 


specific force measured by the accelerometers of the vehicle. This 
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figure shows that the force level changes by two orders of magnitude 

in traversing the Intermediate Phase and that the specific force level 

at transition between phases is a sensitive function of flight path angle. 
The horizontal specific force level cannot give accurate information 

on the relative rate of transfer of angular momentum from the vehicle 

to the planetary atmosphere without accurate specification of flight 

path angle, which is tantamount to requiring precision measurements of 


vertical velocity (altitude rate). 
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Fig. 8.10 : Level of Horizontal Specific Force at Boundaries 
Between Operational Phases as a Function of 
Flight Path Angle. 
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Chapter 9 


APPROXIMATE ANALYTICAL SOLUTION OF GUIDANCE PARAMETERS AND CONSTRAINTS 
FOR THE DIRECT ENTRY PROFILE; RANGE SENSITIVITY TO ERRORS IN CONTROL 


SYSTEM OPERATION 


gam “antroduction 

The investigation described in Chapter 8 demonstrated that the 
trajectory of a vehicle entering a planetary atmosphere from an initial 
point beyond the sensible atmosphere traverses three distinct operational 
regimes; these regimes are called in this thesis, for convenience in 
identification, the Keplerian, Intermediate, and Gas-Dynamic Phases. 

Studies of entry vehicle dynamics in the past have largely ignored 
the existence of the Intermediate Phase as a separate entity. In some 
trajectories, however, the Intermediate Phase may be the most important 
Single phase in specifying the ultimate destination of the vehicle. As 
an example, consider the degenerate orbit of Fig. 9.1. The trajectory 
of Fig. 9.1 consists of a series of "braking passes" through the outer 
reaches of the atmosphere; energy is transferred to the planetary atmos- 
phere during each perigeal passage. During most of each orbit, the 
vehicle obeys Kepler's laws of planetary motion to a high degree of 
accuracy. However, in the vicinity of perigee, the vehicle enters the 


Intermediate Phase of operation momentarily and returns again to the 
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Orbit #1 


Orbit #2 


Orbit #3 


Planet 


Planetary Atmosphere 


Fig. 9.1: The Degenerate Orbit of a Vehicle About 
a Planet Which Possesses an Atmosphere 


el 


Keplerian Phase. As a result of the impulses in the Intermediate Phase 
each successive apogee occurs at a lower altitude and eventual entry is 
assured. If the pilot desires to control the number of orbits remaining, 
he may change the drag characteristics of the yehicle during these 
impulses in the Intermediate Phase to bring about increased or decreased 
total energy transfer per orbit, and therefore shorten or lengthen, as 
desired, the flight time remaining. The fact that the vehicle enters 
the Intermediate Phase during trajectories such as sketched in Fig. 9.1 
makes eventual entry possible without the use of external thrust; other- 
wise, the flight would continue indefinitely. 

The Conservation Parameter | = | was selected as convenient for 
defining boundaries between the Keplerian, Intermediate, and Gas-Dynamic 
Phases of the trajectory. § was defined as the negative of the ratio of 
percentage change in radial distance of the entry vehicle from the 
planet center to the percentage change in horizontal componemt of 
vehicle velocity with respect to inertial coordinates. 

The phases of the entry trajectory were defined in Chapter 8 in 
terms of the Conservation Parameter as follows: 


(1) Keplerian Phase: 


A true Keplerian trajectory exists when & = 1.0. The 
trajectory is near-Keplerian when 


ral < don 


(2) Intermediate Phase: 





(3) Gas-Dynamic Phase: 
- 
g 


Ae 


The phase boundary between the Keplerian and Intermediate Phases was 


defined as ee | ; and between the Intermediate and Gas- 


= 0.1 
Dynamic 7s a = 10.0. It is noted that the Conservation 





Parameter changes by two orders of magnitude in spanning the Intermediate 
Phase. 

Methods of altering the trajectory of the entry vehicle in accor- 
dance with guidance commands may be vastly different in each operational 
regime. Ina pure Keplerian orbit, perturbations may be introduced by 
the use of rocket or ion thrust or by means of external pressure tech- 
niques, such as through streaming a solar radiation sail. Thrust forces 
may be used equally well in the Intermediate and Gas-Dynamic Phases. 
Varying the lift and drag coefficients”, however, has no effect on the 
Keplerian trajectory but may have a striking influence on the trajectory 
in the Intermediate and Gas-Dynamic operational phases. 

The investigation described in subsequent sections of this chapter 
considers each phase of the entry trajectory separately. Solutions 
obtained are continuous or piecewise-continuous within each phase; the 
Conservation Parameter serves as a unifying quantity for matching solu- 
tions at boundaries and for relating the solution of one phase to the 
solutions of other phases. The accuracy of approximations made in each 


= Lift and drag coefficients may be controlled by changing the 
angle of attack or by introducing auxiliary high-lift or high-drag 
devices. Since lift and drag coefficients are coupled through angle 
of attack, auxiliary devices must be used if independent lift and/or 
drag rol is necessary. It is important to note that the coupling 
of lift and drag with angle of attack is different depending on which 
side of the lift curve the vehicle is operating. For example, the 
vehicle may be flown in the "high drag configuration"; i.e., angle of 
attack near 90°. Decreasing the angle of attack reduces drag and 
increases lift in this configuration. In the "low drag configuration" 
(small angles of attack), however, reduction in angle of attack causes 
both lift and drag to decrease. 
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phase are discussed, and the accuracy of the solutions derived are 
compared to numerical solutions obtained from the non-linear equations 
of motion. Limitations of the approximate solutions were established 
as a result of this comparison. : 

The feasibility of using the Conservation Parameter as a switching 
function for the guidance system and as a prediction function for range 
control in degenerate orbits is discussed. In-flight computation of 
the Conservation Parameter from navigational data is shown to be possible. 

The solutions of this chapter are oriented toward the determination 
of guidance quantities. It may be noted that guidance considerations 
of this chapter are relatively unsophisticated, hovever, consisting 
primarily of determination of range capabilities and sensitivity in range 
to control system errors. The hazards of exceeding tolerable accel- 


eration and temperature levels are considered. 


9.2 The Conservation Parameter as a Switching Function: 
Considerable simplification is possible in the design of the 
guidance system if: 
(1) Computations are based on simplified” equations of motion. 
(2) An adequate means is available to shift from one simplified 
guidance mode to another as phase boundaries are crossed. 
The switching parameter should be based on in-flight 
measurements and/or computations. 
At any given time as the entry vehicle passes through the planetary 
atmosphere, the dominant terms in the equations of motion depend upon 


* Provided reasonable accuracy is not sacrificed in the process of 
simplifying the equations. 
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the phase of the trajectory in which the vehicle is operating. The 


parameter - | may be a convenient switching function to shift the 





guidance computer from one set of simplified equations (Keplerian Phase), 
through the region in which few simplifications in the dynamical equa- 
tions of motion are permitted (Intermediate Phase), into the region where 
another set of simplified equations provide an adequate description of 
the trajectory (Gas-Dynamic Phase). 


It was shown in Chapter 8 that: 


e (m)o (1-8) vo. 
fag = - : aro _ = sin} cos / (9-1) 
where fag = tangential component of specific force (in Earth G's) 
measured by accelerometers carried by the vehicle. 

Vv = velocity of the vehicle with respect to coordinates 
rotating with the planet (non-dimensionalized with 
respect to circular satellite velocity at the surface 
of the planet). 

ve = flight path angle; angle between the local horizon and 
v -- positive for climbing flight. 

re = distance from center of planet to vehicle, non- 
dimensionalized with respect to mean planetary radius. 

[(m)o 


a = ratio of mean gravitational acceleration at the surface 
(m)E of the planet 0 to that at the surface of the Earth. 


Equation (9-1) shows that if flight path angle, velocity, altitude and 
tangential specific force are measured or computed by the guidance 
system, it is an easy matter for the system to determine (2:8), which is 
a measure of the degree to which energy and angular momentum are being 
conserved and is indicative of whether the vehicle is in the Keplerian, 
Intermediate, or Gas-Dynamic Phases of the trajectory. Equation (9-1) 


is significant because it demonstrates that the specific force level at 
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transition between phases is independent of the lift-drag characteristics 
of the vehicle, the frontal loading of the vehicle (M/S), and the atmos- 
pheric density variation with altitude. 

The variation of atmospheric density with altitude is not known 
or measurable in advance to a high degree of accuracy, particularly for 
the first-time entry into the atmosphere of a strange planet. Atmos- 
pheric density may vary more than two to one between day and night, may 
vary significantly with latitude, and may change considerably with the 
seasons of the year. The specification of the aerodynamic properties of 
the vehicle may also be subject to error. Determination of aerodynamic 
characteristics at near-orbital velocities is presently not experiment- 
ally feasible; these properties are generally specified on a theoretical 
rather than an experimental basis. 

Equation (9-1) has serious disadvantages as a suitable means for 
computing CA). The direction of the horizon must be known in order to 
measure the horizontal component of specific force. An extremely trouble- 


some factor is sin?’ , which is small for small flight path angles. If ~ 


1-8 
= 


is measured as 0.1°, but is actually 1°, then the computed value of ( 
is off by a factor of ten. Strong dependence on flight path angle is 
the most serious handicap on the practical use of equation (9-1). 

A method for computing the Conservation Parameter that is indepen- 
dent of the sine of flight path angle was derived in Chapter 8 by 


relating it to Energy transfer rates: 


Bi (kin) ¢ | (1-8) ,2 


r= 1 
t aT _2 
E" (pot) g I (9-2) 
In equation (9-2): 
vyé 
ean) = dimensionless kinetic energy per unit mass =a 
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Vr = dimensionless velocity of the vehicle with respect 
to inertial coordinates. 
® (pot) = dimensionless potential energy per unit mass = 
= GR 
Gono ave 
(zero level at infinity). 
R = radius from center of planet to vehicle 
G = local gravitational acceleration at vehicle. 


The prime denotes differentiation with respect to dimensionless 
time T. 

Equation (9-2) has all the advantages of equation (9-1), indepen- 
dence of atmospheric density and vehicle aerodynamic characteristics, 
plus the additional strong advantage of being essentially independent 
of flight path angle. The energy quantities in equation (9-2) depend on 
velocity, altitude, and gravitation; these quantities may be measured 
or computed by the navigation system from data obtained either within 


the vehicle or from external tracking stations.” 


9.3 The Conservation Parameter as a Prediction Function 

One of the problems encountered in the guidance of vehicles under- 
going braking passes through the atmosphere is control of the point at 
which final entry is to be initiated. A vehicle entering the atmosphere 
of a strange planet from an interplanetary transfer ellipse may make 
many braking passes through the atmosphere prior to final entry. A 


C—O ed 


* For first-time entry into a planetary atmosphere, a parent satellite 
from which the entry vehicle is launched appears to be the best external 
source of navigational information in the absence of ground-based 
tracking stations. If there is no parent satellite, then a navigational 
satellite may be left in the original orbit prior to initiating entry. 
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similar trajectory may result if the vehicle is launched from the sur- 
face of a planet for the purpose of orbiting a few turns and re-entering. 
The latter trajectory pattern is planned for the Mercury "man in space" 
program scheduled for operational firings in the early 1960's. 

The most significant characteristic of the degenerate elliptical 
orbit is that the perigeal altitude remains essentially constant” while 
the apogeal altitude drops after each pass through the atmosphere. The 
vehicle transfers energy to the planetary atmosphere in the vicinity of 
perigee during each pass until the total energy level is reduced to a 
point where further orbits may not persist. 

General characteristics of the time-varying nature of altitude, 
total energy, and the Conservation Parameter are shown in Fig. 9.2. 

This figure does not represent quantitatively any particular trajectory, 
but it does contain all the features of a typical trajectory. During 
most of each orbit, the vehicle obeys Kepler's laws of planetary motion 
to a high degree of accuracy. A near pulse of energy is transferred to 
the surrounding atmosphere as perigee is approached. Here, the vehicle 
enters the Intermediate Phase momentarily, then returns to the Keplerian 
Phase. 

After each perigeal passage, the total energy level remains essen- 
tially constant for most of the following orbit at a lower level than 
that of the preceding orbit. The elliptical orbit consequently undergoes 


a "circularization" process. The apogee of the elliptical orbit is 


slowly reduced by drag, primarily exerted near perigee, eventually to 


become a circle. Thereafter, the circle decays spirally as the vehicle 


* See Chapter 10 for a discussion of perigeal altitude decay rates. 
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(a) Altitude vs. Time 
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glides through the Gas-Dynamic Phase of flight to a landing. When 
the total energy level is reduced below that required for a circular 
orbit at perigeal altitude, final entry is inevitable. 

It is interesting to examine the nature of the energy transfer. 
Since each succeeding perigee occurs at approximately the same altitude, 
the potential energy is nearly constant. Therefore, the energy trans- 
ferred at perigee is basically kinetic energy. At the following apogee, 
however, kinetic energy is greater than it was at the preceding apogee 
Since the velocity of the vehicle at apogee increases as the apogeal 
altitude is reduced. Therefore, the energy transferred during the 
perigeal passage results in an increase in kinetic energy and a decrease 
in potential energy at the following apogee when compared to the preceding 
apogee. During the circularization process, the changes in levels of 
energy when viewed at perigee and apogee may be summarized as follows: 

(1) At perigee: Potential energy level remains constant; 
kinetic energy level decreases with each 


impulsive energy transfer. 


(2) At apogee: Potential energy level drops as altitude of 
apogee decreases; kinetic energy level rises. 


The length of time that the vehicle is in the Intermediate Phase 
increases in successive perigeal passes; this is due to the fact that: 


(1) The velocity of the vehicle is slower near perigee during 
each succeeding pass (less kinetic energy). 


(2) The trajectory continually becomes more shallow. 
The energy transferred per orbit becomes greater as the orbit becomes 
more circular. The width of the Intermediate Phase pulses on the graph 


* More precisely, the area enclosed by the pulses. Practically all of 
the energy is transferred in the vicinity of perigee. 
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of Conservation Parameter vs. Time is an indication of the total energy 
transferred per orbit. 

if the total energy level at the first perigeal passage is computed, 
the amount of energy in excess of circular orbital energy (corresponding 
to the perigeal altitude) is easily determined. The magnitude of the 
impulse of energy transferred during the next orbit (i.e., during next 
perigeal passage) shows the rate of energy decay. Energy decays as a 
series of successively largessteps (see Fig. 9.2-b). For given vehicle 
drag characteristics, and for a particular initial energy level, the 
measurement of the magnitude of one energy step uniquely specifies all 
future steps (for constant altitude at perigee and constant drag character- 
istics). Therefore, the number of orbits remaining until the energy 
level has decayed to the circular orbital level may be predicted from 
measurements and computations during the first complete orbit. 

Energy decay per orbital period may serve as a useful method for 
predicting how many orbits remain before final entry; i.e., the number 
of orbits until the onset of the Gas-~Dynamic Phase of flight. With the 
relation between energy transfer rates and the Conservation Parameter of 
equation (9-2), the time-varying characteristics of the parameter oS) 
may serve equally well as a prediction function for determining the 
number of orbits remaining. 

Control of the number of orbits remaining may be implemented either 
by the use of thrust or by changing the drag characteristics of the 
vehicle. If no control action is taken by the pilot, the width of the 
Intermediate Phase pulses observed during successive orbits is an indi- 
cation of how long the flight will last. The flight may be extended or 


shortened only through suitable control action within the limits of 


Zou 


thrust or drag modulation available. 

A framework therefore has been established for the utilization of 
the Conservation Parameter for prediction and control of degenerate 
orbits based on its relation to energy transfer rates (Eq. 9-2) and 
predicated on the fact that knowledge of the energy decay per period 
and the initial energy level uniquely specifies the number of orbits 
remaining. Control of the number of orbits remaining is vested either 
in the engines of the vehicle or in the pilot's facility to vary the 
magnitude of the energy transferred during each braking pass by means 
of drag adjustments. 

The range capability of the vehicle in the Gas-Dynamic Phase is 
small compared to the total range of the orbital phase. If the pilot 
changes the number of orbits remaining by one, through suitable control 
action, he has changed the total range of the vehicle more than the 
distance normally traversed throughout the entire Gas-Dynamic Phase. 
Range corrections in the Gas-Dynamic Phase may be adjudged vernier 
corrections when compared to range corrections in terms of changing 
the number of orbits. 

The function of the guidance system in a vehicle which is under- 
going the circularization process such as sketched in Fig. 9.1 may be 
summarized as follows: 

(1) To determine in advance a suitable geographic point for 
entry; i.e., the point for initiating the Gas-Dynamic 
Paase or flight. 


(2) To adjust the orbital characteristics of the vehicle in 
order tO hit this entry poimer. 


(3) To make range corrections in the Gas-Dynamic Phase to 


correct for position errors in the initial entry point 
and to correct for range errors which may arise as a 


EOL 


result of perturbations from the nominal trajectory in 
the Gas-Dynamic Phase. 


The geographic location of the entry point is determined by solving 
the trajectory backwards from the geographic location of the predetermined 
landing site. The maximum and minimum range capabilities of the vehicle 
in the Gas-Dynamic Phase may be determined from the approximate range 
expressions derived in section 9.6. The entry point is selected as a 
mean position between the limits imposed by backing off from the landing 
point the maximum and minimum ranges obtainable. If the actual entry 
point lies within the corridor thus prescribed, adequate control is 
available to the pilot to prevent overshooting or undershooting his 


destination, 


9.4 The Keplerian Phase 
Any method for guiding the entry vehicle to a selected geographical 


point on the surface of the planet must necessarily involve a pertur- 
bation of its original orbit. Because of the tremendous energy possessed 
by the vehicle in orbit, it is difficult to perturb the orbit greatly 
Without investing a large portion of the total mass of the vehicle in 
rocket fuel. 

The vacuum phase of the trajectory that follows the initiation of 
entry from a satellite orbit beyond the sensible atmosphere is a part 
of the general problem of transfer between orbits that has received 
considerable attention in the literature. The special nature of the 
entry problem, however, requires solutions that are not available in 
general treatises on the subject of transfers between orbits. 

Excellent examples of published papers concerning investigations of 


(49). 


transfer maneuvers between orbits are the works of ieitmann a Lawden 
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and pect 20), Hohmann and Lawden have shown that the optimum transfer 
maneuver between elliptical orbits involves the application of impul- 

Sive forces at departure from one orbit and at arrival in the other 

orbit. Battin has shown that several important features of the transfer 
problem are basically three-dimensional in nature and that two-dimen- 
Sional models are inadequate. Very little has been published concerning 
guidance requirements for the special problem of transfer from a satellite 
orbit to a particular landing point on the surface of the planet. 

The development carried out herein is unique in that the transfer 
ellipse is solved in terms of geometric quantities at the trajectory 
modification point”. These quantities are: 

(1) Inertial flight path angle (see Figs. 9.3 and 9.4) 
(2) Altitude 
(2) Meloedsby 
The solution is valid for entry from either circular or elliptical orbits. 

In addition to the general solutions described above, special 

solutions were obtained for entry from circular orbits in terms of: 
(1) Velocity Impulse 
(2) Engine Gimbal Angle 
(3) Circular Orbital Velocity. 

It should be noted that the relation between | inertial il ohe 
path angle, altitude, and velocity } at the trajectory modification 
point and | velocity impulse, engine gimbal angle, and orbital velocity} 
are uniquely defined only if the entry vehicle is launched from a 


and ap cc: = —_ ana = i=] an ea aa =D ap oz t= =] -—_™ 


* The trajectory modification point is that point in the original orbit 
at which entry is initiated through generation of thrust forces, see 
Pig. 9. 
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circular orbit. For transfer from an elliptical orbit, the relation 
between these two sets of quantities depends on the particular point in 
the original elliptical orbit where perturbations are introduced. The 
first set of quantities are more useful for generalized study of transfer 
from elliptical orbits, the second set for specialized study of transfer 
trem Circular orbits. 

Results are iehecnten in dimensionless form for general application 
to any planet. A generalized range expression is developed in terms of 
the foregoing quantities, and range sensitivity to errors in these 
quantities is discussed. 

Various methods are available for perturbing the satellite orbit 
to bring about controlled entry: 

(1) Impulsive application of forces at departure point by 
means of chemical rockets or other high thrust pro- 
pulsive systems. 

(2) Continuous application of low thrust (e.g., ion rockets). 

(3) Multiple impulses at intermediate* thrust levels. 


(4) Drag modulation, if portions of the orbit pass through 
the planetary atmosphere. 


The analysis of this section is restricted to the first method. The 
analysis is simplified to the two-dimensional trajectory near a spherical 
planet. 

The Keplerian Phase of the trajectory was defined previously as 


that portion of the trajectory where 7 | < O.1. The equations of 





motion developed in Section 6.4 may be written for the Keplerian Phase 


aS: 
x = 0 (9-3) 
N Yr 


eae Fe Fe FSP Pe FSF SF. Fe EF Fe SF FP el ee] SlU 


* Thrust levels between those obtainable with chemical rockets and 
with continuous ion propulsion. 
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hn? "rT! =e =av 


— (9-H) 


ll 
< 
ie 
1) 
i 


Wat r ~ 2 (9-5) 
Vg! = == (vg = Zz Vig) ‘ (9-6) 


Fig. 9.3 describes the geometry of the entry trajectory. Quantities 
at the trajectory modification point are denoted by the subscript "m", 
The transfer ellipse is solved in Derivation Summary (9-1). A brief 
summary of the solutions obtained are tabulated in sections 9.4.1, 
9.4.2, 9.4.3, and 9.4.4 following. 

9.4.1 Description of the Intry Transfer Ellipse in Terms of Geometric 


Quantities at the Trajectory Modification Point (Applicable to 


Entry from Elliptical and Circular Orbits). 


a) Eccentricity: 
e= |i - ae cos’ ¥ (2 - ve cm) 2 (9-7) 
Im m Im Im m 


b) Dimensionless semi-major axis: 


r 
a = e-em) (9-8) 


c) Dimensionless semi-minor axis: 

by = ay ft -E? (9-9) 
d) Dimensionless latus rectum: 

yy eget €~) (9-10) 
e) Dimensionless area of ellipse: 


Area = T[ayby (9-11) 
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Apogee of Transfer Ellipse 





Trajectory Modification 
Pomnt 


=o 


Iw 


Apocenter 


Pericenter 


VV “Ss Original Satellite 
Orbit 


— Entry Trajectory (Transfer Ellipse) 
Perigee of Transfer Ellipse 


Ww, = dimensionless circular orbital velocity prior to modification 
dar = Gimensionless velocity impulse imparted 
Oy, = true anamoly of trajectory modification point 


Wz, = Gimensionless velocity of vehicle with respect to inertial 
coordinates immediately after modification 


%,, = inertial flight path angle immediately after modification 
Ae = engine deflection angle 
ky = dimensionless distance traversed from modification point. 


ANGULAR RELATIONS 


(9.3-1) sin(-77,) 
(9.32) cos(- 9y,,) 


(§ v)sin A,/vq), 
Ve-(3 v) cos Ae | [Vom 


(} w)sin A,/ [vq-(8 v)cos Ae | 


Ut 


tl 


(9.3-3) tan(-7 7,) 


VELOGLIY RELATIONS 
2 2 1/2 
(9.3-4) [(3 v) + Ve 200 ec ose 


GR iE 
(9.3-5) ve = Ginyo®imo = J r 


Fig. 9.3: Geometry of the Transfer Ellipse; Definitions of 
Quantities and Geometric Identities. 
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f) Dimensionless period of orbit: 


" 3/2 
Te ill (9-12) 
g) Dimensionless total energy: . 
E = -1/2 o1 
h) Dimensionless angular momentum: 
Pp = F,V7,008 ae (9-14) 
i) Dimensionless radius at perigee: 
a ay-€ ) (9-15) 
j) Dimensionless velocity at perigee: 
lane 
“Igy a eile (9-16) 
k) Dimensionless radius at apogee: 
rq — Syitre ) (9-17) 
1) Dimensionless velocity at apogee: 
I = 1=€ (9-18) 
a aa 


9.4.2 Guidance Quantities in Terms of Geometric Quantities at the 


Trajectory Modification Point (Applicable to Entry from 
Elliptical and Circular Orbits) 
a) Time of flight (T = 0) 
3 1 /(€)tan 9/2)\ ~ 
s €é sin oe n- tan 9/2 
ai ey + ( ) (9-19) 
a“ 1L-€ Go +€cos 9 - 1- mass oT = ze 9 
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b) Dimensionless velocity: 


_ |2_ (2- “In‘'n) 
mle 


c) Inertial flight path angle: 
Ra = arc cos “n “In Posi, 
I yay in 


Using equation (9-20), this may be written: 


Caen Ke 
1. = are cos! im_cos " Im 
I r _ Za, i 
rm 
d) Range: 
Ky = 9 - On 
where 


2 

cos 9 =(1k) ee - 1) 
z 

cos O =(1/.) So =a) 


Range may be expressed in closed form as: 


cos Ay opel 


m 





where: 


A 
2 


. (0? - r)(p* - me) + [€?r? ~ (pone r) 
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as 


(9-20) 


(9-21) 


(9-22) 


(9-23) 


(9-24) 


(9-25) 


(9-26) 


- (p* - 7,)°| ‘ 


(9-27) 


9.4.3: Description of Transfer Ellipse in Terms of Velocity Impulse 
and Engine Gimbal Angle. Applicable to Entry from Circular 
Orbits Only. 


The variables used in the following equations are defined in 
Fig. 9.3. Equations (9.3-1) through (9.3-5) given in Fig. 9.3 give 
various relations for vehicle velocity and inertial flight path angle 
at the trajectory modification point in terms of velocity impulse and 
engine gimbal angle. 

The following quantity is defined: 
velocity impulse non-dimensionalized with respect (9-28) 


to circular orbital velocity at the trajectory 
modification point. 


Sv. 


tt 
PS 
| 


a) Eccentricity: 


a 
2 
C= {2 alee) Vv. cos in) (2 - ov, (dv. - 2 cos +2] (9-29) 
b) Dimensionless semi-major axis: 
_ i 
i oeve = § Vv, ()v, - 2 cos Ag)) (9-30) 


All other quantities listed in Section 9.4.1 are the same except the 
following: 
h) Dimensionless angular momentum: 


al - §v.cos A.) 


Ee 


p (9-31) 


9.4.4: Guidance Quantities in Terms of Velocity Impulse and Engine 
Gimbal Angle, Aplicable to Entry from Circular Orbits Only. 


a) Time of Flight: 


Same as Equation (9-19) 
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b) Dimensionless Velocity: 


Ua E re ve [2 = dv. Cove = 2 cos 4.))] i (9-32) 


c) Flight Path Angle: 


1 = arc cos (1 - d veces = (9-33) 
Vi ovyr 
f 
or 
6 Gli § vcos A.) i 
= a O-)~C pe pS ee 2 = 
a. CC (2v,.“r SY eta ae Ov. (Sv, = 2 cos A. | )| ae 


d) Range: . 
Equations (9-23) through (9-25) apply. An alternate expression for 


range is; 


Nh 


Vv 
cos X _ { 


N € (9-35) 


se 


where: 
1 ac 2 3 
ol 2 2 2 lee 
was {oF - 9G - % 2) + [6x - (p - r)?| [, - ( - 45)7] 
c V V 
G cj 
Fig. 9.4 shows the geometric relation between flight path angle, 0 , 


and inertial flight path angle, 4%. Equations (9.4-1) through (9.48) 


T° 
given as part of this figure are useful for relating velocities and 
angles. The solutions for the Intermediate and Gas Dynamic Phases given 
in later sections of this chapter are generally written in terms of 1 
because this angle is more convenient when drag and lift terms are 
important (since these quantities involve velocities with respect to the 


atmosphere of the planet). Fig. 9.4 and the equations listed therein 


enable conversion of angles at phase boundaries such that the most 
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Vr = Vz Cos % =v cos V+rtLcos ¥ 
Weh’=z sin TY sn 7 

p=Yrvyz cos F 

sin am ie 

v= -2arz cos Hr tL cos Vtr* cos *¥ 


ue tea cos vr 1 cos Yt 7’ cos? p 


f= Vrwosh-vfcost 
Ww 


COS 


S 
aoen = wes tr LL CO Ww 
Vy 


(94-1) 


(9item) 


(9.4-3) 


(9 4-1) 


(9.4-5) 


(9.46) 


(9 4-7) 


(9.48) 


Fig. 9.4: Geometric Relations Between Flight Path Angle and 


Inertial Flight Path Angle. 
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convenient definition of flight path angle may be used in the analytical 


study of each phase. 


9.4.5: Optimum Engine Gimbal Angle: 

A thorough graphical representétion of equations (9-7) through 
(9-35) was not undertaken in this thesis, An example of a solution to 
the circular orbital equations (given in subsections 9.4.3 and 9.4.4) 
in order to determine optimum engine gimbal angle is presented graphically 
in Figures 9.5 and 9.6. The circular orbital altitude was chosen as 
nn = .07575, corresponding to an altitude of 300 miles above the Earth. 
The final altitude was selected ta be hp = 0, corresponding to impact 
with the planet. No atmospheric effects were considered. 

The following velocity impulses were chosen: 

Ne = 0.0385; 0.0772; 0.1156; and 0.1542 
These correspond to velocity impulses of 1000, 2000, 3000, and 4000 ft/sec, 
respectively, for the planet Earth. | 

Figs. 9.5 and 9.6 show that there is one value of engine deflection 
angle for which range is minimized with a given magnitude of velocity 
impulse.” At this engine deflection angle, any errors in alignment of 
the engine gimbals result in minimum range error. For example, if 
5 y = .0385, the optimum engine angle is A, = +340: anos = 1542, 
optimum engine angle is A, = +60°. 

It is interesting to note in Fig. 9.5 that at velocity impulses of 
large magnitudes, range is relatively insensitive to errors in engine 


gimbal angle from the optimum angle. The range sensitivity to errors in 


a anew ws wvmame a sete ef: F2-F Heo CeO 2 COO  kSSlllUWCCOCUllUllCC 


* Labeled "Optimum Engine Deflection Angle" on the graphs. 
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Fig. 9.5 : Vacuum Ground Range vs Engine Deflection Angle. 
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engine gimbal angle increases as the magnitude of the velocity impulse 
Ss reduced. It might be concluded that large thrust capability for 
extended periods (i.e., big engines with large quantities of fuel) is 
the cure-all for improving impact accuracy of the entry vehicle. A 
compromise is generally dictated in. fuel and engine weight allowances, 
however, by payload requirements. 

It may be observed that with a given velocity impulse and range, 
the engine deflection angle is a two-valued function (except at 


A, ). A given range may be obtained with either of two values of 


optimum 
engine deflection angle. 


It is significant to note in Fig. 9.5 that the range interval 
between curves decreases as Sy is increased** Therefore, the magnitude 
of reduction in range decreases for a given step increase in } y at high 
impulsive velocity levels. 

Graphical determination of range sensitivity to errors in either 
engine deflection angle or velocity impulse exemplified by Figs. 9.5 and 
9.6 is cumbersome and inefficient. A particular set of initial con- 
ditions must be assumed; then the trajectory must be solved repeatedly 
uSing systematically chosen values of engine gimbal angle and velocity 
impulse. A more effective analytical method of obtaining the same 
information was derived in Derivation Summary 9.1. The results of this 


derivation are discussed briefly in the next two sections. 


* That is, the curves are more pointed in the vicinity of A, 


for small § ,. 


optimum 


** On fig. 9.5, the curves are more closely spaced at high values 
of ae. 
Vv 
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9.4.6: Analytical Determination of Range Sensitivity to Errors in 
Operation of Propulsive System: Range Sensitivity to Errors 
in Velocity and Inertial Flight Path Angle (Applicable to Ent 
from Elliptical and Circular Reconnaissance Orbits). 
A single equation relating altitude h, altitude at which entry is 
initiated h,, velocity and inertial flight path angle at the trajectory 
modification point, and range was derived in Derivation Summary 9.1 as 


follows (note: r<=1 +h): 


en 2 = 
= Ba C03" Tig 1 Pm Vin cos © yon COS (Xt) cos XN (9-36) 





This equation was written in other forms in equations (28) and (29) 

of Derivation Summary 9.1. Implicit partial differentiation of equation 
(9-36) yields equations which may be used to evaluate range sensitivity 
to errors in either altitude, velocity, or flight path angle. 


Equation (9-37) gives range sensitivity to variations in: 


wn 
1 


Xn ein VI, cos* oes [= ~cos(Xy + Cen)) 





Ov, [sin i an Vz COS* Sn (Kut % mn) al 


Equation (9-38) gives range sensitivity to variations in Ge 


Y 
oXy Viton vaeaecOS Mit {2 sin tt [cos (Xy+ Cade =| + cos Tn sin Oy) 
Ot, [sin Saar Vo cos’ %,, SIN K%,, | 


(9-38) 
9.4.7: Analytical Determination of Range Sensitivity to Errors in 
Operation of Propulsive System; Range Sensitivity to Errors 
in Velocity Impulse and Engine Gimbal Angle (Applicable to 
Entry from Circular Reconnaissance Orbits Only). 
A single equation relating altitude, altitude of the circular orbit, 


velocity impulse, engine gimbal angle, and range was derived in the 


eo 


following forms by the methods of Derivation Summary 9.1. 


(Sa) COS ‘ae cos Xyt sin Xy tan A .|- da, cose a(= —cos Xy)+sin Xn tan Ae|t 2-1=0 


(9-39) 


or 


Sar = {else Kw) + sin Xy tan Ae| ~ [sin® Me tan’ Aet 4 (1-cos Xy(—2— 


%, 
cos Ay t tan Ae sin xy} 
Z2cos A .(~ cos Xyt+ sin Xy tan Ae) 


(9-40) 
Range sensitivity to retro-rocket system operation was determined by 


implicity differentiation of equations (9-39) or (9-40): 


spas sin X { =~ Yr 
OX y ie (sa) in Ae (cos ie pe SIN Xn tan Ke) + cos ee cos*Ae sin a a 7 | ) 
dAe eae Sie eM Re lemm ercos Xa ton ‘ (9-42) 


| 
dXy 2,cosAe (=~ cos ky T sim Xy tan he) oa (z (= COS Xt sin Xv tan Ae| 
Spy = A 
o Sng) ie sim wa Gos Xn tan ANe)-Snr% sin (XytAe) 


(9-42) 
With the aid of equations (9-36) through (9-42), a description of range 


sensitivity to errors in retro-rocket system operation is available 


for analysis of entry into any planetary atmosphere from elliptical or 


circular reconnaissance orbits, 


9.5: The Intermediate Phase 


The Intermediate Phase presents the most complex computing phase 


for the guidance system. In this phase accelerations due to gas-dynamic 
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terms are of comparable magnitude with other terms in the dynamical 
equations of motion. Some of the simplifying assumptions made by most 
authors* ‘analyzing entry trajectories are not permissible in the Inter- 
mediate Phase. A partial list of typical assumptions and their accuracy 
with respect to the Intermediate Phase are summarized below: 

(1) Lift and Drag coefficients are independent of Mach Number and 
Reynolds Number. This assumption is generally accurate at the 
high Mach mmbers of the Intermediate Phase. 

{2) Gravitational acceleration is constant. This is clearly a poor 
assumption in the Keplerian Phase, and a questionable assump- 
tion in the Intermediate Phase. For example, the gravitational 
acceleration of the Earth decreases about 1% for each 20 miles 
increase in altitude. It was shown in Chapter 8 that, for 
constant 7, the Intermediate Phase spans an altitude band of 
20 miles for the Earth's atmosphere. If a value of gravitation 
acceleration is used that corresponds to the average altitude 
of the Intermediate Phase, then variations of gravitational 
acceleration of less than 1% from this local average value may 
be expected in the Intermediate Phase. 

(3) An isothermal atmosphere is assumed; therefore, atmospheric 
density decays exponentially with altitude. This assumption 
is reasonably accurate for Earth below altitudes of 80 miles 
(see Fig. E.2); the accuracy of this assumption for other 


planets cannot be predicted at this time. It was shown in 


ga _., = se =F Se 2  & = =] 2 8&2. = = 


(2) (31) | 


* See, for example, Chapman , Eggers, Allen and Neice and 


Allen and Becera >). 
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Chapter 8 that the onset of the Intermediate Phase may be 
higher than 80 miles. In general, the altitude of the Kepler- 
ian-Intermediate Phase boundary increases as: 

(1) Drag coefficient is increased; 

(2) Surface loading (M/S) is decreased; 

(3) The magnitude of flight path angle is decreased. 
Since the altitude of the Intermediate Phase may be above the 
region where the exponential approximation of the atmosphere 
is accurate, either a power approximation or a variable decay 
parameter k may be more accurate representation for analytical 
studies of this Phase. 

(4) Planetary rotation is neglected; therefore, Coriolis forces are 
not included in the equations of motion (even though these 
equations are written in a coordinate system that rotates with 
respect to inertial space). This assumption is reasonably 
accurate in the Intermediate and Gas-Dynamic Phases for 
terrestrial planets. 

(5) (1+ L/D tan7) ~ 1.0. This assumption is accurate in all 
phases for moderate lift-drag ratios and normal flight path 
angles (except in steep terminal phases near the end of the 
Gas-Dynamic Phase). 

(6) Quantities in the equations of motion involving = 
are neglected. This is equivalent to assuming ra | > 10 


Gites. 





&| << 1.0). It was shown in Chapter 8 that this 


assumption is accurate only in the Gas-Dynamic Phase. This 


* Or terms equivalent to this quantity. 
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assumption cannot be made in the Intermediate and Keplerian 
Phases. Solutions to the simplified equations of motion in 
the Keplerian and Gas-Dynamic Phases are discussed in other 
sections of this chapter. Few simplifying approximations are 
warranted in the Intermediate Phase, however, and closed form 
solutions were not obtained for this phase. 

In powerless flight”, the two-dimensional equations of 
motion may be written in terms of velocity, altitude, and the 


Conservation Parameter as follows: 


x - £ : (9-43) 
N 
h' =ric ly 
— , (gate) 
a (1 Yr “I¢ (i-8) (1- 4/p) vz) 
Vv. = -=- == -4 
r =» oe vgr  g (FE) =) (9-45) 
Vv 
i (1+8) 
Vg =e Ivg - g vag | (9-46) 
where 
Vig = %¢ +rfcos P (9-47) 
Equation (9-45) may be simplified slightly by assumption (5) discussed 
previously: 
Ly ail L 
(1+ “/p up = (1+ “/p tanY) = 1.0 


The Intermediate Phase may be traversed impulsively by using certain 


* Thrust 


terms are zero in powerless flight, i.e., Vy = 0. 
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trajectory control techniques at the onset of this phase; i.e., at the 


point where pe = 0.1. Examples of such techniques are: 


(1) 


(2) 


Application of thrust: 
When the phase boundary between the Keplerian and Inter- 


mediate Phases is approached, the Conservation Parameter rises 
from near zero to 0.1. At the point where this parameter 
reaches 0.1, impulsive retro-thrust may be applied in such a 


manner as to cause | to step to 10 or greater, corresponding 





to flight in the Gas-Dynamic Regime. 
Use of Auxiliary High-Drag Device” 

Fig. 8.8 showed upper and lower phase boundaries of the 
Intermediate Phase as a function of drag coefficient. ‘It may 
be observed from this figure that a step change in drag 
coefficient of approximately two orders of magnitude is re- 
quired to move the operating phase of the vehicle from the 
Keplerian-Intermediate Phase boundary to the Intermediate- 
Gas-Dynamic Phase boundary (with constant altitude). This step 
drag increase may be obtained at the onset of the Intermediate 
Phase, for example, by streaming a drogue chute or umbrella. 
If the horizontal component of external specific force steps 
two orders of magnitude as the drag device is extended, then 


the Intermediate Phase is traversed during the step drag increase. 


Certain lift-drag programs** may be performed in the Intermediate 


_-_ «- - ape |. 2 >] lUWDDeSlClleSllWeeSSlUeSlUlmltCOF 


Such as a drogue chute or umbrella. 


Programmed lift and drag, such as discussed in this section, are 
not easily realized in practice. 
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Phase to induce a particular trajectory pattern. One example of such a 
lift-drag program is that of holding constant flight path angle through- 
out the Intermediate Phase. 

The required lift program to hold constant flight path angle, J , 
may be determined by setting {= O in equation (6-41). The drag 
program is determined by integrating from equations (6-40) and (6-43) 
with 7 constant. It may be seen that the flight path angle remains 


constant if lift and drag are programmed as follows: 
-khi 
e 


7 cot 7 (1 - BENG -~ eh ) 


iE 
eg. eo 
al 


(9-48) 


Peo, (1 - v~) 
CL = (OSDIR ado) ,2g-eh Jatt) 
2M/S 


_k sin V1n = 
"(PSE F (mo) (j-h_g- ky dale 


2M/S 


It is noted that the exponential atmospheric model was assumed and that 
terms involving planetary rotation were dropped because they are higher 
order effects in comparison to the quantities retained. 

By switching to the above lift-drag program as ral risesetonOr 
and maintaining it until 3 = 10, the Intermediate Phase is traversed 
at constant flight path angle. Closed form expressions for range and 
other quantities developed in section 9.6 for the "Ballistic Trajectory" 


apply in the Intermediate Phase if this lift-drag program is maintained. 


9.6 The Gas-Dynamic Phase: 


The final portion of any entry mission to the surface of a planet 


through its atmosphere necessarily involves flight in the Gas-Dynamic 
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Phase, The analysis of motion of vehicles in this operational regime 
has been a substantial percentage of man's total scientific effort during 
the past half-century and resulted in the evolution of an entirely new 
branch of science and engineering, viz. "aeronautics". With each advance 
in propulsive system design, the increased velocity and altitude capabil- 
ities of the vehicle has brought on problems of increasing complexity. 
The German V-2 program was the first substantial step up the ladder 
toward flight beyond the sensible atmosphere at hypervelocities. This 
beginning has been greatly extended in the IRBM and ICBM programs in the 
United States and the Soviet Union. 
The character of the atmospheric phase of flight at near orbital 
velocities is, in many ways, far removed from that of flight at sonic 
or subsonic velocities. The vehicle possesses tremendous energy and 
angular momentum at near-orbital velocities. Most of this energy must 
be transferred to the atmospheric envelope surrounding the planet prior 
to initiating a safe landing. Two of the primary dangers encountered 
during the process of transferring this energy are manifested as heating 
of the vehicle and deceleration of the vehicle and its occupants. It 
was not until the feasibility of space activities was clearly demon- 
strated in the post World War II period that analytical investigation of 
hypersonic flight was carried out on a scale more extensive than the 
occasional publication of a paper touching the subject™. Since 1955, 
the interest in this problem has grown rapidly, as evidenced by the 
steadily increasing frequency of published work. Among the papers 


ee ee ee > ~~ - ee — ef 


* One of the earliest papers of significance is that of Sanger?) . 
published in 1933. 
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published in recent years, some of the best are the works of Sanger and 
pepe >? (1944); Eggers, Allen, and Neice (31) (1957); Allen and iepeme so” 
(1957); ea (1957); and Chapman‘15) (1958). Chapman considers 
atmospheric entry into not only Earth tut also other planets of the solar 
system; the other authors cited above restrict their analysis to the Earth 
system. Much of the results presented in the following pages are based 

on conclusions of the foregoing papers. The results that follow, however, 
go beyond the scope of the above treatises to the derivation of guidance 
quantities in terms of certain independent variables (such as time and 
atmospheric density ratio) not previously considered. 

Three basic trajectory patterns are possible in the Gas-Dynamic 
Phase. The trajectory profile which is traced out on entry depends 
strongly on the initial flight conditions and on the lift-drag character- 
istics of the vehicle. These trajectories are classified as follows: 


(1) Ballistic trajectory: Non-lifting vehicle with large initial 
flight path angle. 


(2) Glide trajectory: Lifting vehicle with zero initial 
flight path angle. 


(3) Skip trajectory: Lifting vehicle with finite initial 
flight path angle. 


Approximate analytical solutions were derived in this thesis for 
each of these classes of trajectories. Fig. 9.7 shows the regions in 
which these approximate analytical solutions are an accurate representa- 
tion of the trajectory when compared to numerical solutions obtained 
from the non-linear dynamical equations of motion. 


Chapman‘~?? 


discussed trajectories in the Gas-Dynamic Phase by 
transforming the equations of motion to a single, ordinary, non-linear 


differential equation of second order. A transformation similar to 
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an Phase 


Intermediate Phase 


Gas-Dynamic Phase 





Y= constant = a 


Initial Conditions: Zero lift vehicle with steep flight 
path ; or programmed lift-drag to 
maintain constant flight path angle. 


Limitations of Approximate Solution: 


Approximate solution is accurate 
(1) Throughout the velocity spectrum for 


texas (this corresponds to entry 
at \gal> 5° forwarth). 


(2) For flight path angles less than listed in 


(1) above, approximate solution is accurate 
down to v equal to about 0.8. 


(a) BALLISTIC TRAJECTORY 


Fig. 9.7: Regions Where Approximate Solutions of This Thesis are 
Accurate. 
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a" Ula Phase 


Intermediate Phase 





Gas-Dynamic Phase 


Initial Conditions: Lifting vehicle with zero initial flight 
path angle. 


Limitations of Approximate Solution: 


Approximate solution is accurate: 
(1) For v> 0.2 if vehicle characteristics are: 


1/2 
k L/D > 30 (i.e, (L/D) p14, > 1-0) 
(2) For v >.0.5 if vehicle characteristics are: 
o> fa - w/D > 15) (i vespe dec ey e S705) 
ae ee Barth ~ ° 
(3) For v > 0.7 if vehicle characteristics are: 


15>” L/D> 7.5 (dees, 045 > (L/D) pang p> +25) 


(b) GLIDE TRAJECTORY 


Fig. 9.7: Regions where approximate solutions of this thesis are accurate. 
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_— | Keplerian Phase | 







Intermediate Phase 


Gas-Dynamic 
Phase 





Initial Conditions: Lifting vehicle with finite initial flight 
path angle. 


Limitations of approximate solution: 
Approximate solution is accurate if 2k a L/D| > 1.0 


(i.e., 93 L/D| parth > 5-6 X 107") 


Approximate solution can be applied only to a single skip. 
Subsequent skips may be analyzed by considering each skip 
individually with a new set of initial conditions. 


(c) SKIP TRAJECTORY 


Fig. 9.7: Regions where approximate solutions of this thesis are accurate. 
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Chapman's Z-transformation may be applied to equations (9-45) and (9-46) 


with the following assumptions: 


1+ (L/D) tan®? ¥ 1.0 
ai 

€«1.0 

r = 1.0 


The equivalent Z-transformation is: 


7 rg aw, cos 7 sint (9-51) 
E k?y k? & 


The single equation derived by performing this transformation is: 








d*Z aie (GZ wueos “74 RL 
co, — i 28 a z — 2 3 MED 
P dag Es Bs (| 18) Zar k* cos % (9-52) 
+ 7 ————__— | 
Vertical Vertical Component Gravity minus Lift Force 
Acceleration of drag Centrifugal Force 


Alternate forms of expressing the Z-function are as follows: 


x 
cos 
Tm a a — ig (9-53) 


k*y k?y 
Chapman's numerical solutions are compared to the approximate 
solutions of the dynamical equations of motion derived in this thesis 
in Figs. 9.8, 9.9, and 9.10 for the ballistic, glide, and skip trajec- 
tories respectively. The abscissa on each of these graphs is the 
horizontal component of dimensionless velocity with respect to the 


atmosphere; The ordinate is the horizontal component of specific force 
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Horizontal Specific Force 


% 


in Equivalent Surface Gs of the Planet 


20 


Fig. 9.8 : 





(Avi earth® =o 


v; in radians 


Key. 


. -—— Chapman's 
Numerical Solution 


iO. \e 





Approximate 


Solution of this Thesis 


Note : 


\ 
See Fig 6la) of Ref I5 \ 


0.2 0.4 0.6 0.8 1.0 


Dimensionless Horizontal Velocity a¢ 


Comparison of Avproximate "Ballistic" Solution to Machine 
Computed Numerical Solutions. 
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Horizontal Specific Force 
m Equivalent Surface G’s of the Planet Times Lift-Drag Ratio 


1.0 L/D fy es: 


L/ De aeth= |.0 —-—-— Chaperamie 


Numerical Solution 


Ye = 
on k* L/pD=30 “= Op GeNimate 


y Solution of this Thesis 
VY - 
0.5 \ RY ae DAIS 


© 
oo 
\ 


© oO 
& ~) 
~\ © 

Ee) 

On 

a a 
Or 


be 20 
/ ; \ 
| 
0.5 
\ 
0.4 \ ' 
OS \ 
%,=0 

VG: = 1,0 \ 

G- “| 
\ 
Note ; 
See” rt 4(b) of Ref. 15 
\ 
O 0.2 0.4 0.6 0.8 1.0 2 


Dimensionless Horizontal Velocity v¢ 


Fig. 9.9 : Comparison of Approximate "Glide" Solution to Machine 
Computed Numerical Solutions. 
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Horizontal Specific Force in Equivalent 


Surface G's of the Planet 


0.6 


O25 


0.4 


Oss 


Cz 


0.1 


0.2 
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'| Note: \ 


| / See Fig.G (e) of Ref (5 ' 
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| 
/ Key \ 


\ / co 5 aaa am Chapmans 
Numerical Solution \ 
/ ‘ 
\ ———= Ao proximate Solution 
\ i of this Thesis \ 
a \ 
\ 
0392 094 O36 098 1.0 


vg Dimensionless Horizontal Velocity 


Fig. 9.10 : Comparison of Approximate "Skip" Solution 


to Machine Computed Numerical Solutions. 
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in surface g's of the planet concerned.” 


9.6.1 Solution of the Ballistic Trajectory 
The ballistic trajectory is characterized geometrically by a 


constant flight path angle. Prerequisites for this trajectory may be 
summarized as follows: 

(1) Zero-Lift Vehicles: 

If entry is made along such a steep path that the sum 
of vertical components of lift, gravity, and centrifugal forces 
are small in comparison to the vertical component of drag force. 

(2) Programmed Lift and Drag:. 

If lift and drag are programmed in such a manner as to 
maintain constant flight path angle. (See equations (9-48) 
through (9-50). 

The limitations of the approximate solution developed in Derivation 
Summary (9-2) when compared to more exact computer solutions of the 
non-linear system of equations are briefly: 

(1) The solution derived here is an excellent approximation for 

Jk 7 |> 2.5 throughout the velocity spectrum. This corres- 
ponds to entry into the Earth's atmosphere at angles equal to 
or greater than 5°. 

(2) The approximate solution is accurate for flight path angles 
less than | (x fv |= 2.5 over the limited velocity spectrum 
vw> 0.8. Therefore, the solution is accurate for all flight 
path angles down to velocities corresponding to maximum 


* Ordinate is (L/D) times horizontal specific force in Fig. 9.9 for 
glide vehicle. 


mu 


stagnation point temperature.” 


Comparison of the approximate solution derived herein to numerical 


(15) 


solutions computed by Chapman is presented graphically in Fig. 9.8. 


The derivation carried out in Derivation Summary (9-2) and summarized 


below is based on the same approximations made by Chapman in deriving 


(31) 


his Z_ solution and by Eggers, Allen, and Neice for their "Balliegie 


I 
Trajectory". The solutions of this thesis, however, extend the work of 
the cited papers to include guidance quantities expressed with either 


density ratio, range, velocity, or time as the independent variable. 





e 
(1) Density ratio 0 as independent variable: (C= aa ~e 


(a) Altitude: h(o) = -(1/k)1no (9-54) 


(b) Velocity: ” vig) = ea nae (9-55) 


Veelk sin 1%, 
(c) Range: 
Xy(O) = Xy. + cot %, In eee ~ Xs a 
| (9-56) 
(2) Altitude as independent variable: 
(a) Density ratio: o(h) = ae (9-57) 


* Maximum stagnation point temperature occurs when v = 0.846 v, 
(See Derivation Summary 9-2). 


** For convenience in writing this and all subsequent equations, the 
following quantity is defined: 


y -C(SE)*(m)o ° 


2M 
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(b) Velocity: (9-58) 
Uc, (e~kh_g-khj ) 





(c) Range: 


X,(h) = Kus + cot TC in ee) = Xs + (h-hj) cot tes 





(9-59) 
(3) Velocity as independent variable: 
(a) Density ratio: 
k sm © 
O (ar)= — ia = ap ce (9-60) 
UCpy Ww} 
(bo) Altitude: 
{ k sin Tb VU te 
--1y, (SSP), Zee (9-61) 
A) ln ( hy 
(c) Range: 
{ k sim %, ww 
=e ln = TO} 
Xy lard =Xy.t cot % In ale 
(1+h,) 


Kul) ky; cot Yo | f In (sate ln yt Si) + hi| 


(4) Dimensionless time Tas independent variable: 








Note: Solution is valid if |UCp « 26 
k sin 7, 
(a) Density ratio: 
—UCDS6 








Bei A(c k sin Yr) = In|e +). eee nen! 


n=1 


odLe, 


(9-62) 


(9-63) 


(9-64) 


(5) 


(ob) Altitude: 





UlCp -khi = “UCR WM (enkh_ enki) 
—(T-T; ) kar; sin Mp eK singh © =kthi-h)+ ) | 2 + (9-65) 
Kk sm % 


nen | 
N= { ’ 





(c) Velocity: 














“UCD G; ¥ a) ae 
Be eee. k cin % =), | Sees » ae (9-66). 
(T-T,) ka; sinh e In UC 6; In ll | = 
(d) Range: 
Cae) 
, in 4 -67) 
~(T-Tkny sine P= te 
nee UC hor yk tan t (Xs Xq.) 
—k (XyrXy, ) tan i, + WS (canis) a 
= nn | 
Range Xy as independent variable: 
(a) Density ratio: 
~k(Xx, = X. 
o (Xy) = 6 se he Lament (9-68) 
(6) Altitude: 
; f 
n(X,) = (12 +h, en ~ Ana) tan? a — 
. 


(c) Velocity: 


up fg ehlinthpe Ker am 





ar (X,)2 ai e ns (9-70) 
UC -k [tan % (Xy-Xy. + hi 
D b\ AN SN l 
S k 2 | )-5;| (9~70a) 
ar \ Kn) ar, e 
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As a typical example of a ballistic trajectory solution, consider a 


vehicle entering the karth's atmosphere with: 


hs = POO =i Les 
ana 
Ay, = 0 


The solution for total range from Equation (9-59) is: 
oe = 760 statute miles 


From Fig. 9.6, it may be seen that the vehicle would experience about 
18 g's maximum specific force at v = 0.6. The velocity spectrum and 
time of flight depend on the drag and density characteristics of the 


vehicle. 


9.6.2 Solution of The Glide Trajectory 


The glide trajectory is characterized geometrically by a monotoni- 
cally increasing negative flight path angle. The flight path angle 
initially starts at zero and remains small throughout most of the trajec- 
tory. The flight path angle increases to appreciable magnitudes in the 
terminal phases of the trajectory after velocity has decayed to subsonic 
values. The gliding trajectory occurs when the flight path angle is 
shallow enough and velocity is small enough that vertical accelerations 
and the vertical component of drag force are negligible in comparison 
to other terms in the dynamical equations of motion. In equilibrium 
gliding flight, there is a balance of lift, gravity, and centrifugal 


forces. 


oan 


The limitations of the approximate solution presented here when 
compared to more exact computer solutions of the complete system of 
equations are summarized below: 

(1) The approximate solution derived here is accurate for 

[x L/D = 30 for dimensionless velocities greater than 0.2. 
For the planet Earth, this corresponds to (L/D) equal to or 
greater than 1.0 and velocities greater than 5200 ft/sec. 
(Mach > 4.5). The assumption that lift and drag coefficients 
are independent of Mach number is generally accurate over this 
velocity spectrum only”. 

(2) The approximate solution is accurate over the velocity spectrum 


vy cee bet on Jk L/D = 15. This corresponds to (L/D) Oe. 


Earth — 


(3) The approximate solution is accurate for v)> 0.7 for 


j/k L/D = 7.5. This corresponds to (L/D) garth = 0.25. 


Comparison of the approximate solution derived herein to numerical 
solutions computed by Chapman is presented graphically in Fig. 9.9. 
The derivation carried out in Derivation Summary 9.3 and summarized 


below is based on approximations originally made by Sanger ‘>*) . These 


Cl. 


approximations were also used by Eggers, Allen, and Neice 


5), 


discussing the "glide trajectory", and by Chapman in deriving his 


* The worth of generalized numerical solutions of entry trajectories 
for velocities less than Mach 5.0 may be questioned since Mach 
Number variations of lift and drag coefficients become significant 
at these lower velocities. The variation of ©, and C, with Mach 
Number for Mach less than 5.0 is difficult to ~ generalize because of 
the strong dependence on design characteristics of the vehicle. 

Most authors assume lift and drag coefficients are independent of 
Mach Number effects throughout the entire velocity spectrum. 
Exceptions, of course, are studies of a particular vehicle config- 
uration using Cr and Ch data from experimental tests. 
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ir solution. The derivation summarized here, however, extends the 
solutions of the above authors to expressing guidance quantities in 
terms of independent variables such as dimensionless time, density 
Ratio, etc. 
(1) Density ratio J as independent variable: 
(a) Altitude: 
Equation (9-54) 
(6) Velocity: 


"§ = Grove, Bayt (9-71) 


(c) Flight path angle: 


nae = -- (1+UC) (L/D ) 


(9-72) 
~ k(L/D) 


(d) Range: 
aale)- ic 6 ; [tUCy EO 
ny bo =Ay tT op nein TFUCp & a | (9-73) 


(2) Altitude h as independent variable: 


(a) Density ratio: 
Equation (9-57) 


(bo) Velocity: 


al: 
“khy 3 (9-74) 


vgh) = Gl +UC = e 


(c) Flight Path Angle: 


L -kh 
ie ee a: + UCB ewer) 


I 
KD 


(9-75) 


aa 


(3) 


(4) 


(d) Range: 


web 


Ky(h)=Xi+ . [thm -: ln oa 





[E 
I+ UCp 5 € 


Flight Path Angle as independent variable: 


(a) Density ratio: 





kv { 
OFC UCD & 
(bo) Altitudes 
1 Fu S42) 
h(D=— In ow 


(c) Velocity: 


(2 _ 
avg 1 )= kk c(E Eytan | 


(d) Range: 


tp 


} 
+ —-=— 
In Ks tan v 


(1-179) 


1k 
2 D 


Velocity (vg) as independent variable: 


(a) Density ratio: 


(bo) Altitudes 


h (nrg) =e in| Gece vet, | 
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(9-76) 


(9-77) 


(9-78) 


(9-79) 


(9-80) 


(9-81) 


(9-82) 





(c) Flight Path Angle: 


as 2 ie 
Xn) = Tan — (9-83) 


(d) Range: 


ae) | 
2 nae 
AN (mp)= Ky. +5 p lice . (9-84) 





(5) Dimensionless time as independent variable: T = 0 
i 





(a) Density ratio: oT \ 2 
7 Ata) — (1-ag,) e 
ae (I+ ay, ) + (1~ 0g, e Ub 
io : Ge (9-85) 
ae (Ita )— (-ag: )e ® 
a 
(tag. ) + A-ag; de 1D 
(bo) Altitude: 
n(v) = - dan [o(7)) (9-86) 
(c) Velocity: 
Ps 
U0.) % 
(I+ vg; ) 
A) (9-87) 
(-1¢)) : =a 
(d) Flight Path Angle: 
2T v2 
-2 | Uta.) +C-ag,)e 
Van = =. | — rr (9-88) 


L 2T 
D (+ ng)- U-ag, Je L/D 


jar 


(e) Range: 


x| sls 


(1+ Up, )-(1-ng, Je * “| 
(: ({ + arg. \t (arg, ) c We 


wo Ky t 5 = In (0g) (9-89) 


(6) Range as independent variable: 


(c ueDensisy ratio: 
(1-2) HH? 


olXy)= SEE Lar aE 
UC» D A ai vg) e 1D (Xy =n) 


(9-90) 


(b) Altitude: 


h(%y) = - Ean [2 cK) (9-91) 


(c) Velocity: 


5 (Xn-Xn)) * 
Xn 7Xn. 
nig (Ky)= |I-(-ay*)e (9-92) 
(d) Flight path angle: 
2 
bon Y= — 
sat K [1 C-age ) eH ANN? (9-93) 


9.6.3 Solution of The Skipping Trajectory 


The glide solution discussed in section 9.6.2 is valid only when a 
lifting vehicle enters the atmosphere under very special initial 
conditions, viz. zero initial flight path angle at sub-orbital velocities. 


The glide solution does not apply to trajectories with non-Zero initial 


S22 





flight path angles because a finite initial vertical velocity induces 
a trajectory in which the vertical acceleration term in the equations of 
motion is not small compared to lift force. 

A skipping trajectory is the consequence of a lifting vehicle 
penetrating the atmosphere at high velocity with a finite initial flight 
path angle. The magnitude of the skip in terms of the spectrum of 
altitude, velocity, and specific force spanned during the skip is 
primarily a function of the initial conditions. A number of skips in 
the Gas-Dynamic Phase may be anticipated; the number of skips and the 
severity of the skips depend strongly on the initial conditions and the 
lift-drag characteristics of the vehicle. 

The derivation of an approximate analytical solution of the Skip 
Trajectory is carried out in Derivation Summary (9.4). This derivation 
is based on the following assumptions: 

(1) The difference between the components of gravitation and 
centrifugal force in the lift direction is negligible when 
compared to lift force. 

(2) Coreolis force is negligible when compared to lift force. 

(3) The difference between the components of gravitation and 
centrifugal forces in the drag direction is negligible when 


compared to drag force. 


2 
(4) sing = ;cos7% 2 1 - & 


The approximate solution of this thesis is applicable only for 
flight in the Gas-Dynamic Phase and may be applied only to one skip with 
a given set of initial conditions. Subsequent skips may be analyzed by 


considering each skip individually with a new set of initial conditions. 
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Fig. 9.7 shows the regions in which the approximate skip solution is 
applicable, The approximate solution is accurate for lift-drag ratios 


and initial flight path angles that satisfy the following: 
L 
2K | 15 | 77 1.0 : 


A comparison of the accuracy of the approximate solution of this 
thesis to a numerical solution for the first skip of a vehicle entering 
the Earth's atmosphere is presented graphically in Fig. 9.10. Initial 


conditions for this trajectory were: 


LD B40 
Planet : Earth 
{5 = 2° 
Vgi = 0 
Note: 2k|72[= 62 


The fact that analytical approximate solutions to the skipping 
trajectory must be performed in a series of individual piece-wise 
continuous ‘solutions With proper matching of end conditions for each 
skip may appear to relegate such a solution to the "academically 
interesting but impractical" category. This is not entirely true. In 
general, a skipping trajectory for the manned entry vehicle is 
undesirable because of a number of aerodynamic and thermodynamic reasons.* 


* One of the disadvantages of the skip vehicle is the relatively high 
lateral loads that the vehicle would be required to withstand during 
a skip. These loads, coupled with high thermal stresses due to high 
convective heating rates, would require the vehicle structure to be 
stronger and heavier than that of a comparable glide vehicle. 


Eggers (21) concludes that skip vehicles with L/D = 4.0 and 6.0 
cannot radiate heat at rates comparable to the maximum convective 
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There may be guidance advantages, however, to inducing a partial skip 
early in the Gas-Dynamic Phase in order to reduce a large flight path 
angle”* to zero, or near zero, then "turn off" lift to maintain either a 
constant altitude or slowly descending flight path in order that an 
appreciable percentage of the total energy be transferred at high 
altitudes where heating and deceleration loads are less severe. The 
skipping solution derived in Derivation Summary (9.4) and summarized 
below may be used to analyze this partial skip with a reasonably high 
degree of accuracy. 
(1) Density ratio 6 as independent variable 
(a) Altitude: . 
Equation (9-54) 
(b) Velocity 


\ 


Khe eee, | 2 
&-( ” kK (%) (s -6)) 7) 
aw (o)=nre 





heating rate. The skip vehicle with lift-drag ratios in the neighbor- 
hood of 2.0 absorbs less heat than vehicles developing higher lift- 
drag ratios; the former vehicle, however, still absorbs more heat 
than a glide vehicle or a comparable high-pressure-drag ballistic 
vehicle, 


** In order to control accurately the geographic position at which the 
vehicle enters the Intermediate Phase, it is desired that the 
Keplerian transfer ellipse intersect the atmosphere at as large an 
angle as possible. This may be visualized on the physical grounds 
that if the flight path is nearly horizontal, the slightest error 
in altitude or a perturbation in the density of the upper atmosphere 
may cause a large position error at the penetration point. 

It was shown in Chapter 8 that the effective altitude of the 
sensible atmosphere decreases for steep flight paths; therefore, a 
bonus advantage of entering with a large flight path angle is that 
the exponential atmospheric approximation is a more accurate 
representation of actual conditions. 
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(c) Flight Path Angle: 





Y 
21 aa 2 
ro) #0 = +A | (9-95) 


(d) Range: 


F-O)(%+0) 
Xv (o=Xy, t Nahe 


KO] | (740)(% -@) (9-96) 


where 7 (co) is given by equation (9-95) and the following 
quantity is defined: 


le _ 
DG uae 
Q i “4 aeons (9-97) 


(2) Altitude h_as independent variable 


(a) Density ratio: 


Equation (9-57) 
(b) Velocity: 


Gl Hy _ 22d (kh kh | 
arlh) =a i Ce aot | ii 


(c) Flight Path Angle: 


\ 
ee 
eve) * 


rth) fr2te™e*%) : 


(9-99) 


(d) Range: 


Same as equation (9-96) with equation (9-99) ae 
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(3) 


(4) 


and equation (9-97) for @ 


Flight Path Angle as independent variable 
(a) Density ratio: 


k 
o(1)=o6,+ SUC e (%-4°) 
D 


(b) Altitude: 


n(Y) =-<£1n [5 (¥)| 


5 


(c) Velocity: 
ha 
viv) = vse L/D 


(d) Range: 


Same as equation (9-96). 


Velocity v_as independent variable 


(a) Density ratio: 


, Pee) ee: 
6(w)=6+ oe RK 5 nm) 


(b) Altitude: 


- 1 
hGae= - = iia [o (v)| 
(ce) Filaigat Path Angle: 
_, An Ww 
F(w)= 7-5 Ina 
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(9-100) 


(9-101) 


(G=NG2) 


(9-103) 


(9-104) 


(9-105) 


9.6.4 Specific Forces Acting on the Mntry Vehicle: 


The specific force in Earth g's measured by accelerometers carried 

by the entry vehicle is: 
4 
f, = ~(ado ve c >o[2 + (L/D) si (9-116) 
S(m)E 

Specific force may be written as a function of any one of the indepen- 
dent variables given in equations (9-54) through (9-115) by substituting 
the appropriate expressions for v and 6 as functions of the single 
independent variable selected. As an example, the specific force may be 


written as a function of velocity by substituting 6 (v) as follows: 


(1) Ballistic Trajectory (L = 0; Y= constant = >) : 





sin % 
fin)e= = ite Une coat Te a eeu "a 


Com)e 


(9-117) 


(2) Glide Trajectory: 
rl 


_ Fm)o (1-v6*) 7, eo 
tvge = soe ED + wD) (9-118) 





(3) Skip Trajectory: 
Gumn)o Ly2]3 ky ar Ly ow 
£ (ar), = ee it { H{2)*|*)6, UCptS INS (2q-= nZ) (9-119) 
E 

Since 6 43 <<6 during most of the ballistic trajectory, Equation 
(9-117) shows that the total specific force is essentially independent 
of vehicle drag characteristics, vehicle mass to area ratio, and 
planetary surface satellite dynamic pressure. It is a function primarily 
of the exponential decay characteristics of the planetary atmosphere, 


velocity, and flight path angle. 
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It is interesting to note that the total specific force for the 
glide trajectory, on the other hand, is independent of the decay 
characteristics of the planetary atmosphere. Specific force decreases 
as (L/D) is increased for the glide vehicle. 

The approximation 6; <<6 cannot generally be made in the case 
of the skipping vehicle since altitude may not change by orders of 
magnitude during any one skip. At the bottom of the skip, however, 
this approximation is reasonably accurate; hence, the specific force 
experienced at the bottom of the skip is independent of Cp» vehicle 
frontal loading, and surface atmospheric density. It is a strong 
function of the planetary ec teGeneric decay characteristics, velocity, 


and flight path angle at the beginning of the skip. 


9.6.5 Maximum Specific Force 
It may be seen from equation (9-116) that the specific force is a 


maximum when: 


a(v*o ) _ 
See = 0 (9-120) 


Therefore, specific force is a maximum when: 


Z 
vi= a sin Y (9-121) 


Using equation (3) of Derivation Summary 9.2, it is seen that 


maximum specific force of the ballistic vehicle occurs when: 


v = 0.607 vs (9-122) 


The magnitude of the maximum specific force for the 760 mile ballistic 
vehicle discussed in section 9.6.1 is readily determined by substituting 
equation (9-122) into equation (9-117) (assuming v; = 1.0). A maximum 


pee 


force of 17.3 g's at v = 0.607 was computed in this manner; this result 
compares well with a value of 18.3 g's at v = 0.59 determined from 
numerical solution of the non-linear equations of motion. 

Solving equation (9-121) for the glide vehicle results in maximum 
specific force at v = 0. This result could have been predicted by 
noting on Fig. 9.9 that the approximate solution of this thesis has its 
maximum ordinate value when the abscissa equals zero. Numerical 


calculations, on the other hand, show that: 


(1) For Jk L/D 


ape Vg at fai 0.23 (Corresponds to 
(L/D) po eth = 2-9) 


(2) For [x L/D 15, Vg i — = 0.29 (Corresponds to 


(L/D) porth = 0-5) 


= 0.35 (Corresponds to 
vw 
(L/D) garth = 0.25) 


(3) For Jk L/D = 7.5, vg at fag 


x 


12 


0.45 (Corresponds to 
= 0.2) 


(4) For {k L/D =—om0, Vg at A 
(L/D) arth 
In general, the velocity at maximum specific force for glide vehicles 


in excess of (L/D) porth = 1.0 occurs at velocities less than 0.2. 


It was shown in Derivation Summary (9.4) that flight path angle 
is a convenient parameter for determining the point of maximum specific 


force for the skip vehicle, 


mk |y_fpafel)*(y? a 2 Re) 2 
aie A A+@2)*(x, a 6 Fe a ian (9-123) 


For small initial flight path angles and large lift-drag ratios, the 


. © fe od ° e 
maximum specific force occurs when (C= QO. The velocity at maximum 
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specific force for the skip vehicle is: 


Be % \2 
(nr) = © {3% 5 (75) I (9-124) 
max 


9.6.6 Stagnation Point Temperature 


The stagnation point temperature was written in equation (7-21) 


as Pol leus: 


eS a il 
8 


T, = 1,392 X 10*(HF) - (VF) oom) °Rankine (9-125) 


where (HF 59 is a "heating function ratio" of the atmosphere of planet 0 
with respect to that of the Earth,’ and (VF) is a vehicle function 
depending on the emissivity of the skin structure and on the radius of 
curvature in the vicinity of the stagnation point. 

Stagnation point temperature may be written as a function of any 
one of the independent variables given in equations (9-54) through 
(9-115) by substituting the appropriate expressions for v and 9 as 
functions of the single independent variable selected. As an example, 
the stagnation point temperature may be written as a function of velocity 


by substituting 6 (v) as follows: 


(1) Ballistic Trajectory (L = 0; Y = constant = Y,) : 


7 diners; w[kinthre a \S , 
Ty (ar)= 1.592 X10 (CHF) © (VF) eres In R (9-126) 
(2) Glide Trajectory: 
\ arnt ll—-n; 28 
Te (7)=1.392X10° (HF) 6 (VE) “WG, )* -R (9-127) 


Be 


(3) Skip Trajectory: 


®@ 


| 1 a ls 4 
T, (wr) = 1.392 x10* (HF) (VE) * ar ai Buc, eh 5 In uh ) 


oR (9-128) 
The stagnation point temperature for all three trajectory profiles 
is,aetunction of Cr and/or Cp » vehicle mass to area characteristics, 
and surface satellite dynamic pressure. It is interesting to note, 
however, that the stagnation point temperature of the glide vehicle is 
independent of the exponential decay characteristics of the planetary 


atmosphere. 


9.6.7 Maximum Stagnation Point Temperature 


From equation (9-125), stagnation point temperature is a maximum 
when: 


d 
Selva ) = 0 (9-129) 
This corresponds to: 


(me 
iy ey etnd (9-130) 


For the ballistic vehicle, stagnation point temperature is a maximum 
when 

y =o .046 V5 (S213) 
This result may be used in equation (9-126) to give the maximum tempera- 


ture that may be expected during the ballistic trajectory. For example: 


(1) a 0 
Ch = 1.0 
Y= -6° 


oa 





(VF) = 1.0 
M/S =.1,0 
Planet Earth 


Te for this situation is 4850° Rankine. 
max 


(2) If all conditions are identical as in (1) above 


except Ch 


= 10, then Ts wax = 3650° Rankine. 

By differentiating equation (9-127) with respect to T and setting 
the result equal to zero, the velocity at maximum stagnation point 
temperature of the glide vehicle is determined. This occurs when: 

vg = 0-815 7 (9-132) 

The maximum temperature for the skip vehicle is incurred at the 

following flight path angle: 


(1). - 45 1- A +(SP)" (G+ 2 Yea) (9-133) 


Max 


(Y) ab (een) | 


eran > (9-133a) 
This corresponds to a velocity of: 
vr, =. ote eM } 
geen t (9-134) 


Comparing equations (9-133) and (9-123) it is seen that for the skip 


véhicle: 


or, (9-135) 


max max 


Equation (9-135) shows that in a typical skip trajectory, the maximum 


temperature level is incurred before the flight path angle reduces to the 


boo 


value corresponding to maximum specific force level. Comparison of 


equations (9-124) and (9-134) shows that: 


V, - e *Céonan| -%s Uc was) 
ae 


Equation (9-136) demonstrates that the velocity at maximum temperature 
is greater than the velocity at maximum specific force level during the 


skip. 


9.7 Summary 

Atmospheric entry trajectories were examined in Chapter 8 from a 
unique and instructive vantage point, namely by studying the behavior 
of the Conservation Parameter 73 | . Examination of the conduct of 
this parameter with respect to altitude, specific force level, energy, 
and angular momentum led to the conclusion that the entry trajectory 
could profitably be analyzed in three distinct operational regimes; 
these regimes were named the Keplerian, Intermediate, and Gas-Dynamic 
Phases. 

The investigation discussed in Chapter 9 was oriented toward 
utilizing the Conservation Parameter for determination of guidance 
quantities. Two dissimilar entry profiles were considered: 

(1) The direct entry profile results from perturbation of a stable” 

satellite orbit by means of retro-rocket thrust. In this 
profile, the vehicle generally traverses in sequence the 


Keplerian, Intermediate, and Gas-Dynamic Phases to a landing. 


* Or relatively stable. 
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(2) The degenerate orbital profile consists of a series of braking 
passes through the outer reaches of the atmosphere in order to 
reduce the energy level preparatory to final direct entry. In 
this profile, the vehicle enters the Intermediate Phase in a 
series of near impulses until the energy level is reduced 
sufficiently for final transition to the Gas-Dynamic Phase. 

The feasibility of using the Conservation Parameter as a switching 
function for profile (1) and as a prediction function for profile (2) 
was discussed. 

Range capabilities and range accuracy in the Keplerian Phase were 
Beisciosed in terms of geometric quantities at the trajectory modification 
point for entry from elliptical satellite orbits, and in terms of 
velocity impulse and engine gimbal angle for entry from circular 
satellite orbits. Examples of range sensitivity to errors in operation 
of the retro-rocket system were computed and presented in both 
graphical and analytical form. 

Approximate solutions of various guidance quantities were derived 
for ballistic, glide, and skip profiles in the Gas-Dynamic Phase. 
Limitations on the approximate solutions were established after comparing 
them to numerical solutions obtained from the non-linear equations of 
motion. Analytical determination of instantaneous and maximum decel- 


eration loads and temperature levels was discussed. 


Bo / 


ps 


Derivation Summaries 


AL1 Derivation Summaries referred to in the previous 


of this chapter are appended. 


Derivation 
Summary No. Title 
Pie Solution of Keplerian Phase of Entry Trajectory 


Oe 


9.3 


934 


in Terms of Geometric Quantities at the 
Trajectory Modification Point. 


Solution of the Ballistic Trajectory 


Solution of the Glide Trajectory 


Solution of the Skip Trajectory 
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Derivation Summary 9.1 


Solution of Keplerian Phase of Entry Trajectory in Terms of 
Geometric Quantities at the Trajectory Modification Point. 


Angular momentum is conserved in a pure Keplerian trajectory; 
p is a convenient constant of integration to use in solving the 


equations of motion for the Keplerian Phase of the trajectory. 


oe *YI1¢ (1) 


Assuming planar motion” , the horizontal component of velocity with 


4 


respect to inertial coordinates is: 


YI¢ 5 rg! (2) 


Therefore, angular momentum is written: 


2 
p=r @' (3) 
Using equations (3) and (9-4), equation (9-5) is written: 
‘ aL 
roe +26 o (4) 
Yr p a 


This equation may be written in a form suitable for integrating by 
making the following transformations 

We i (5) 
Therefore: du = -dr/r* (6) 


edu 


r' = -r“u! = "a Gt = PS (7) 


C-- —  ~ —  —  —-  ~ - ~~ - a - ~~ — | 


* See Section 6.1 for a quantitative discussion of the effects of 
non-spherical gravitational components in limiting the assumption 
of planar motion. 


oo 





Derivation Summary 9.1 (cont.) 


Equation (4) is written with the aid of equation (8) as: 


Be L 

OUs,y = + = constant (9) 
2 2 

dg p 


The solution of equation (9) is as follows: 
= a 
u(¢) = Cc) cos(@ + Co) i (10) 
It is convenient to express u in terms of the true arlomaly: 


1) =n +0 (11) 
Y ds the angle measured from the ascending line of nodes to the 
line of apsides; 8 is true angmaly. If the entry transfer ellipse 
is fixed in inertial space, then 2” is constant. Therefore: 
= ae 


p 
Constant C2 is determined by noting that the time rate of change of 


radius (altitude) is zero at perigee and apogee. From equation (7): 


r' = at = 1 C, sin(@ + C3 )0! (13) 


Therefore: 


@+C3=0+n7 (n= 0,1,2...) (14) 


From Fig. 4.1, it may be seen that 6 =7[ at apogee and 6 = 0 
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Derivation Summary 9.1 (cont.) 


at perigee. The equation of the ellipse is given in Fig. 4.1 in 


terms of R, a, 8, and €. Writing this in dimensionless form: 


88 (15) 


In equation (15), ay is the semi-major axis non-dimensionalized with 


respect to Remo Equation (15) may be written: 


a 
= = —s (1 + € cas@) (16) 


Bn 


where by is dimensionless semi-minor axis of the ellipse. With 


Co = 0, equation (12) and equation (16) are used to eliminate r: 


pe = BN (17) 
N 
An 


The constant C, can be determined by differentiating CZK 
i v 
a a (19) 
Using equations (9.4=-2) and (9.4—3) given in Fig. 9.4, this may 


be written: 





c, sine = tan Fz (20) 
Tt 
From equation (12): 
a i 
C, cos@ == - —————_-— (21) 
1 if révzScos* 7 I 
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Derivation Summary 9.1 (cont.) 


Squaring and adding equations (20) and (21) makes it possible to 


eliminate 8. Using the result in equations (17) and (18) gives: 
ae v2 peos* I y(2 = vrr ) (22) 


Equation (22) expresses the eccentricity of the entry transfer 
ellipse in terms of velocity, radius, and flight path angle. The 
magnitude of these quantities at the trajectory modification point 
is denoted by the subscript m. 

The semi-major axis is determined from equations (17) and (22) 


and the following relation: 


enn) ae (23) 


The resulting equation in terms of quantities at the trajectory 
modification point is given as equation (9-8). Equations (9-9) 
through (9-18) are readily determined from the foregoing results 
and from the equations listed jm Fig. 4,1. 

Time of flight may be derived in terms of © by noting from 


equation (3) and (11), with n = constant: 
@' = p/ré (24) 
Using equations (16) and (17), this is written: 


2 
Ai (1 + Ecos8) (25) 
p 
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Derivation Summary 9.1 (cont.) 


Integrating equation (25) gives the time of flight listed as 
equation (9-19). 

The velocity at any instant in terms of quantities at the 
trajectory modification point is written from equation (9-8) with 
ay = constant. Flight path angle may be written in terms of 
quantities at the trajectory modification point from equation 
(9.4-3) with angular momentum equal to a constant. 

Range is defined by equation (9-23). With the aid of equations 
(17) and (18), equations (9-24) and (9-25) are written. The true 
anomaly may be eliminated from the range expression by taking the 
cosine of equation (9-23) and substituting equations (9-24) and 
(9-25) into the result. The range expression so derived is given 
by equations (9-26) and (9-27). 

With equations (9.3-1) through (9.3-5) given in Fig. 9.3, the 
foregoing results are written in terms of velocity impulse $ar, 
and engine gimbal angle A, for the special case of a circular 
reconnaissance orbit. The results are summarized as equations 
(9-29) through (9-35). 

Equations (9-36) through (9-42) are derived by writing for 
equation (9-3): 


cos 8 = cos Xy cos @, - sin Xy sin @, (26) 


© and @, are eliminated from equation (26) by using equations 


(9-24) and (9-25): 
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Derivation Summary 9.1 (cont. ) 


ane | Za 
= [p -r| = cos Xy (oe Y ,) 7Sinky [EP rae (p - ry) 7) 
Substituting for eccentricity from equation (9-7) and for angular 


momentum from equations (9.4-3) gives equation (9-36). This may be 
written in the following alternate forms: 


2 _ (1 - cos Xy) 


Vv (28) 


cos Ty [x €08 Vom -TY ,cos (Xy+ Vr) 


xX 
Peos(2 Vom) - cos(2 ~ + Xy) = 20 -og fs - cos ry (29) 


in Ym 
Substituting equations (9-29), (9-31) and (9.3-5) into equation (26) 


gives equation (9-39) or its alternate form (9-40). Equations (9-41) 


and (9-42) were derived by implicit differentiation of equation (9-39). 
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Derivation Summary 9.<¢ 


Solution of the Ballistic Trajectory 


The most convenient set of planar equations for solving the 
ballistic trajectory are equations (6-40) through (6-44); these were 
written in terms of components in the ly and ly directions. 
Assuming powerless flight (Ty = 0), it may be seen from equation 
(6-41) that if the centrifugal force and Coreolis force terms in’ the 
ly direction are balanced by gravitational force, then: 


3 


y' = 0 (1) 
Therefore: 
Y = constant = ies ("ballistic" flight path angle) 
(2) 
If the gravity force in the 1. direction is small in comparison 


X 


to drag force, equation (6-40) gives: 


V'Gs.- Np = = UC.v“o (3) 
where: 
PS Was 
y= —D _ia)o (4) 


The time rate of change of density ratio is: 
g' =- ke ane (5) 
Therefore, from equations (3) and (5): 


eee we (6) 
dg k sing 


a 


Derivation Summary 9.2 (cont.) 





Integrating equation (6) gives: 


UEC 
Dane — sinnerman ane a leergees: 6.) (7) 
Vv, k sind, i 


The solution for range is readily determined from equations 


(6-43) and (6-44) with Cai: 


KX, = X 


ny ng = cot ty in = = cot Tp (n-h, ) (8) 


Bi 
Atmospheric density ratio may be solved as a function of time by 


using equation (5) and (7): 





a k sin % 

Tero = —kay sin fp e (9) 

ksin % 

Ove 
Integrating equation (9) gives: 
~VE no; 
k sin Vp Z 
—(-T,) kui sn % € =m|Z|— Yoo (¢_¢.) 
: l Oj l 
k sin % 


+(2S \ (o*- 6") | (2 ) (6°-6;") 
k sin % a. k sin % nen | 


(10) 
Equations (7), (8), and (10) are three fundamental equations from 


which all of the equations (9-54) through (9-70) are derived. 
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Derivation Summary 9.3 


Solution of the Glide Trajectory 


If the flight path angle, 7 , at the onset of the Gas-Dynamic 
Phase is zero, then the vertical velocity at this point is zero. 
Basic assumptions required for equilibrium gliding flight to exist 
throughout the trajectory are: 

(a) Vertical accelerations are small, i.e., v,'—» 0. 

(b) Vertical component of drag force is small; 

=. UC hv vo 0 


With the above assumptions and the following approximations: 


|E|< 1.0 (1) 

V 
Gi eID s) ~ 1.0 (2) 
vig = vg (3) 
r= (1 +f) = 1.0 (4) 


equation (9-45) is written as follows: 


ZC 
V 
22 £) -be vou (5) 


Substituting equation (5) into equation (9-46) gives: 


2] 
Les) 
Vg" = “TED (6) 
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Derivation Summary 9.3 (cont. ) 


For L/D equal to a constant, equation (6) may be integrated to give: 


2(T-T;) 
(7) 


=1 , +1 
(Vg ) (Vg ) . 
Solving this explicitly 


Equation (7) holds only when Vy x a 


for Vg gives equation (9-87) 
Flight path angle may be computed by differentiating equation 


(5) with respect to Vg 
) (8) 


L dary | do 
5 UCpew(s. $F O dary 
Using equation (6) above and equation (5) of Derivation Summary 9.2 


(9) 


=a ( + 


Comes K(=) arg tan © 
2 


1 dé _ 


Also, the following may be written 
(10) 


aha lee ve At ) 
g 


The second term in the parenthesis of equation (10) is the product 


small quantities and can be neglected in comparison to l 


of two 
Using equation (9) and (10) in equation (8) gives 


tan EY = REDDY? 


Range is determined as a function of altitude from 


(Ce 


equations (9-43) and (9-44); 


Derivation Summary 9.3 (cont.) 


ax 
N _ i v 1 (ae)) 


— 


dh rtany tany 


Using equation (11) for tan Y and equation (5) for vg6h) gives: 


oy -k L/D 1 es 
an . (1+vC, (L/D) ao) 


Integrating equation (13) gives equation (9-76). 
Range as a function of velocity is determined from equation 


(9-43) and (6): 


ax Vv Vv 
N L/D... (Lip) 
dv, = 2 rot) - ee a) 


Integrating equation (14) gives equation (9-84). Velocity as a 
function of range may be written by solving equation (9-84) 
explicitly for Vg ° 

Equations (5), (7), (11), (9-76), and (9-84) are fundamental 
equations from which equations (9-71) through (9-93) were derived by 


algebraic operations. 
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Derivation Summary 9.4 


Solution of the Skip Trajectory 


i. 


The glide solution discussed in Derivation Summary 9.3 is valid 
only when a lifting vehicle enters the planetary atmosphere under 
very special initial conditions, viz. zero initial flight path angle 
at velocities equal to or less than orbital velocity. The glide 
solution does not apply to trajectories with non-zero flight path 
angles because a finite initial vertical velocity induces a trajec- 
tory in which the vertical acceleration term in the equation of motion 
is not small compared to lift force. 

The following basic assumptions are made in @riving the 
solutions for the skipping trajectory; i.e., the trajectory of a 
lifting vehicle entering the planetary atmosphere with a finite 
initial flight path angle: 

(1) The difference between the components of gravity and 

centrifugal force in the lift direction is negligible 
when compared to lift force. 

(2) Coreolis force is negligible when compared to lift force. 

(3) The difference between the components of gravity and 

centrifugal forces in the drag direction is negligible 
when compared to drag force. 

Using assumptions (1) and (2) in equation (6-41), and assuming 
powerless flight ({y = 0), gives: 


ve" = ny (1) 
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Derivation Summary 9.4 (cont.) 
Using assumption (3) in equation (6-40) gives: 
v' = -np (2) 


Dividing equation (2) by equation (1) gives: 


Di (3) 


1b 
aia 


Assuming (L/D) = constant and integrating gives: 


“(7 te) 
m= 2? L/D (4) 


Equation (4) gives velocity as a function of flight path angle. 
Dividing equation (1) by o' given by equation (5) of 


Derivation Summary 9.2 gives: 


av --«;UGp (L/D) 


6.6 KK sin (5) 


Integrating equation (5) for a constant lift-drag ratio gives: 


UC 
cos J = cos ie = = (igi) (Oo - 0) (6) 


For small flight path angles: 
Z 
cost = lee Le (7) 


With equation (7), equation (6) is written: 


6=S8, tage (Ty - 7") (8) 


Jol 


Derivation Summary 9.4 (cont.) 


Equation (8) relates density ratio (altitude) and flight path 
angle. The rate of change of range with respect to density ratio 
is determined by dividing equation (6-44) by density ratio given by 


equation (5) of Derivation Summary 9.2: 


), 





ON cot 1 
ao 


At altitudes where skipping may occur, the distance from the center 


of the planet to the vehicle differs little from the mean planetary 


as 


radius, hence r=1.0. With this assumption and the small angle 
assumption on cot ¥ , equation (9) is written: 
ONy lL __ 
do Oke 


tie 


(10) 


Equation (8) is used to eliminate % from equation (10). Integra- 
ting the resulting equation gives range as a function of density 
ratio. This integrated expression is listed as equations (9-96) 
and (9-97). 

Flight path angle as a function of time can be @étermined by 


writing equation (1) as: 
ao Uve, o (11) 


Assuming that the velocity change during a single skip is small and 
substituting ao from equation (8) into equation (11) gives an 


expression for Y' which may be integrated. The resulting solution 
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Derivation Summary 9.4 (cont.) 


for Y(t) is listed as equation (9-107). 

Equations (4), (8), (9-96), and (9-107) are fundamental 
equations from which equations (9-94) through (9-115) were 
determined. 

The flight path angle at maximum specific force level is 


determined by substituting equation (2) into equation (9-121): 


20UC)h = - k sin 7 at fy (12) 


Thus, at maximum specific force level: 


2 ae ie (13) 





Substituting equation (8) for 6(% ) into equation (13) and 
solving explicitly for v gives equation (9-123). The velocity 
at maximum specific force level is determined by substituting 
equation (9-123) into equation (4). 

The flight path angle at maximum stagnation temperature level 
is determined by substituting equation (2) into equation (9-130). 
The velocity at maximum stagnation temperature level is determined 


by substituting the resulting solution for f into equation (4). 
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Chapter 10 


APPROXIMATE ANALYTICAL SOLUTION OF GUIDANCE PARAMETERS AND 
CONSTRAINTS FOR THE DEGENERATE ORBITAL ENTRY PROFILE 


ion. sintroduction 

The investigation described in Chapter 9 was devoted largely to 
derivation of guidance parameters and constraints for the direct entry 
profile. This atmospheric entry trajectory is characterized by 
sequential transition from Keplerian flight through the Intermediate 
Phase into the Gas-Dynamic flight. The Conservation Parameter ral 
was convenient for defining the boundaries existing at transition points 
between the three separate operational regimes encountered during the 
course of entry. 

The degenerate orbital profile is examined in Chapter 10 from an 
entirely different standpoint. This entry profile has a time history 
far removed from that of the direct entry profile. A series of braking 
passes through the outer reaches of the atmosphere are employed in this 
concept of entry for purposes of reducing the energy level of the 
vehicle until final direct entry is assured. The Intermediate Phase is 
entered in a series of near impulses in the vicinity of perigee; during 
each of these pulses, energy is transferred to the planetary atmosphere 
until ultimately the energy level is low enough that final transition 
to flight in the Gas-Dynamic Phase evolves. 
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The characteristics of powerless flight in the outer reaches of the 
atmosphere are examined in Chapter 10 for vehicles initially in cir- 
cular and elliptical reconnaissance orbits. For this analysis, the 
entry mission was assumed to start at the first perigee. Analytical 
techniques are presented for predicting the range and time of flight 
of the craft, the rate of circularization of elliptical orbits, the rate 
of decay of circular orbits, the number of orbits remaining under 
specified initial conditions, perigeal and apogeal decay rates, and the 
effect of the vehicle's aerodynamic and mass characteristics on these 
quantities, . 

It is shown in this chapter that a true circular orbit or a linear 
decaying circular orbit cannot exist, even under idealized conditions 
of injecting a vehicle exactly at circular orbital velocity in the 
vicinity of a spherical planet. Because of the influence of the 
planetary atmosphere on the dynamics of energy transfer, the altitude 
and flight path angle in the high altitude circular orbit oscillate at 
orbital frequency with such small damping as to be essentially undamped. 

The analysis of the circularization phase of elliptical reconnais- 
sance orbits described in this chapter shows that the drag character- 
istics of the vehicle are important in specifying the resulting trajec- 
tory and that the lift characteristics of the vehicle are relatively 
unimportant. The altitude at apogee decays in proportion to the sum of 
two large constants which are almost equal in magnitude* while the 
altitude at perigee decays in proportion to the difference of these two 


* The constants are described by modified Bessel Functions of the first 
kind; i.e., Bessel Functions of the first kind with pure imaginary 
arguments, 
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near-equal constants. To a first order, therefore, perigeal altitude 
remains constant during the circularization process. 

Guidance parameters important in the conceptual phases of entry 
systems are derived in closed form in this chapter. The solutions 
obtained are compared with more exact machine-computed numerical 
solutions. As a result of this comparison, accuracy and limitations 


of the analytical solutions of this thesis are determined. 


10.2 The Altitude Differential Equation 


The planar equations of motion in the 1. and 1 directions were 
written as equations (6-47) and (6-48) in Chapter 6 of this thesis. In 
powerless flight (Ty = 0) and assuming (1 + L/D tan¥ ) = 1.0, these 


equations reduce to: 


Vic 
tt — i 1 ° 
r gel. Zi Ny sintY + ny, cose Ci0=ay) 
Vv 
Vig" =. 12. - Ny cos (10-2) 
Vv 
Xi" = 3 (10-3) 


Equation (10-1) shows that the acceleration in the ig direction is 

equal to the centrifugal specific force minus gravitational acceleration 
plus lift and drag specific force components in this direction. For 
small flight path angles, the drag specific force in the radial direc- 
tion may be neglected and cos ¥ =1.0. Substituting for ny and ny, from 
equation (6-36) and (6-37) and assuming v = vig , these equations are 


written: 


2 2 
rv Xe 4 +(L/o)o & * (10-41) 
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2 
ar (10-5) 


eg tn 
zs 


nNi< 
| 


© 
} 


Xy' = v/r (10-6) 


In equations (10-4) and (10-5), C is a dimensionless constant for a 
particular vehicle at a given initial altitude above a particular 
planet: 


C is the atmospheric density in slugs per cubic foot corresponding 
to the initial flight altitude. . 
It is convenient to write equations (10-4) and (10-5) in terms 
of Or and Av, which are defined as changes in radius and velocity 


from given initial values 5 and ee? 


Ar ee (10-8) 


ee (10-9) 
Since rj; andv, are constants in equations (10-8) and (10-9), the 
rate of change of Ar is equal to the rate of change of r and the 
rate of change of Av is equal to the rate of change of v, i.e.: 

Ar'=r' (10-10) 

Av' =v' (10-11) 


With equations (10-8) through (10-11), equation (10-4) is written: 





arte ME [ieee 5. (ary ] f 


L 


+H & [14 Caw (ary) 


(10-12) 


oot 


In the Keplerian and Intermediate Phases of the trajectory and through 


Z 
much of the Gas-Dynamic Phase of flight, Con and (S2)* and all 
a 


higher order terms can be neglected; therefore equation (10-12) is written: 


whet ote at +28 bene 


Pc 
+ hc at p 
D Zz OC: 





(10-13) 


v in equation (10-13) represents the dimensionless circular orbital 


ci 


velocity at the initial altitude: 


a (10-14) 


ei 2 
+ 


There are three time varying quantities in equation (10-13): Av, Ar, 


and @ . Equation (10-5) may be used to give an expression for Av: 


| 
a: = . e (lowals) 


Integrating equation (10-15) with Ts = 0 gives: 


Z ic 
1 V 
Av=- Arc 4 fo ~ ALL ar (10-16) 
rs ee Ps, 


The first term on the right hand side of equation (10-16) represents 
the change in velocity resulting from kinetic-potential energy trade- 


off; the second term represents velocity loss due to drag. 
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Substituting equation (10-16) into equation (10-13) gives: 


_ Vi “ec (Mey bs | fe ic 
ar Fesh eae Sate + pcre) sees 
ge G 
™ 


ot ) f(t e)f dr 


(10-17) 
Equation (10-17) is the basic equation which must be solved in order to 
determine the behavior of Ar with time. It is non-linear because of the 
dependence of P onr. If the exponential model of the planetary atmos- 


phere is assumed, then: 


C= ek4r (10-18) 
Cy 
Solutions to the non-linear equation (10-17) with the atmospheric 

model (10-18) are difficult to obtain by conventional techniques. It 
is shown subsequently in this chapter that a surprisingly accurate 
analytical solution may be obtained which is valid over reasonably large 
time intervals by a sequence of approximations. This method may be 
described as follows: 

(1) Assume a reasonable approximation for MAr(T) based on the 
initial conditions established for the trajectory. 

(2) With this approximation of Ar(7T) the atmospheric density 
behavior is specified as a function of time for the parti- 
cular atmospheric model chosen, 

(3) Substitute e (7) from step (2) into equation (10-17) and solve 
for Ar (T). This gives a refined approximation for Ar(T) 
which should be more accurate than the assumed solution of 


step (1). 5 


(4) Repeat steps (2) and (3) as necessary. Provided the solution 
converges, the degree of accuracy at any stage may be esti- 
mated by comparing successive approximations. 

This method for obtaining analytical solutions is carried out in 
subsequent sections of this chapter for both circular and elliptical 
reconnaissance orbits. Accuracy of the analytical solutions obtained 
by this procedure is confirmed by comparing them with machine-computed 


numerical solutions for various vehicles and initial conditions. 


10.3 Circular Orbital Entry 

If the vehicle is in a circular reconnaissance orbit at zero time, 
then certain simplifications may be made in equation (10-17). 
Specifically: 


wie dF (10-19) 


Equation (10-17) therefore reduces to: 


Ar ‘= -ar Ct “effi 5Ce Slt + ECE 


coe + $5) ta) 8 dr 


(10-20) 
Since the vehicle is initially in a circular orbit, a reasonable 
approximation to Ar(7) required in the first step of the method of 


solution outlined in section 10.2 is: 


Ar =0 (a10=21,) 


Equation (10-21) is the equation of a circle. If Ar does not change, 
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then regardless of the atmospheric model”: 
Ci ey =i (10-22) 


Substituting equation (10-22) into equation (10-20) gives: 


a's (Ter ben) cote Sb Lee) = pS 


(10223) 
The following constants are defined: 
ei ee 
We = (==) (1 +(2/D) Cry) (10-24) 
i 
be (L/p) Cv (10-25) 
Cc a 
B. = Cv, [ 2/ry +(1/0)$] (10-26) 
Cc = Pe (10-27) 
Cc ry 


The subscript "c" in definitions (10-24) through (10-27) is used for 
purposes of identifying the constants as applying to the circular 
orbital case only. In terms of these constants, equation (10-23) is 


written: 
T 
Ar" +w,“ Ar - C, fazer = A. - BLT (10-28) 
O 


Cc 
It is noted that at the initial point (T = 0): 


* Equation (10-22) assumes a homogeneous atmosphere around the planet; 
i.e., the density at a specific altitude is constant. This, of 
course, neglects diurnal and latitude density variations. 
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LO y= A, 
Ar'(0) = 0 (10-29) 
Lr (O}*="0 . 


Equation (10-28) is solved by means of the Laplace transform: 


i (ar(n)} = RG) =f Ee" ar(r) dT (10-30) 


The transform of the first and second derivatives is given by: 


£ {ar} = 5h(s) — ar(0) 





(10-31) 
i far'cn| = $Rrsy — sar(o) —ar‘o) (10-32) 
Applying the transformation to equation (10-28) gives: 
R(s) = Ae — Fe 
x Paes 16: (10-33) 
Equation (10-33) may be solved by setting: 
so +u,2s - G,= (s+ b,)(s + bg)(s + b;) (10-34) 
where: 
oS Sos (10-35) 
bo a oa - Ly, (10-36) 
bg =X + dy, (10237) 
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The following constants are defined: 


oa. ih etna 3/3 
X= Fe sinh|Z sinh |) We (10-38) 
Yi = uk cosh [sinh 35.) (10-39) 
fe, (10-40) 


Substituting equations (10-34) through (10-40) into equation (10-33), 
expanding by partial fractions, and taking the inverse transform of the 


result gives: 


£4T yr 


_ Ac IX . op 
ar) = Garey eM cour tHE sin yr) 


Be 


ae ee ae 
v 2x, (x2 + y2 UIxt + ys) 9% YE (+42) 


26" [Ae cos T + 22(3x?-y?)si 7] 
1 vf Y, hy Yh PEs 


(10-41) 
Equation (10-41) is the solution of equation (10-17) for circular 
orbital motion. Considerable simplifications may be made in this 
equation by examining the magnitudes of the quantities involved for 
reasonable initial reconnaissance orbital conditions: 


Se: 1.0 (10-42) 
2 to. — a 
ale 


= 


TY. 


363 





x} xz yy (10-444) 
ox. T rw ws 3 BoT 
mw Gi t+ ext Hl re; (10-45) 
iW, 
BY 
eT 2a - 2 C (10-46) 
r.W 2 
Lae 


Substituting approximations (10-42) through (10-46) into equation 


(10-41), and noting that 3x, < w, and 4x, € vw, gives: 


A 
Ar(T) = ae [ (arexyr) Z (1-x,T) cos wr] (10-47) 
C 


Be Saucy. | 
- | 7- (1-41) —— 
Wr, : Wo 


Equation (10-41) and its simplified version, equation (10-47), 
correspond to oscillatory motion of Ar; the frequency of the oscil- 
lations is equal to orbital frequency. These oscillations are damped 
at such a small rate” that, from the engineering standpoint in guidance 
system design, the damping terms can be ignored. Therefore, equation 
(10-47) reduces to: 


A B ; 
eee) = (1 wT) See je 10-48 
r(T) an cos WT) oe (10-48) 


Equation (10-48) is the fundamental solution for the circular orbital 
trajectory. It may be observed that equation (10-48) is obtained 


directly by assuming for the integrand of equation (10-20) that: 


: x7 is generally an extremely small number. 
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(lies <) = 1.0 (10-49) 
i 


Therefore, equation (10-48) is an accurate solution to equation (10-23) 
for changes in altitude which are small in comparison to the initial 
distance of the vehicle from the center of the planet. This does not 
necessarily mean, of course, that it is an accurate solution of the 
true trajectory because of the approximations made in arriving at 
equation (23). Proof of the accuracy and limitations of equation 
(10-48) is carried out in Figs. 10.1 through 10.7. Before discussing 
these results, it is desirable to examine equation (10-48) in greater 


detail. 


It is noted that both A. and B. in equation (10-48) contain lift 
_ ¢ vs? 


Ps 


terms. For the zero-lift vehicle, A, = O and Ip) The alti- 


tude of the zero lift vehicle is characterized by a linear decay on 
which is superimposed a sinusoidal oscillation at orbital frequency. In 
the case of a lifting vehicle, 7 O, the trajectory also has a cosine 
term. Lift, from equation (10-48), has no effect on Ar at the termin- 
ation of each complete orbit (i.e., at W\T= n27i ). It may be seen 
from equation (10-47), however, that if 2x,T becomes significant in 
comparison to 1.0, lift causes succeSsive perigees to be higher than in 
the zero-lift case, 

Equation (10-48) is plotted on Fig. 10.1 for one complete orbit of a 


zero lift vehicle with Cps/M = al ft.“-slug 7. Superimposed on this 


plot 16 Nielsen's**!) 


numerical solution for the same vehicle. It may 
be observed that the solution of this thesis and Nielsen's numerical 


solution are practically identical throughout the orbit. The total 
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altitude loss during the orbit is approximately 760ft. Also shown in 


Fig. 10.1 are mlutions to equation (10-48) for lifting vehicles with 


ft.) 
slug’ * 


L/D = 1.0 and L/D = 2.0 having the same drag parameter = = 1.0 
It is seen that the L/D = 1.0 vehicle climbs about 50 ft. before descen- 
ding while the L/D = 2.0 climbs about 200 ft. prior to descending. In 
any case, Fig. 10.1 and equation (10-48) show that a true circular orbit 
or a linear decaying circular orbit cannot exist because of the dynamics 
of energy transfer even if the problem starts under precise circular 
orbital conditions. 

Fig. 10.2 compares the solution given by equation (10-48) with the 
machine computed numerical solution for the first ten orbits under the 
identical initial conditions used for Fig. 10.1. It should be noted 
that the altitude at the completion of each orbit alone is plotted in 
Fig. 10.2; the sinusoidal oscillations during each orbit have been 
omitted. It may be observed that the analytical solution and the 
numerical solution are remarkably close; at the completion of 10 orbits 
the analytical solution shows an altitude loss of 7,600 ft., while the 
numerical solution shows about 7,800 ft. altitude loss. 

Fig. 10.3 compares the analytical solution of this thesis with 


numerical solutions for one orbit of a zero-lift vehicle with 





Crp mG 
er © 10 ain launched in a circular orbit at 120 miles above the 


Earth. The trajectory has the same general shape as that of Fig. 10.1 
except the total altitude loss is considerably greater here because of 
the ten to one drag increase of the vehicle. The analytic solution 


derived a total altitude loss during the first orbit of approximately 


7,000 ft. while the numerical solution computed the altitude loss to be 
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Numerical Solution for High-Drag Zero-Lift Vehicle 
Initially in Circular Orbit at 120 Miles Altitude 
above the Earth. 
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about 7,750 ft. 

The limitations of the analytical solution in the high drag case 
is demonstrated by Fig. 10.4. This figure represents the same infor- 
mation as Fig. 10.3 except the number of orbits has been extended from 
one to eight. The altitude at the completion of each orbit is plotted 
in Fig. 10.4; the sinusoidal oscillation which is seen during each orbit 
has not been represented in this graph. The analytical solution was 
carried through the first four orbits, at which time the atmospheric 
density was revised to conform with the altitude computed to exist at 
the completion of the fourth orbit. At the completion of two more 
orbits, the atmospheric density was again revised to agree with the 
altitude level derived at the completion of orbit number six. It may 
be seen that the analytical solution and numerical solution are 
essentially identical after eight complete orbits. Maximum difference 
between the two curves during the first eight orbits is less than two 
miles while the total altitude lost at the completion of eight orbits 
is approximately 27 miles. 

Equation (10-48) was derived with the assumption that an expon- 
ential model atmosphere exists. The only density value required in 
obtaining quantitative answers from this equation, however, is the 
initial atmospheric density eC; - Since the machine-computed numerical 
data for Figs. 10.1 through 10.7 assumed the ARDC model atmosphere, 
the values chosen for Ga for use in the analytical solutions plotted 
in these figures was also taken from the ARDC model. 

Fig. 10.5 shows a much more severe test of the accuracy of 


equation (10-48). The initial altitude of the vehicle is 80 miles above 


ft.* 


Cp 
the Earth and drag parameter is = 1.0 slug 


M . Fig. 10.5 shows that 
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Fig. 10.4 : Comparison of Analytical Solution to Machine Computed 
Numerical Solutions for 8 Orbits of High-Drag Zero- 


Lift Vehicle Initially in Circular Orbit at 120 Miles 
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pee 


the analytical solution and the numerical solution are essentially 
identical for the first 2200 seconds of the trajectory. The total time 
of flight computed numerically is in the vicinity of 6000 seconds; hence 
the analytic and numerical solutions agree for more than one-third of 
the total flight time. During this time, the total loss in altitude is 
approximately 30,000 ft. Also plotted in Fig. 10.5 are two curves 


(21) 


representing Nielsen's approximate solution of the trajectory. His 
approximate method involves solution of the differential equations of 
motions over suitably selected time intervals based on predictions of 
mean atmospheric density during this time interval. His approximate 
method.is represented in Fig. 10.5 in two separate curves, one for a long 
time interval and the other for a short time interval between piecewise 
continuous trajectory segments. The solution of this thesis shown in 
Fig. 10.5 follows the numerical trajectory solution more closely than 
Nielsen's approximate solution over the time interval plotted. A 
distinct advantage of the solution advanced in this thesis by equation 
(10-48) is the piecewise continuous solutions are not required for 
reasonable accuracy except after much larger time intervals than 
required by the @proximate methods of reference (21). 

Fig. 10.6 shows the accuracy of equation (10-48) when the atmos- 
pheric density ratio C5 is revised after 2000 seconds for the entry 
trajectory previously considered in Fig. 10.4. The one-step solution 
of Fig. 10.5 is shown in this figure for comparison. The curve of 
Fig. 10.6 was plotted for the first 2000 seconds by using equation 
(10-48) with Cs corresponding to 80 miles altitude. This solution 
was extended to 3000 seconds in order to predict the mean altitude 


during the interval 2000 sec. * +t © 3000 sec. With the mean altitude 
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thus predicted from the initial solution, the second step was plotted 
over the second time interval. The simple analytical solution and the 
machine-computed numerical solution are clearly very close for the first 
3200 seconds of the trajectory. During this time the vehicle has com- 
pleted more than half of its total flight time and lost approximately 
one-fourth of its total initial altitude. 

The method of determining values of Ar during the time interval 
greater than t = 2000 seconds in Fig. 10.6 requires some explanation. 
After carrying forth the original solution to 3000 sec., the average 
altitude during the interval 2000-3000 sec. was derived. (4, was 
determined for this altitude from the ARDC model atmosphere and the 
constant B, in equation (10-48) was adjusted accordingly’. It should 
be noted that in the derivation of equation (10-48), the initial altitude 
rate Ar' was assumed zero at zero time. At time 2000 sec., Which is 
zero time for the second step, M4r' is not zero; on differentiating 
equation (10-48), it is seen that A4r' is in the neighborhood of 
~30 ft./sec. at time 2000 seconds. This value of Ar' must be 
carried through the second step of the solution. For example: 

(1) At t = 2500 sec., equation (10-48) is used to determine Ar 
by substituting 7 equivalent to 500 seconds into this 
equation (i.e., t - by = t, = 2500 - 2000 = 500 sec., where 
wii is initial time of the second segment and t2 is time 
measured in the second segment of theslution.) This gives 
the value of Ar at 2500 seconds if and only if Ar' were 
zero at 2000 seconds. 

* A. = 0 in equation (10-48) in this case since lift is zero. 


a i, 





(2) Sieg Am 'i0) = 3 30 ft./sec. for the second phase of the 
solution, approximately 30 X 500 = 15,000 ft. must be added 
in the negative direction to the Ar computed in (1) above 
in order to determine the total Ar at t = 2500 sec. This 


may be summarized as: 
ore) = Ar(t5) een, ke to (10-50) 
2 


Eq. (10-48) with 
T corresponding to 
to =t- tis 
The time variation of Ar' is plotted in Fig. 10.7. This is 


easily determined by differentiating equation (10-48): 


Ar'(T) = ae sinw,T - ze (1 - coswT) (10-51) 


It should be noted that the value of Ar' at 2000 seconds was carried 
through as an initial value for the second step in the two-step solution. 
It may be seen from this figure that the difference between the computer 
solution and the analytical two-step solution of this thesis differ by 
less than 5 ft./sec. over the first 3000 seconds of flight. During the 
first 1500 seconds of the trajectory, the numerical and analytical 
solutions agree exactly. 


Velocity may be determined from equations (10-16), (10-48), and 
(10-49) as follows: 


o 
; a 
Aa = Mt Ae [1 cos ur|+ ee Be Lr sinuk _ Oy [ear 
O 








ia Lue We Ze 
(10-52) 
Using equation (10-18) and equation (10-48) without the oscillatory 


terms, the last quantity in the above equation may be integrated. The 
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Z£RB (10-53) 
Kinetic-potential Velocity loss due to 
energy trade-off. atmospheric drag. 


Range rate may be determined from equations (10-6), (10-8) and (10-9): 


KI< 
< 


== (1 + SV ous 


~.= (10-54) 
i 


a. 


a ~ 


In equation (10-54), Ar(T) is given by equation (10-48) and Av(T) 
by equation (10-53): 


AB 
/ WU: t Cc 
Ket ae 7 Br | — 1atk?) =A. (iE costtt)|- Caiek fe es = (10-55) 





This is integrated to give range as a function of time: 


Z. ¢ Cc , A y #BT 
Xn = au Be (cos WkT— t)+ Acsinuc| a eer (c 7 1 


Ve Zar:Ac mC . 
+e ta 2) + BL 








(10-56) 

The total range to impact determined by numerical ethods 
forta zero-lift vehicle with = “6 ft.2-slug™> initially at 80 
miles in a circular orbit above the Earth is approximately 26,000 miles. 
Time of flight was determined to be 5991 seconds. The total range of 
the same vehicle computed from equation (56) in a one-step solution of 
the trajectory is approximately 28,000 miles. 

Equation (7-19) showed that the convective heating rate at the 


stagnation point is: 





, =-kh 
(Q.)5 = 18,000 (HFho 3. 2 BTU (10-57) 
JR ft.<-sec. 
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where (HF) 50 is the heating function ratio tabulated in Table 7.3 for 
the terrestrial planets and @ is the radius of curvature at the 
Stagnation point. 


Equation (10-57) may be written as follows: 


-khi. -k Ar 
2 2 


~ 18,000(HF) 9 vs (1 + 288 a é (10-58) 
EL 


(Q,) - 
Qo)5 JR 
Substituting for Av(T) from equation (10-53) and Ar(V) from 
equation (10-48) gives analytical expression for convective heating 
rate as a function of time. 


Stagnation point temperature in degrees Rankine as a function of 


time is represented analytically in similar fashion from equation (7-21): 
-khi -kAr 


\ de a. 3a - = © % 
Ts = eos X TO (HF) 59 (VF) V5 (1+ ~~ e (10-59) 
Vs ; 


An approximate time history of vehicle accelerations is given by 


differentiating equation (10-53). 


10.4 The Degenerate Elliptical Entry Trajectory 
The initial point for the analysis of the elliptical trajectory 


is at the first perigee point. An elliptical trajectory results if 

the velocity of the vehicle at its first perigeal passage is greater 
than that required to generate a @ircular orbit; i.e., the first perigee 
is characterized by an energy excess over circular orbital energy. In 
the specific case that a vehicle is launched from the planet for the 
purpose of orbitting the planet a certain desired number of orbits and 
re-entering, then an elliptical orbit results if the velocity at engine 


cut-off in the launch phase exceeds circular orbital velocity. In the 


(7 


case of interplanetary operations in which the high energy level of the 
vehicle in the interplanetary transfer ellipse is to be reduced by 
atmospheric braking, the degenerate orbital profile follows the first 
perigeal passage if the energy transfer is sufficient during this 
passage for the planet to capture the vehicle. 

The basic equation that must be solved in order to determine the 
time behavior of Ar was given as equation (10-17). It was shown in 
section 10.3 that the approximation given by equation (10-49) may be 
made in this equation. Therefore equation (10-17) is written for the 
exponential atmospheric model as: 
artecdor = © [fel] Seo 


i 


te 


+ 
2 
SE (Lt Sn EO) [eM ar 
O 
(10-59) 

The angular velocity term, W, in equation (10-59) represents a 
Slowly varying quantity during the circularization phase of the 
degenerate orbital profile. For the purpose of this analysis, W is 
assumed constant; however, by writing it in terms of vehicular velocity 
and the distance of the vehicle from the planet center, then W may be 
changed at each separate initial point chosen for piecewise continuous 


solution of the trajectory. Ww is defined as follows: 


w= £0 


(1.0260) 
Ty 


where Ny is the dimensionless period of the orbit. Using equations 


(9-8) and (9-12), the average angular velocity W is written in terms 
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of velocity and radius at the initial point as follows: 


5 

uw = & - V “ (10-61) 
rs 1 

It is possible to make considerable simplification in equation (10-59) 

ifm 


-kAr 


He 


J 


Equation (10-62) is automatically satisfied for zero-lift vehicles. 


[a +(L/DI\(C /2) r,e 1.0 (10-62) 


This approximation is accurate if: 


N< 1.0 
-kAr 
where N =(L/D\C/2) r;e (10-63) 
= 
The quantity L Cps is plotted as an approximation parameter in Fig. 
2D M 


10.8 versus the minimum altitude for which the elliptical solution 
developed herein is accurate in the case of lifting vehicles entering 
the Earth's atmosphere. 

As an example of the use of Fig. 10.8, consider a vehicle with the 


following characteristics’: 


L/D» = a5 
Che 
= = 0.4625 ft.*/slug 


N = 0.1 for this vehicle at an altitude of 54.5 miles above the surface 
of the Earth. By the time the vehicle has descended to this altitude, 
orbital flight conditions no longer exist, consequently the limitations 


on the solution as a result of assuming approximation (10-62) are 


* This particular vehicle is used later in this section to compare the 
analytical solution derived herein to machine-computed numerical solutions. 
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generally not severe. 


With approximation (10-62), equation (10-59) reduces to: 


Y 
Ar" + we Ar=A. eK4Ti ae A, . eK Argr (10-64) 
1 Z J 


where the constants are defined as follows: 


2 e 
A, = (L/D) Cs (10-65) 
V : V, 
es eR Cue 
A, = = [a a (——) (10-66) 
i 1 
, 
A= - (10-67) 
i 


In accordance with the method of solution outlined in section 10.2, a 
first estimate is made for r(T) in equation (10-64) by assuming a 
vacuum trajectory. With (; = 0, then C = 0 and equation (10-64) 


reduces to: 


Ar" + oe ae = Ap (10-68) 
The solution to equation (10-68), with Ar(0O) = 0 and Ar'(0) = 0 is: 


Ag 
Ar = 2 (1 - cos WT) (10-69) 


At apogee of the vacuum trajectory: 


Arg = 2 (ay - r;) (10-70) 
where a is dimensionless semi-major axis of the ellipse. Comparing 
equations (10-70) and (10-69) at apogee (WT = TY ) shows that: 


Ao = WF Gay = ap (10-71) 
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Therefore, by substituting for a, from equation (9-8) A, may be written 


in terms of initial velocity, radius, and angular velocity as follows: 


Zo 


2 

Veet: =)! 

A, = Pr [eta] Congas 
zz = we, rs 


With A, defined by equation (10-72) and the first estimate of 


Ar(T) given by equation (10-69), equation (10-64) is written: 


7 
Pee Sur Ave A, = Gh, eet atcoswT _ A ent qf (ehtSOS4T) Gy (10-73) 
O 


where A, is defined as follows: 


4 
kA5 v.é a 
we 2 =k ry ( ar (10-744) 
2 = Vs rs 


Equation (10-73) is the fundamental differential equation of the 
degenerate elliptical entry trajectory. 


Taking the Laplace transform of equation (10-73) gives: 


1 Au 
Ris) = Gee) [Ae + @ (sA, ~ Az) Fs) (10-75) 


where R(s) is given by equation (10-30) and F(s) is the Laplace transform 


t fh cos WT 


O Before the inverse transform of this equation can be 


taken, F(s) must be determined. A search of the mathematical literature 


failed to disclose a representation of eco a 


suitable for this 
purpose. Derivation Summary 10.1 analyzes this function in terms of 
modified Bessel Functiors of the first kind. This investigation led 
to the following representation: 


coclUlUmel ellCmDESS!lUlUlUllUrWDSsSUwU COE eS lL Se See Se Sie i Ea e—s 


* Note thet Ar(0) = 0 and Ar'(0) = 0. 
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co 
WY 
ett COS T (Ay) + 2 > I (Ay) cos nwt (10-76) 
hal 


where I, is the nin Bessel Function of the first kind with pure 


imaginary argument. 

Fig. 10.9 is a graph of I for small values of Ay taken from 
reference (80). Engineering analysis of most entry trajectories leads 
to fairly large values of Ay; for example, A), for the vehicle plotted 
in Fig. 10.10 is in the neighborhood of 58 for the first segment and 
19 for the second segment. The following asymptotic aoe?) is 


convenient for large values of Ap : 


TA) oh Gro B), Gre-P)(4n*- 3°) 
(= eae iaa, * BIA | 


Ie 


(10-77) 


The magnitude of I, for small n is generally extremely large (e.g., 
of the order of 10!" or greater) due to the eT usm in equation 
(10-77). 


The Laplace transform of equation (10-76) is: 


es 


dig 
ae). = = c 2s) s* 4 (nw)* (10-78) 


n=1 


Substituting equation (10-77) into equation (10-75), expanding by 


partial fractions, and taking the inverse transform of the result gives: 


Brcr) = (Ae a Aa be™*)(1- COS wr) 


UJ)” 


ie 





~A , Y : 
Se [As] lr 882) +, Lor ssrA], (S08 —reosu) 


‘yi a A loos wT — Cos nwt] = As [sin wT — sau) 


Ve 2 





(10-79) 
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Equation (10-79) is the solution for altitude as a function of time 
for the degenerate orbital profile. Following is a summary of constants 


in this equation: 


I, : Fig. 10.9 for small A); equation (10-77) for large A). 
uw: Equation (10-61) 
A, : Equation (10-65) 
A, : Equation (10-72) 
Az: Equation (10-67) 
Ay : Equation (10-74) 


s 


The term in curly brackets appears at first glance to be formidable 


for engineering computations. Fig. 10.9 shows, however, that this 


Zikin ) 


n*—4 





solution converges rapidly due to the fact that terms involving ( 
become relatively unimportant for n greater than 2 or 3. For all 
reasonable initial problem conditions the asymptotic equation (10-77) 
is an accurate and easy method for computing i's It should be noted 
that the exponential term in equation (10-77) eantels the exponential 
term in equation (10-79). 


At perigee: 


sin nwT = 0 
(10-80) 
cos nwlT= 1 
Equation (10-79) reduces to the following for the time history of 


perigeal altitude: 


_ As -A 
ae =e > ee fel [)v (10-81) 
At apogee: sin nWT = 0 
(1.0=82) 
cos nw? = -1 
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Equation (10-79) reduces to the following for apogeal altitude as a 


function of time: 


At, (1) = = (A, eee ll, es) _ He r(Io+ Lie” 


uj* 





(10-83) 


m 


Equations (10-81) and (10-83) were plotted in Fig. 10.10 for a 
lifting vehicle launched in an elliptical orbit with initial velocity 
3% greater than circular orbital velocity. Initial altitude was 
350,000 feet above the Earth and the vehicle parameters were arbit- 


rarily selected as follows: 


ea. 5 
Cas 
— = 0.4625 ft.“ /siee 


The analytical solution is represented in two segments with a revised 
atmospheric density Cs chosen for the second segment; the first segment 
covered 103 orbits and the second segment used the atmospheric density 
corresponding to the altitude computed for the llth perigeal passage. 
Total flight time during circularization was computed to be 20.3 hours. 
The perigeal decay rate was computed from equation (10-81) to be 
approximately 0.14 ft./sec. during the first segment and 0.61 ft./sec. 
in the second segment. Perigee decayed 8030 ft. over the first 960 
minutes and an additional 8780 ft. between 960 and 1200 minutes. 

Total loss in perigeal altitude during the 20 hour flight was about 3 
miles. Apogeal decay rate averaged 31 to 32 ft./sec. during the first 
segment and increased to 46 ft./sec. during the second segment. Apogee 
decayed more than 360 miles during the first segment and 170 miles 


during the second segment. 
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Fig. 10.10 : Comparison of Perigeal and Apogeal Altitudes of 
Analytical Solution with Fhose Determined Numer- 
ically with IBM-704. Lifting Vehicle Launched 
with Initial Velocity 3% Greater than Circular 
Orbital Velocity 
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In order to estimate the accuracy of the analytical solution 
discussed in the preceding paragraph, Fig. 10.10 also shows decay of 
perigee and apogee computed numerically with the IBM-704. The numerical 
runs represented a trajectory initially inclihed 30° to the equatorial 
plane about the oblate, pear-shaped Earth with an exponential atmos- 
pheric model. Atmospheric density was assumed zero above 600,000 ft. 
in the computer run. Total flight time in the circularization phase 
was computed to be about 19.9 hours. 

Some rather interesting facts may be deduced from equations (10-81) 
ena (10283). First, it may bevseen that lift has little effect cu 
perigeal and apogeal altitudes during the circularization phase. 
Equation (10-81) shows perigee decay rate to be independent of lift. 

The only term involving lift in equation (10-83) is the term containing 
Al. For the particular example plotted in Fig. 10.10, the Ay term was 
orders of magnitude smaller than the A> term. The apogeal decay rate 
given by the second term in equation (10-83) is independent of lift. 

The major difference between the decay rate at apogee and that at 
perigee is the fact that perigeal decay rate is proportional to the 
difference of two near-equal large numbers (I, - I,), while apogeal 
decay rate is proportional to the sum of these two large numbers. It 
may be seen from Fig. 10.9 that the difference between these two numbers 
is essentially constant for Ay > 5; the sum is very sensitive to the 
magnitude of Ay. 

For most degenerate orbital trajectories, Ay is large enough 
that the first two terms in equation (10-77) are sufficient for accurate 
estimation of perigeal and apogeal decay rates. Differentiating 


equations (10-81) and (10-83), and substituting the first two terms 
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from equation (10-77) for I, and 1, gives: 


} _ As 
Alt = — lew of AE _ 


e ae 


lle 


Other guidance parameters may be computed from the foregoing 
results. Some of the quantities of interest are summarized below: 


(1) Dimensionless semi-major axis ay: 


ly + la 


eee? cee | ea 10-86 


Using equation (10-81) and (10-83), equation (10-86) becomes: 


Qn 


vy 
Man + (Ata) bs As eh To (10-87) 


The semi-major axis changes almost as a step function in the vicinity 
of perigee; hence equation (10-87) should not be viewed as a continuous 
time solution. Recomputation of ay should be performed at each perigeal 
passage with a near constant value during the following orbit. 


(2) Period of orbit: 


Y 
T, = 2 a (10-88) 


where ay,(T) is given by equation (10-87). 
(3) Average angular velocity in orbit: 
-% 
WT) = ay (10-89) 


where a,(T) is given by equation (10-87). 


ol 


(4) Instantaneous velocity: 


In a Keplerian elliptical orbit, velocity is given by: 


ae 
Vv = 


(10-90) 


ry tO 
1 

a 
g 


N 
Velocity may be approximated during the circularization phase from 
equation (10-90) as follows: 


ey ae ary) i 
eS (10-91) 


Ar(T) is given by equation (10-79) and ay (T) by equation (10-87). 
Velocity at perigee is derived from equation (10-91) as follows: 
ae = 11) — (Ar, -3ar) alomey>) 
1 G Ae : 
where Ar, (7) and Ary (T) are given by equations (1026) 
and(10-83). 


Velocity at apogee is as follows: 
Ze en SAT 
oS = Ly A i =) (10-93) 


(5) Total energy: 
Total energy changes approximately as a step function 
in the vicinity of each perigeal passage. The time behavior of energy 


level is given by: 


aL 


E a 
(tot) Zan 


(10-94) 
where ay(T) is given by equation (10-87). 


(6) Eccentricity 


Eccentricity, like the semi-major axis and energy level, 
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decreases in a series of near steps in the vicinity of perigee. The 
value between perigeal steps is reasonably constant. With this step 
behavior recognized, the time behavior of eccentricity is given by 
equation (10-94) (which was derived from equations (9-16), (9-18), 


(10-92) and (10-93)): 
San = = (At, — At) (10-95) 


(7) Heating Rates, Stagnation Temperatures, and Accelerations: 


During the circularization phase of the degenerate orbital 
profile, heating and accelerations are not of major concern in guidance 
analysis since maximum levels are generally encountered after the 
trajectory has degenerated to a circle. The methods of Chapter 9 and 
Section 10.3 may be used to predict these quantities during the Gas- 


Dynamic Phase of flight. 


10.5 Summary 
Chapter 10 examines the trajectory of a vehicle in the degenerate 


orbital profile. Initial conditions assumed in this analysis were that 
a lifting or non-lifting vehicle was injected in horizontal flight at 
edi’ initial altitude above an arbitrary planet with A4r(0) and Ar'(0) 
equal to zero. The results were described in dimensionless form suitable 
for analysis of entry into the atmosphere of any planet. The exponential 
atmospheric density model was assumed. 
Two separate trajectory phases were examined: 
(1) The circularization phase for the degenerate orbit; 
(2) The entry phase which follows the circularization process. 


In the special case wherethe vehicle is initially at circular orbital 
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velocity, the first of these two phases obviously does not exist, 

The solution for altitude during the circularization process was 
derived as equation (10-79). This solution was written in terms of a 
number of constants which depend on the initial conditions of the 
problem and on the characteristics of the vehicle. Some of the important 
constants in this equation were written in terms of Bessel Functions. 
Since these functions are tabulated, numerical answers to the resulting 
equations are easily obtained. The time variation of perigeal and apogeal 
altitudes was given in equations (10-81) and (10-83) and the rate of 
decay for most practical problems is described by equations (10-84) and 
(10-85). Analytical representation of other quantities which are important 
in the conceptual and preliminary design stages of guidance systems are 
summarized in equations (10-86) through (10-94). 

It was shown in the analysis of the circularization phase that 
the drag characteristics of the vehicle are important in specifying the 
resulting trajectory and that the lift characteristics of the vehicle 
are relatively unimportant. It was shown that, within the approximations 
made in this analysis, lift does not enter in the specification of the 
decay rates of apogee and perigee. Perigeal altitude, to a first order, 
remains essentially constant. The major difference between the decay 
rates of apogee and perigee is that perigeal decay rate is proportional 
to the difference of two near-equal large numbers” while apogeal decay 
rate is proportional to the sum of these same numbers. 

The solution for altitude in circular orbital entry was given by 


* Perigeal decay is proportional to (Ip - I,), where I. and J, are the 
¢ * * e QO } 
first two Bessel Functions with an argument that is generally greater 
than 10, 
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equation (10-41). This solution was written in simplified forms suited 
to almost all practical problems as equations (10-47) and (10-48). It 
was shown in this chapter that a true circular orbit or a linear 
decaying circular orbit cannot exist even under idealized conditions of 
injecting a vehicle exactly at circular orbital velocity near a spherical 
planet. The influence of the atmosphere on the dynamics of energy 
transfer result in undamped oscillatory motion in altitude and flight 
path angle. 

The analytical solutions derived in this chapter were compared to 
machine-computed numerical solutions under a variety of initial 
conditions. It was shown that the analytical solutions presented an 
accurate picture of the resulting motion except after long periods of 
time. The analytical solutions depart from more accurate computer 
solutions cnly after there is a significant increase in atmospheric 
density from the assumed initial value as a result of altitude loss 
during entry. This limitation on the solutions may be greatly 
alleviated by starting the solution over again under a new set of 


initial conditions. Examples of this method of solution were given. 
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Derivation Summary 10.1 


a CosX sn Terms of Modified Bessel Functions of the 


™ 


Expansion of e 


First Kind. 


In order to solve equation (10-75), the Laplace transform of 
a cos xX . ; 
& is required. An expansion of this function in terms of 
Bessel Functions of the first kind with pure imaginary arguments was 
carried out in this thesis. This particular expansion is not given in 


the mathematical literature that the author has searched. 


The following is written: 


2 2 Z 3 
QCos x a“ cos x Q’ cos~ x 
e — socom of ll i ee (1) 


th power of cos x: 


c nx +n cos{(n-2)x] + ntn=t cos|(n-A)x) 


= cine cos|(n-6)x| oa ooo } 


The following equation defines the n 


rl 
COS Xx = 


1 
Fe 











(2) 
Substituting equation (2) into equation (1) gives: 
2 3 
pe _ 1 oes a Fi oes ax a 1) 4 = s- (08 5x + 3cos x) 
A 5 

a a eae, Gael 
1D dh eosdnrdeas2e+$) 4.2 (cosSx- Scosdx + lose 
af ial 4 (cos 6x + 6cos 4x + I15cos 2x + =) 

6! 32 fs 

a 4 

i eA co? fx * [cos 5x + 2icos Sx + 3) c0s x) + o00 


(3) 
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Derivation Summary 10.1 (cont.) 





Equation (3) is written: 


Qcge% lie f o te + 7 I, cos 2x as 2» T COS NX (4) 
n=3 


After a rather difficult process of mathematical bookkeeping is 
performed on equation (3), the coefficients in equation (4) are derived 
as follows: 


ee 


Cy 
Equation (5) is the series representation of Bessel Functions of the 
first kind with~pure nowthady eummange?” - 
~n 
aj=i J,(ia) 


where D are Bessel Functions of the first kind. 
I (a) are plotted in Fig. 10.9 for small values of a. For large 


a, the asymptotic series given as equation (10-77) is adequate. 
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Appendix A 
COORDINATE FRAMES USED IN ENTRY TRAJECTORY ANALYSIS; 


GLOSSARY OF SYMBOLS, CONSTANTS, AND DEFINITIONS 


A.1 Coordinate Frames 
The guidance process consists of measurements of vehicle position 
and velocity, computation of control actions necessary to properly 
adjust position and velocity, and delivery of suitable adjustment 
commands to the vehicle's control system. Guidance of vehicles through 
or above the atmosphere of a planet involves the twofold functions of 
indication and control. Navigation is the indication aspect of guidance. 
Inherent with an investigation of guidance requirements of vehicles 
entering planetary atmospheres is the selection of suitable coordinate 
frames of reference”* and trajectory parameters. Many sets of coordinate 
frames are required to completely specify the three dimensional nature 
of the entry trajectory under the following general conditions: 
(1) The figure of the planet is non-spherical; i.e., it may 
be an oblate spheroid with local surface irregularities 
and large scale harmonics leading to a figure resembling 


that of a pear. 


* In this thesis, "frame of reference", "frame", "reference space", 
etc., are used synonymously. 


oo? 


(2) 


©) 


(4) 


(5) 


(6) 


(7) 


The planet is rotating about an axis through its center 
ef gravity. 

The gravitational mass attraction of the planet is not 
spherically symmetric. . 

The atmosphere surrounding the planet is rotating with 
the planet. Relative motion between the atmosphere and 
planet may exist; i.e., winds may be present. 

The entry vehicle may have the capability of generating 
variable lift, drag, and thrust forces in order to control 
the shape of its trajectory. 

The vehicle may have the ability to bank; i.e., rotate 
the direction of the lift vector out of the plane of 

the trajectory. 

The propulsive system may be mounted on gimbals in order 


to vary the thrust direction as desired. 


Seven sets of coordinate systems were chosen in order to analyze 


the three dimensional dynamical equations important in the study of 


guidance. One of these coordinate systems is the guidance grid, for 


which two alternatives are presented herein. The coordinate systems 


selected were as follows: 


ae) 


(2) 


(3) 


Planet centered Inertial Frame: Newton's laws are valid 
in an inertial frame. The use of a planet centered 
inertial frame is justified because gravitation and 
acceleration effects are indistinguishable. 

Planet Reference Frame: Establishes the orientation of 


the planet with respect to the inertial frame. 


Instantaneous Trajectory Plane Frame: Establishes the 
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position and orientation of the instantaneous plane of 
the trajectory with respect to the inertial frame, 

(4) Position Reference Frame (or Guidance Grid): ‘Two alter- 
native methods discussed in subsequent paragraphs of 
Appendix A are: 

(a) The Latitude-Longitude Triad 
(b) The Great Circle Triad. 

(5) The Wings-Level Triad: This frame is basically orientated 
by the velocity vector of the vehicle with respect to the 
atmosphere when the lift vector is in the instantaneous 
plane of the trajectory. 

(6) The Vehicle Coordinate Triad: Required for a vehicle that 
has the ability to "bank"; i.e., to rotate the lift 
vector out of the plane of the trajectory in order to 
generate trajectory curvature as desired. 

(7) The Engine Gimbal Triad: Required for a vehicle that has 
an engine capable of being rotated with respect to the 
vehicle in order to control the direction in which thrust 
forces are generated. 

Relations among the angles involved are presented for ready reference 
in converting quantities specified in one frame into components in any 
of the other frames. 

Fig. A.l shows a vehicle located at the point P entering the atmos- 
phere of a planet. The destination of the vehicle is some landing 
Site on the surface of the planet that was selected prior to initiation 
of entry. The instantaneous plane of the trajectory is not fixed with 


respect to an inertial reference. Following is a brief definition of 


Ol 


Instantaneous Plane of 
Trajectory 





Equator 


Line of Nodes / 
Xo 


Planet Centered Inertial Frame: Centered at the center 
of the planet and sidereally nonrotating relative to the 
"fixed stars". Xr and Y,; are in the equatorial plane 
and Z, is directed along the polar axis (North). 


Planet Reference Frame: Centered at the center of the 
planet and nonrotating with respect to the planet. This 
frame rotates about the planet's polar axis relative to 

the inertial frame at the planet's daily sidereal rate. 

X54 and Yo are in the equatorial plane and fixed with 
respect to the planet. Zo is directed along the polar axis. 


Instantaneous Trajectory Plane Frame: Centered at the 
center of the planet. X, is the line of intersection of 


equatorial plane with the instantaneous plane of the 
trajectory; Y+7 is in equatorial plane perpendicular to 
X- ; 4; is directed along the polar axis. 


Fig. A.1l: Definition of the Inertial, Planet, and Instantaneous 


Trajectory Plane Coordinate Systems. 
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the coordinate frames shown in Fig. A.l. 


Xr Vy Zy Planet Centered Inertial Frame: Centered at the center 
of the planet and sidereally nonrotating relative to the 
"fixed stars", Xy and Yy are oriented in the equatorial 
plane and sidereally nonrotating; ay is directed along 
the polar axis (North). 

Kor Yor 2 Planet Reference Frame: Centered at the center of the 


planet and nonrotating with respect to the planet. This 
frame rotates about the planet's polar axis relative to 
the inertial frame at the planet's daily sidereal rate. 
Xo and Yo are oriented in the equatorial plane and fixed 
in the planet; Z. is directed along the planet's polar 


O 
axis (North). 


Xn Los Zn Instantaneous Trajectory Plane Frame: Centered at the 


center of the planet. Xn is the line of intersection of 
the planet's equatorial plane with the plane of the 
trajectory (ascending line of nodes); Yp axis is in the 
planet's equatorial plane and is perpendicular to the 
ascending line of nodes. Ln axis coincides with the 
planet's polar axis (North). 

In an analysis of the trajectory of entry vehicles, it is conven- 
ient to study the motion of the vehicle in spherical polar coordinates. 
The reference frame is usually selected with one axis along the radius 
vector from the center of the planet to the point at which the guidance 


is taking place. This radius vector, for the planet Earth, is the 


geocentric radius. In this thesis, the term "geocentric" is used in 
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a more general sense to identify quantities associated with the spheri- 
cal polar coordinates of any planet”. 


The following coordinate axes systems are shown in Fig. A.2: 


lal slip 7 Geocentric Latitude-Longitude Triad: Orthogonal set of 
unit vectors centered at the entry vehicle center of 
gravity and directed in radial, longitude, and latitude 


directions respectively. 


‘laps 1g Ly Instantaneous Great Circle Triad: Orthogonal set of unit 
vectors centered at the vehicle center of gravity and 
directed in radial, along-track, and cross-track direc- 
tions respectively. 

The following quantities are defined in Fig. A.2: 

Line of Nodes: The line of intersection of the plane of the instan- 
taneous trajectory with the equatorial plane. 

Ascending Node: The point of intersection of the line of nodes 
and the northerly-directed segment of the instantaneous 
ellipse describing the trajectory. 

Descending Node: The point of intersection of the line of nodes and 


the southerly-directed segment of the instantaneous ellipse 
describing the trajectory. 

R : The position vector directed from the center of the planet 
to the instantaneous position of the vehicle. 


* The interpretation of "geocentric" in this thesis may be considered 
to be a shortening of "geometric-centric" rather than the more precise 
“geoid-centric", 


HO 





re 
i 6 


Z. la (North) sc 


lg (Along Track) 










| 


(Across Track) 


Plane of 
Trajectory 


ais 
S ee 


Equatorial Plane 





Iepaios VAN Geocentric Latitude-Longitude Triad: Orthogonal 
set of unit vectors centered at the entry vehicle 


center of mass and directed in radial, longitude, 
and latitude directions respectively. 


l-,'g,ly Instantaneous Great Circle Triad: Orthogonal set 
of unit vectors centered at the vehicle center of 
mass and directed in radial, along-track, and 
across-track directions respectively. 


Fig. A.2: Definition of Latitude-Longitude and Great 
Circle Guidance Grids. 
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I? The instantaneous velocity vector of the vehicle with 

respect to inertial coordinates. 

TS; Inclination of the instantaneous trajectory plane with 
the equatorial plane (measured: from equator to trajectory 
plane). 

Art : "Inertial" geocentric longitude of ascending line of 
nodes; i.e., angle measured in the equatorial plane from 
Xr iO Kne 

App : "Inertial" geocentric longitude of the vehicle; i.e., 

angle measured in equatorial plane from X,. to entry 


vehicle P, 


A: (not shown in Fig. A.2) geocentric longitude; i.e., angle 
measured in equatorial plane from Xo to the vehicle. 
A: geocentric latitude of the vehicle. 

ie (not shown in Fig. A.2) geographic latitude of entry 
vehicle. 

QO : angle measured in the plane of trajectory (in direction 
of vehicle motion) from ascending line of nodes to the 
vehicle. 

ba bearing angle of entry vehicle as seen by an observer 


directly under the vehicle on non-rotating planet; angle 
between North ( 1, ) and the horizontal component of 


velocity (directed along 1g ye 


The following angular relations are useful in conversion among the 


sets of angles defined in Fig. A.2. 
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Zable A-1 


Angular relations among coordinate systems given in Fig. 2% 


sin A sin @ sin 


sin( Arp - Arr) sin 9 cosV_ 


sin A 


COS 


cos( Azp - Aq) ~ CosA 


sinW cos ¢ cos/\ cos 6 


cos A sin 6 


tan A 


sin Ayp - Agr! 
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The drag vector is directed anti-parallel to the velocity vector of 
the vehicle with respect to the atmosphere. The lift vector is perpen- 
dicular to this velocity vector and acts in the plane of symmetry of 
the vehicle. In the analysis performed herein, it is assumed that the 
entry vehicle has the capability of rotating the lift vector (i.e., 
banking) in order to produce curvature of the trajectory in a controlled 
manner to enable the vehicle to reach landing sites at some distance from 
the "no-bank" trajectory. For military vehicles, such as a reconnaissance 
platform, the banking capability may permit adequate military coverage 
of an entire hositle nation from a single parent satellite. 

Two additional vehicle-centered coordinate triads are defined in 
Fig. A.3; the "Wings Level" triad a Ty? be 2m the "Vehicle" 
triad 1,, ly: 1,- it is noted that both the "Wings Level" triad and the 
"Vehicle" triad are basically oriented by the velocity vector of the 
vehicle with respect to the atmosphere. The left wing unit vector in 
the zero bank condition yo always lies in the 1, = i plane . 
Rotation of the vehicle about Vv am) (coincident with Ty? 1,.) causes 
the lift vector to move out of the vertical plane, hence causes a side 
force which tends to curve the trajectory. The amount of rotation of 
the vehicle from the wings level condition is defined as the bank angle, B. 

Following is a brief description of the coordinate systems shown 


merie. &. 3: 


1 , iE Wings Level Triad: Orthogonal set of unit vectors 
uu Fu’ 


centered at entry vehicle center of gravity. 1 


nae allem Dell ell Dell OllmlmlCCLDOSSP SSS eee Sl 


* 1, #£4«+9185 positive in the geocentric East direction. l,j is positive 
in the geocentric North direction. The la - 1, plane is the 
geocentric horizontal plane, or "level plane" -- hence the term 
"ings Level" coordinate system. 
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Bank Angle 







(along 
left wing) 


A.3: "Wings Level" and "Vehicle" Coordinate Triads. 
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is along the velocity vector of the vehicle with 


respect to the atmosphere. ly is along the left 


Wing of the entry vehicle in the zero bank condition. 


a is along the lift vettor in the zero bank 


“WL 
condition. 
Vehicle Coordinate Frame: Orthogonal set of unit 
vectors centered at entry vehicle center of gravity. 
a. is along the velocity vector of the vehicle with 
respect to the atmosphere. 1, is along the left 
wing of the vehicle; 1, is in the plane of symmetry 
of the vehicle and is a unit vector in the lift 


direction. 


It is noted that ik, and 1 are always coincident. l, differs 


XT, Jy 


from Lyin and 1, differs from Lot by the angle of bank, B. In zero 


bank condition, both triads are coincident. 


meal 


The following angles are defined in Fig. A.3: 


=A ( Laie) : angle between "East" ay ) and horizontal component 


of vehicle velocity with respect to the atmosphere. 
ae) angle between the horizontal component of vehicle 
1A; ~1, 

velocity with respect to the atmosphere (i.e., the component in 


the 1) - In plane) and the total velocity vector of the vehicle 


with respect to the atmosphere. 


5 Sa & = =A (= —\ : bank angle of vehicle. Corres- 
Tig “i Ton 5 Te) _ 
ponds to rotations about the 1, axis from the wings level condition. 
Following are definitions of A and Y in terms of velocity 
components: 
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A,, = arctan vam = are sin sth =are cos __, amir ‘ (A-1) 
a (Van + Varn 4 (Viarnn + Veen )* 
os 
¥ = are tan Mare = arc sin Vamte (Viawn + Viana) 
= —— 3 _ 
(Venn + Vim) mas arc co Van (A-2) 


Rotations from the 1, , =” 1, triad to the i, 1, i. triad are 
performed in the following order: 
(1) Rotate Gr, sly ae i, ) to GF; lA» Lys by angle A,, about 


= axis. 2 


(2) Rotate (1,, 1A, Lt) to (le Lae Lyn) by angle ¥ 


about —— axis. 


(3) Rotate Oy da 
aly axis. 


“WL 


hvu Rom 1) by angle B about 


Direction cosines between io ar : Th, triad and fi, sl ne triad 


are as follows: 


te Le ijing Sea 
x H H 
H 


1, 
4 — in A, cos B cos A cos 6 
1, cos? sinB 















-~ sin? cos AL sinB - sin? sin A. sin B 


(A-3) 


sin Ay sin B - cos A, sinB 








- sinYcos Ay cos 6 -~ sin? sin Ay cos B 
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The coordinate systems defined above are not those used in 
conventional aircraft analysis where it is common to define a set of 
axes fixed to the airframe. It is necessary to define a set of vehicle- 
fixed axes if the stability and control charaoteristics are to be 
studied -- a subject beyond the scope of this thesis. Lift and drag 
forces are the primary gas-dynamic forces affecting the shape of the 
trajectory in the atmospheric phases of flight; secondary forces and 
moments are of lower order. 

Thrust forces generated by the entry vehicle propulsion system 
will, in general, be applied in such a direction as to slow the vehicle 
down” (retor-thrust). It was assumed that the thrust vector passes at 
all times through the vehicle center of gravity in order to avoid 
excessive moments. The engine gimbal angles shown in Fig. A.4 were 
selected to describe the thrust direction. 


Rotation from the lige ty 
xXx V Z 


triad to the thrust triad (ys 


At 


yo? 12, ) is performed in the following order: 


(i) Router (ie; L_, a) to (Ly, » aA 


y ; 22) by angle A, about 


- 1 axis, 
J 
(2) Rotate Ads * ly,» 42? to Cx, ly, 12) by angle A, 


about 1, axis. 
ae 


Direction cosines between 1, ; a ne 
2 Z 


Yo" triad and ee dys LF triad 


are as follows: 


* Except possibly in the landing phase. 
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Fig. A.4: Engine Gimbal Angles 
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The following angles are shown in Fig. A-4: 


A : Engine gimbal angle generated by rotations about the 
negative a axis. 
A ; Engine gimbal angle generated by rotations about the 


displaced z axis (after rotating through angle Ae). 


Thrust components may be written in component form in the 1, 
ly ; a. triad by using the direction cosines given in (A-3). These 
components may, in turn, be transformed to the 1,, 1g, ly triad by 
the table of cosines and angular conversion relations given in Table A-l. 

Once the velocity requirements of the entry vehicle are established, 
the propulsion system must be designed in such a way that thrust may be 
applied in the proper direction with the proper magnitude and for the 
proper length of time; i.e., thrust vector control. Thrust vector 
control is also important for generation of orientation torques for 
attitude stability; problems connected with thrust vector control for 
attitude orientation and rotational moment stability are not considered 


in this thesis. 


Wily 


A.2: The Two-Dimensional Trajectory 


Most of the analytical study @scribed in this thesis was concerned 
with a two-dimensional approximation of the entry trajectory. Following 
is a brief summary of the principal assumptions made in postulating a 


two-dimensional trajectory: 


(1) The gravitational field of the planet is spherically 
symmetric and is free of anomalies. 

(2) The atmosphere rotates with the planet, and no wind 
exists. The velocity component of the rotating atmos- 
phere in the ly direction is negligible. 


(3) The vehicle flies in the zero bank condition with Ag = 0. 


Figure A.5 defines the important forces and geometric quantities 


of the two-dimensicnal trajectory. 


A.3: The Instantaneous Ellipse 


It is common in celestial mechanics to describe the motion of 
heavenly bodies in terms of parameters of conic sections. There are 
many definitions of conics; the following definitions are convenient: 

(1) Ellipse: The locus of points the sum of whose distances 

from two fixed points (foci) is constant. 

(2) Circle: The locus of points equidistant from a single point 

(special case of the ellipse). 

(3) Hyperbola: The locus of points the difference of whose 

distances from two fixed points (foci) is constant. 

(4) Parabola: The locus of points equally distant from a fixed 


point (focus) with a fixed straight line (directrix). 


N15 


| ii 
Initial Pomt 
(X,=0) / 
L 





Definitions : 


= distance flown 


<i > 
il 


velocity of vehicle with respect to coordinates rotating 
with the planet (same as velocity of vehicle with respect 
to the planetary atmosphere, Voam) » for no-wind condition) 


= altitude of vehicle above planet 


= flight path angle with respect to local geocentric 
horizon - positive as shown. 


mass flow rate of propellant - 


= equivalent exit velocity of rocket propellant (assumed 
constant) 


on = x 
il 


= instantaneous mass of vehicle 
drag 
= ee 


= engine gimbal angle of retro-rocket; angle the thrust 
vector makes with negative 1, axis. 


= & OR 
il 


Fig. A.5: The Two-Dimensional Trajectory 
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When a body is in motion under the action of an attractive central 
force that varies as the inverse square of the distance, and if no 
external forces act on the body, the path described will be a conic 
whose focus is at the center of attraction. A particle moving accor- 
ding to such a force obeys Kepler's laws. Kepler's three laws of 
planetary motion, published around 1610, were the result of his pioneering 
analysis of planetary observations, and laid the ground work for many 
of the important contributions of Sir Isaac Newton. Kepler's laws may 
be summarized as follows: 

A particle moving under the action of an attractive central 

force that varies as the inverse square of the distance 

(1) travels in an ellipse or a hyperbola (or their special 
cases, a circle or parabola) with the attracting center 
at one of the foci; 

(2) the radius vector from the center to the particle sweeps 
out equal areas in equal time ; 

(3) for the elliptic orbit, which results in periodic motions, 
the squares of the periods of rotation are proportional 
to the cubes of the major axis of the orbits. 

Klliptical parameters are useful for analyzing the motion of the 
entry vehicle during the orbital phase and in the early phases of entry 
where gas-dynamic forces are negligible** The use of elliptical 
parameters may be extended to the study of the entire trajectory by 
using the "instantaneous ellipse" technique; i.e., by describing the 
trajectory as a time-varying ellipse which matches the position and 
velocity vectors of the vehicle at each instant. 


* Kepler's second law, the cap eer tie of areal velocity, is a general 
theorem for central force motion Since angular momentum is always 
conserved. The first and third laws are restricted specifically to 
the inverse square law of force, 


** This segment of the trajectory is called the "Keplerian Phase" in 
this thesis. 


hi? 





Frequent reference is made in this thesis to elliptical parameters, 
particularly in the study of retro-rocket system requirements described 
in Chapter 9. A description of the parameters of the ellipse are given 


in Fig. A.6. . 
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Descending 
Node 


Line of 
Nodes” 


Apogee — 


Ascending as 
Node 


Line of 
Apsides 


Line of apsides: The major axis of the ellipse, connecting the 
apsides (apocenter and pericenter). 


2a 
2b : 
21 1 =e |) 


Major axis 
Minor axis 
Latus rectum 


Eccentricity = € (Always less than 1.0 for ellipse) 
€= /1 - (b/a)* 


True anomaly = 6 
Area of ellipse = 7fab 


n=9- 0 


Equation of ellipse: R= (aaa a) 
1 +€ cos@ 

3/2 

T = period of orbit: T= ££ a7” 


V FpMo 


Fig. A.6: The Instantaneous Ellipse. 
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A.4 Glossary of Symbols 


A fundamental difficulty in the presentation of any mathematical 
description of physical phenomena is that of choosing a complete and 
workable system of symbols and notation. Following is a summary of 


symbols and definitions used in this thesis. 


A.4.1 Mathematical Symbols 


( ) : a line above a quantity designates a vector quantity; e.g., R 





° d 
( ) : a dot above a quantity represents total time derivative {2 


(') : denotes total derivative with respect to dimensionless time T7T, 


i. ie) &(m)0 
i.e., —— Note: Y= /——“~—t 


v( ): del or gradient of the quantity in parenthesis. 

() X (): vector cross product of quantities in parenthesis. 

() * ( ): dot or scalar product of vector quentities in parenthesis. 
The magnitude of a vector quentity is represented oy the letter desig- 
nation of the quantity with the vector symbol removed. Thus, the 


magnitude of R is R. 


a ) 


Sx ¢ partial derivative of quantity in parenthesis with respect to x. 





o( ) : variation of the quantity in parenthesis. 
A( ) : incremental change of the quantity in the parenthesis from its 


initial value; e:c., Al =a ai 


a 
OO 
, a ), 2 Summation of the quantity in parenthesis. 
n=1 
log ( ) = log, () 
im ( ) = tog’ ( ) 
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A.4.2 Subscripts 


ow JA : Components along the ie P 15 ‘ 1 axes, 
r, 0,W : Components along the les ; 1g ; Ly axes, 


i : Denotes initial value of the quantity. 

f : Denotes final value of the quantity. 

O : The subscript 0 refers to the generalized planet 0. When 
specialized to the planet Earth, "O" is replaced by "BE". 

m : Denotes quantities at the "trajectory modification point". 
The trajectory modification point is that point in the 
original satellite orbit at which entry is initiated through 
generation of thrust forces. (See Fig. 9.3) 

I : Measured with respect to inertial coordinates. For example, 

u is the velocity of the vehicle with respect to inertial 

space, while V denotes velocity of the vehicle with respect 
to coordinates rotating with the planet. 

a: Value of the quantity at apogee; e.g., Ry is the radius 

at apogee of an elliptical trajectory. 
Ti: Value of the quantity at perigee; e.g., R, is the radius 


at perigee of an elliptical orbit. 


A.4.3 Symbols 


The author has endeavored to be consistent in the use of notation 
snd symbols in this thesis. For example, the symbol " represents the 
same quantity regardless of where it may appear in the text. 

As is well known, certain particular symbols are used within each 


scientific field to represent quantities frequently encountered. Examples 
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of this are Q for atmospheric density and Cp for drag coefficient 

of aerodynamics, © for "true anomaly" of celestial mechanics, "p" 

or "h" for angular momentum in mechanics or physics, etc. Where possible, 
such conventions are followed in this thesis. * Unfortunately in an analysis 
that spans many fields, one finds a conflict between the common symbols 

of one branch of science or engineering with those of another field. In 
such instances, the author had no alternative except to invent a symbol 
which he considered convenient: the choice may often be unsuited to the 
taste of the reader who is conversant with the more conventional 


representation. 
A.4.3-1 English Letters 


acceleration vector of entry vehicle ,ft./sec.“ 


=| 


A : Engine gimbal angle in 1, - 1, plane (Figs. A.4 and A.5) 


(Ae) opt ; Optimum engine gimbal angle (i.e., that angle for which 
range errors are minimized with respect to errors in 
operation of retro-rocket system). 

As ; Engine gimbal angle; consists of rotations out of the 


18 a i plane (Fig. A-4). 
A : Heading angle. Angle between "Rast" (1> ) and horizontal 
component of vehicle velocity with respect to the atmos- 
phere (Fig. A.3). 
A : Constant defined in circular trajectory study of section 
10.3 as follows: 
Vs 
A= Ges 


Ay : Constant defined for degenerate elliptical orbital profile 
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in section 10.4 as follows: 


2 
v 


+s 


Ay 


(I 
on 
vik 
nd 
R 

it 
i 


: Constant defined in section 10.4 as pllows: 


2 ‘ 1 
wr = 
me | era d 


: Constant defined in section 10.4 as follows: 


2 
—— via 
2 i 


s 


Constant defined in section 10.4 as follows: 


=i hha 


One-half the major axis of the elliptical trajectory 


: Dimensionless semi-major axis of the ellipse: 


—_ — 
N Rimyo 





a 


: Angle of bank of the entry vehicle; angular rotation about 
the velocity vector of the vehicle with respect to the 
atmosphere (Fig. A.3) 


Constant defined in section 10.3 as follows: 
¥) 
z 


One-half the minor axis of the ellipse. 


If 1 
B, = Cu (= 5 


oa 


: Dimensionless semi-minor axis of the ellipse 


be b/ROn)0 
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Constant defined in section 10.3 as follows: 


Constant defined in section 10.3 as follows: 
bs = x, = 2 yy 

Constant defined in section 10.3 as follows: 
ba | ates 5) 

Trajectory constant. Defined as follows: 


C= Rimyo ei See 


Constant defined in section 10.3 as follows: 


Dimensional constant for the planetary atmosphere 


(Btu tt.73/? sec.-}) 

Drag coefficient 

Zero-lift drag coefficient. 
Lift coefficient, 


An acceleration constant defined as follows: 


c= . Log Cro 
= “he leeie 


Auxiliary quantity defined as follows in Eq. (8-24) 


a ale caty Rimio 
a ——— 
. Re tog Ge 4 ( — Be) 
Lau 





ai kin): 
Cas 


© Cue) 


Bin) ° 


= (pot) 


ect) 


Auxiliary quantity defined as follows in Eq. (8-28): 


C= Eine log (2 ) 


Heat capacity at constant pressure. 
A temperature constant defined as follows: 


18.44 
AR. 





Cy = - log (HF). 


Heat capacity at constant volume. 
Maximum allowable specific force level in Earth g's 


Threshold value of accelerations detectable by specific 
force measuring subsystem in Earth g's. 


Drag force (1lb.) 


Kinetic Energy (ft.-lb.) : CKin) z Ow 


Potential Energy (ft.-lb.) : C(pot) —— % et 


Total Energy (ft.-lb.) : Chie a Cnt eee) 


Dimensionless kinetic energy per unit mass. 


Res = Ekin) 
in 


MVS 
Dimensionless potential energy per unit mass. 


_ Evpot) 
“(pot) ~ “My2 


Dimensionless total energy per unit mass: 


(tot) 
B(tot) = “My2_ 
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rxj| 


Specific force; i.e., force per unit mass. ft./sec.* 


Ba) : Drag specific force vector. ft./sec.* 
F(14) : Lift speeafic force véctor. ft. /sec.* 


(th) ; Thrust specific force vector. ft./sec.* 








£ Dimensionless specific force (in surface g's of the planet). 
f= 
(m)O 
if Dimensionless specific force in surface Earth g's. 
ae all 
a GE 
G : Gravitational field intensity. (ft./sec.“) 
Gino : Mean value of gravitational field intensity at surface of 


planet 0. (ft./sec.“) 


Gin)E : Mean value of gravitational field intensity at surface of 


= 2 
Earth. Gon)E = 32.22848 ft./sec. 
- ; Spherical component of gravitational mass attraction. 
a TaMo 
sp oe 
=: Planet's gravity field. (ft./sec.“) 
g=G- Wro * (Wio % R) 
E(m)E ; Mean value of gravity field at surface of Earth. 


E(m)k = 32-17405 ft. /sec.* 


H ; Altitude of vehicle above planetary surface. (ft.) 
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H 


Sollee . 


Koad 


Dimensionless altitude of vehicle above planetary surface: 


H 


h=y>"— 
R 
(m)O 
n&h modified Bessel function of the first kind with argu- 
ment x. 


Specific impulse of the rocket propellant (sec.) 
Total impulse of rocket thrust (1b.-sec.) 


Symbol used in complex numbers; defined as follows: 


{reas oe 

Exponential decay parameter of planetary atmosphere @eg,-+) 
e = ene 

Stefan-Boltzman constant = 4.81 xX 107! Bru £t7* sec7t ne 


where °R represents degrees Rankine. 


Vehicle surface radiation emissivity. 


Radiation emissivity of atmospheric gas. 


Dimensionless decay parameter of planetary atmosphere 


a 


= K Rimyo 


=feulos 
ee eri 


Ratio of rate of total heat input at any point on the 
surface of the entry vehicle with respect to the rate of 
total heat input at the stagnation point. 


K) = GG 
Lift force (lb.) 
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£ : 
ly : 
la 

& 


a 
M 
m : 
7 


One-half the latus rectum of the ellipse: 1 = a(1 oS 


Dimensionless semi-latus rectum of the ellipse. 


1 
1 = =— 
N Rimyo 


.S 


Unit path length of atmospheric gas. 
Laplace transform of f(T) 
Mach number. 


Instantaneous mass of vehicle (slugs) 


° ® e oe M 
Dimensionless instantaneous mass of vehicle: m a0 
initial 


Mass of planet. 
Mass of Earth. 


Constant defined as follows in section 10.4: 


» Ls ~k Ar 
N = 5D r; e 


This number must be less than about 0.3 for degenerate 
orbital solution of 10.4 to be accurate. 


iL 
Dimensionless lift load factor: "L = MEG 
m 


D 
Dimensionless drag load factor: N= MG... 
(m)O 


Angular momentum per unit mass of the vehicle; twice the 
areal velocity of the instantaneous ellipse. 


Tfab z 
P = RVig =e 4: /sec. 





Dimensionless angular momentum. 


P 


aS R(m)o3/2 G(m)ol/2 
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UI 


Hy | 


Prandtl number. 
Total heat input per unit surface area of entry vehicle 
(BTU/ft.*) 


2 = & + @, 

Total convective heat absorbed per unit area (BTU/ft.“) 
Total heat entering unit area by gaseous radiation (BTU/ft.*) 
Total heat reradiated from unit surface area of vehicle 
(BTU/ft.*). 
Radius of curvature of vehicle in vicinity of stagnation 
point (fte 
Universal Gas Constant = 8.31 X 10° ergs/°Kelvin (i.e., 
work per degree Absolute temperature). 
Radius vector from center of planet to vehicle. 
Dimensionless radius vector from center of planet to vehicle. 
om 

R(m)0 





Incremental change in dimensionless radius from initial 


value, Ar=re-=j0<.973xr 


_ V x length 
ak 72 


Reynolds number, 


Equatorial radius of planet O. 
Mean radius of planet QO. 
Polar radius of planet O, 


Reference area of entry vehicle used in drag and lift 


computations (£t.~) 
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Sw 


(Ts ) max 


T On)O 


<I 


ae 


Wetted area of entry vehicle (£t.*) 
Complex variable used in the Laplace transform. 


Period of elliptical orbit: 


3 
en ae . 
a AT 


1% Mo 


Dimensionless period of elliptical orbit: 


oe Ve /R 


(m)O 

time (sec.) 

Burning time of rocket (sec.) 

Stagnation point temperature (degrees Rankine) 


Maximum stagnation point temperature. 


Mean temperature of planet's atmosphere (degrees Kelvin) 
Stagnation temperature of atmospheric gas (degrees Rankine) 


Auxiliary quantity defined as follows: 


U = PGi) Romo 5_ 
Zz M 


independent parameter having no physical meaning 

introduced to describe propellant mass flow. 

Velocity vector of vehicle with respect to coordinates 
rotating with the planet. 

Dimensionless velocity with respect to coordinates rotating 


with the planet. 


< 
H 
ar|> 


Velocity vector of entry vehicle with respect to inertial 


coordinates. 
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wr : Dimensionless velocity of vehicle with respect to inertial 


coordinates. VieG Vc 


Av : Incremental change in dimensionless velocity from initial 


value. AV =V- Fi 
av. =VvV+-Vv 
if L I. 
Voam) : Velocity vector of vehicle with respect to the planetary 
atmosphere (Note: V = Vow) for no-wind condition). 
le tens) : Velocity vector of planetary atmosphere with respect to 


s 


inertial coordinates. 


Vocam) : Velocity vector of atmosphere with respect to coordinates 


totating with the planet (i.e., wind). 


a : Equivalent exit velocity of rocket propellant (assumed 
constant) 
Ve : Dimensionless equivalent exit velocity of rocket propellant. 
_ we 
Ye ~ Ve 


2 
V : Circular orbital velocity: Va = (GR)? 


Ve: Dimensionless circular orbital velocity: ve. = Vc/Vs 


Vg : Circular orbital velocity at surface of planet: 
3 
oe Gn)0(m)o? 
$V Velocity impulse imparted at trajectory modification point 


by retrorocket system. 
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Ov 


= 


IO 


J] 


6 


Dimensionless velocity impulse imparted at trajectory 
modification point by retrorocket system: 


a 


Dimensionless velocity impulse imparted at trajectory 
modification point by retrorocket system; non-dimensionalized 
with respect to circular satellite velocity at the trajectory 


modification point: 
V 
Sv, = 
Cc Vo 
Daily sidereal angular velocity of the planet; a vector 


quantity directed along the polar axis (North) of the planet. 


Vector angular velocity of the center of gravity of the entry 
vehicle with respect to inertial coordinates. 


Distance flown - measured at surface of planet. 


Dimensionless distance flown: 


X 
Rem)o 





Ay = 


Constant defined in section 10.3 as follows: 


x = FE sinh sinh (S38) | 


A function of angle of attack, angle of sweep, and Mach 
number in convective heat transfer equation (7-1). 


A constant defined in section 10.3 as follows: 


C3 5)| 





yy = Ob cosh| 5 sinh (SS 


Chapman's transformation variable. Equivalent to: 


v cosy sinyY 
E lk 


a 
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A.4.3-2 Multiple Letter Symbols 


(AC) : Dimensional acceleration constant for entry into the 


planetary atmosphere: 


_  Q<Sty Gimo Rim)o slugs_ 
(AC) Z Gime £42 


(AC) a0 : Ratio of acceleration constant of the planet to that of 


the Earth: 


(AC) 
(AC) zo = (AC), 


3 





(HR) 29 : "Heating Ratio" of the atmosphere of planet with respect 
to Earth: 
a z Se, , iii to>1) 4 
(HR) = ZECSL) Gim)o a / Rm)o as Fe 5 ( ‘ 
= (suE]\ Game] | Rime Ke PR. ( 5+) 
%e 
(HF) 5: "Heating Function" of the planetary atmosphere. 


L 


es 
(HF), = oan) (Remo \* (HR)... (slug /#t”) . 


(HF) 59 : Ratio of the heating function of the planetary atmosphere 


to the heating function of the Earth's atmosphere. 





at O 
(HF) 59 = (HF), 
(TC) : "Temperature Constraint Constant" -- a constant for a 


particular vehicle entering the atmosphere of a particular 


planet: 
mm 


iil. Bp 2a (ar) (ry* 


(TC) = T. 
ma 


~,, 





(VF) : Vehicle function, 


(VF) = —= 


Krad] 


Values of the above quantities which are constant for the planet are 


listed in Table A-~2 for the terrestrial planets. 


Table A-2: Approximate values of Heating 
and Acceleration Constants for the 
Terrestrial Planets 


2.86 X 10° ss 2.49 X 107 Sue 
ft 


07,0136 





A.4.3-3: Greek Symbols 


OX : angle of attack 


6 : Bearing angle of entry vehicle as seen by an inertial 


Waly 


observer in the vehicle; angle between North (ly, ) and 
tangential velocity vector of vehicle (directed along 14) ; 


6, : An auxiliary real variable (Eq. 7-26). 


¥,¥ : Flight path angle in inertial framework and planetary 
axes, respectively. 


% : Ratio of heat capacity of the atmosphere of planet O. 
Cp 
De 


Ye i Ratio of heat capacity of the Earth's atmosphere. 


eugetie vy at P gm 
uf : Universal Gravitational Constant = 6.658 X 10 ae eme 


a) 


Propellant mass flow rate (slugs per sec.) 


[° : Dimensionless propellant mass flow rate: 





a = a imo 
N Me y¥ Gimjo 


$5( ) 3: Variation of the quantity in parenthesis. 


6( ) : Incremental change of quantity in parenthesis from initial 
value; e.g., Ar= r-T. 


€ : Eccentricity of ellipse E=/l- (b/a)* 


ey : Component of eccentricity along major axis: Tt = €cos 6 
E> : Component of eccentricity along minor axis: ¢, = €sin 6 


Cp ; Eccentricity of the Earth: Cp =  |0.0067226700 
5 : An auxiliary real variable (Eq. (7-26)) 


E> ; An auxiliary real variable (Eq. (7-32)) 


wD 


angle measured in the plane of an elliptical trajectory 
from the ascending line of nodes to the line of apsides. 


An auxiliary variable (Eq. (7-27)) 


True anomaly; angle measured from line of apsides in the 
ellipse to the position of the vehicle. 


Angle defined as follows in equation (9-97): 
E 4 
2 Ca; 
8 = (% ate “Sg DOL | 


Geocentric longitude of the vehicle; i.e., equatorial angle 
measured from planet's Xo axis to the vehicle. 


Inertial geocentric longitude of the vehicle; i.e., equatorial 
angle measured from planet's Xs axis to the vehicle. 





Inertial geocentric longitude of ascending line of nodes; 
i.e., angle measured in equatorial plane between Xy and Xp. 


Lagrange Multipliers 

Geocentric latitude; angle measured between the radius 

vector from the planet center to the vehicle and the 

projection of this vector on the planet's equatorial plane. 
Geographic latitude; angle measured between the normal 

to the reference ellipsoid of the planet and the projec- 

tion of this normal on the equatorial plane. 

Coefficient of viscosity of the planet's atmosphere (slug/ftssec) 


Coefficient of viscosity of the Earth's atmosphere. 


Dimensionless constant describing non-spherical component 
of planet's gravitational potential. 


The negative of the ratio of percentage change in radial 
distance of the entry vehicle from the planet center to the 


percentage change in horizontal component of the vehicle's 
velocity with respect to inertial coordinates. 


R/R 


€ is used in the definition of the "Conservation 





Parameter". Conservation Parameter = at 
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Te * 


Pisuy? 


PS ? 


Pi (3.1416) 

Free stream atmospheric density (slug/ft.) 

Sea Level or surface value of atmospheric density (slug/ft.) 
The intercept of the straight line which best fits a curve 
of log vs altitude (not the same as true sea level 


density). 


Atmospheric density ratio: (we es = ee 
(SL) 


Summation of quantities in parenthesis 
V (a 0 
Dimensionless time Ji. == = ~ t 
(m)O (m)O 


angle measured in plane of trajectory in direction of 
motion from line of nodes to the vehicle. 





Gravitational potential (potential energy per unit mass). 
In this thesis, the gravitational potential of the Earth 
is assumed to be: 


, 2 
g = GaMe (1 - “Rene (1-3¢0s 2A)| 
Angle of inclination of the trajectory; angle measured 


from equatorial plane to the trajectory plane (Fig. A-2). 


A constant for the skipping trajectory solution: 


2 2uCstslo: ‘s 


@ =| 7 


Dimensionless angular velocity of the planet about its 
polar axis: 


1 = vo (em 


Average angular velocity in the degenerate elliptical trajectory: 


2e ye “y 
WwW = 7 
Angular velocity in the decaying circular orbital trajectory: 
2 aye *( ly ) 
= t + ——— 
Wk - 1 Cr D 
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A.5 Summary of Physical Constants of the Planets and Atmospheres 


Table A.3 summarizes values of various physical characteristics 
of the planets and their atmospheres used in the numerical calculations 
in this thesis. New observations and peeiet: che revise existing data 
frequently; hence complete agreement is seldom found between any two 
tables of this type. Much of the physical data concerning many of the 
planets is fragmentary; large scale revisions may be expected when 


additional data is obtained from artificial satellites and planetary 


probes. 


Table A.3: Summary of Planetary Atmospheric and Physical Data 


(a) Data on the planet Earth: 


378,388 + 18 meters 
963.3386 statute miles Ref. (55) 
0,926,428 ft. 


Req) E 


9356,911.946 meters 
949.9941 statute miles Ref. (55) 
0,855,969 ft. 


R mn) E = 2ONGes,104 ft, Ref. (55) 


(eq)E - “(p)E 
Ellipticity —s.: = 1/297 = .0033670034 Ref. (55) 
eq)E 


& = 0,0067226700 Ref. (55) 
= 7,2921159 X 107? rad/sec. Ref. (55) 


= = (1+ 8 x 107°)(3, 986329 x 20% cm.?/sec.“) Ref. (56) 
Suet 0.14077500 X 5,17 


Gn) E = 32.22848 ft./sec. 
E(m)E eeeee 405 ft. /sec.* 


=e x 107- Ref. (56) 
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Table A.3 (cont.) 


(b) Data on other planets 


ee [ae 


G(m)o/G(m)E 


Atmospheric COs, No, CO5 He, Ho, 
Gases CH), 


yO) /78 0 


Mean molecular 
wet. (gm/mole) 


W 
mn 


Th 0 (°K) 170 


= 


slug sea 
(sr) rt. 


Msn) 12.1 0.88 0604 
(SEE 


Ve 


APPENDIX B 


PHYSICAL CHARACTERISTICS OF MAJOR BODIES OF THE SOLAR SYSTEM 


The solar svstem is the collection of bodies, of which” thé Earth 


eee err 


is one, that are distinguished from all other bodies in the universe 


by having the Sun as the center of motion. The bodies of this system 


may be classified under the following headings: 


(1) 
(2) 


(3) 


(4) 


Sun 


Planets 


(a) The four inner planets -~ Mercury, Venus, Earth, and 
Mars -= are relatively small, dense bodies and are 
known as the "terrestrial" planets. 


(bo) The next four in distance from the Sun -~ Jupiter, 
Saturn, Uranus, and Neptune -- are often called the 
"major" or "giant" planets. These are relatively large 
bodies composed principally of gases with solid ice and 
rock cores at unknown depths below the visible upper 
surfaces of their atmospheres. 


(c) Pluto, recently discovered (1930), has relatively un- 
known physical characteristics. 


Natural Satellites of the Planets 


Thirty-one moons have been discovered to date (1959). The 
system formed by a planet and its satellites is distinguished 
by the name of the planet, such as the Martian system, the 
Jovian system, the Saturnian system. 


Asteroids 


Some 2000 minor planets have been discovered; others are being 
discovered each year. The largest, Ceres, has a diameter of 
480 miles. Most of the rest are less than 50 miles wide. The 


LO 





(5) 


(6) 


(7) 


majority move between the orbits of Mars and Jupiter. Most 
families of asteroids seem strongly influenced by Jupiter, and 
move in orbits determined by the gravitational attraction of 
Jupiter and Sun. 


Comets 


About 5-10 new comets are discovered yearly. Comets are very 
loose collections of orbital material that sweep into the 
inner regions of the solar system from space far beyond the 
orbit of Pluto. Some return periodically; some never do. 

The body of the comet probably consists of rarified gases and 
and dust; the heads are thought to be frozen gases or "ices". 


Meteors 


Meteors exist in untold millions. Meteor fragments that 
reach the ground are,known as meteorites. The average meteor 
weights 0.0005 ounce ‘57 and over 100 million syEake the 
Earth's atmosphere daily. F.C. Leonard (1952) 58) estimated 
the amount of meteoric material falling upon the earth, 
including the ashes of burned out meteors and unburned micro- 
meteorites, to be 6000 metric tons daily. This estimate was 
much larger than former ones, but is confirmed by the work of 
Hans Pettersson and Henri Rotschi in Sweden. A more recent 
estimate a ay ee Observatory favors a value of 2,000 
tons per day 97) | Soviet scientists announced in August 1958 E 
that their a5; data indicated an estimate of 8 to 10 x 10 
tons per day ° 


Micrometeorites and Dust 


The smallest dust particles (micrometeorites) are concentrated 
ae in the ecliptic plane* or along the orbits of 
comets . Evidence that cosmic dust is concentrated in the 
plane of the ecliptic consists of observations of a faint 
tapered band of light centered along the ecliptic, called 


the zodiacal light. 


A brief summary of physical data on principal bodies of the solar 


system is given in Table B-l. New discoveries together with more complete 


observations and calculations revise existing data frequently. 


Table B-2 summarizes physical and astronomical data collected 


concerning the principal natural satellites of planets in the Solar System. 


Plane of the Barth's orbit about the Sin. 


yy 


Table B.1: Summary of Physical Data on Principal Bodies of Solar System 
Note: Superscript numerals refer to notes at end of table. 
Mass M1. Volume Density Density 


(Earth = 1.0) (Earth = 1.0) | (EArth = 1.0) | (Water = 1.0) 


© Sun 332000 1,306,000 
O Mercury 0.0543 0.06 
OQ Venus 0.8136 0.92 
@Earth |(4.09 x 10¢3slugs) 
OMars 0.1069-.1076 
4+ Jupiter 
A Saturn 

Uranus 
W Neptune 

Pigiito 
€ Moon 


5 


e e 


dO 


e 


% 
Li NM 0 


PHORPWWL Pr 
e e 
ON DM NW SOW SO © 


i ~we 


Equatorial Gravitational Accel. 
Radius R at Surface G 
(eq)o (m)o 


(nauts “mae. ) ft./sece“ (Earth 1.0) 


375,600 

1,350 

3,348 

3440 

1,787 

Jupiter B77 00 

Saturn Boe 
Uranus 13,800, ) 
Neptune ee a 

Pil wee 1, 5 

Moon 938 


Surface Escape | Minimum launch Transit 
Velocity ~ Velocity to Time 
reach Earth 


(t/sec) (km/sec ) (ft./sec.) 


sun 
Mercury 
Venus 
Earth 
Mars 
Jupiter 
saturn 
Uranus 
Neptune 
Putco 
Moon 


oO 
e e 


= 


win EUIOMW EW 


Nh MW ON 
WIR ND ONE WW 
s @ 





mae 5.906 + 001 x oe! gms . 


Table B-1 (cont.) 




























Orbital Distance 
Perihelion | Apheljon 
(x TOE) | (x 107 €e.) 


Mean Distance from Sun 
ee ee Millions 


Units of miles 










sun 
Mercury 
Venus 
Earth 
Mars 
Jupiter 
Saturn 
Uranus 
Neptune 
Paatto 
Moon 





Orbital Speed (ft. /sec.) Eccentricity Orbital Period 


Perihelion Aphelion of Orbit 


















Sun 


























































Mercury 176,000 143,000 88 days 

Venus 11% ,,.000 114,500 224 a 17 hrs. 
Earth 98,500 96,800 365 d 6 hr. 9 min. 
Mars 83,000 75,700 687 days 

Jupiter 43,900 41,800 4328 days 

Saturn 32%, 200 30,800 29 years 167 days 
Uranus 22,600 21,800 84.01 yrs. 
Neptune 17,900 17.9700 164.8 yrs. 

Pluto | 248 yrs. 

Moon 3,480 35200 27a7h 43m 


Orbital Inclination 
to Ecliptic Plane 






Sun 

Mercury 

Venus 

Earth 23° B7' 
Mars 25 -aai0 ' 
Jupiter om ap" 
Saturn 26° 45! 
Uranus 98° 
Neptune 29° 
Pluto 

Moon 


Earth's equator 
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sun 
Mercury 
Venus 
Barth 
Mars 
Jupiter 
Saturn 
Uranus 
Neptune 
Pasato 
Moon 


Sun 
Mercury 
Venus 
Earth 
Mars 
Jupiter 
saturn 
Uranus 
Neptune 
Pluto 
Moon 


Distance from Earth 
(millions of miles) 




















Pllipticie, 


Fi eq)o - yo 





Table B-1 (cont.) 








Total Radtatrone’] ux 






Intensity of 
Sunlight at 
mean distance 
Harth = 1.0 









8.95 
2D 
1.34 
0.58 
0.0496 
0,0147 
0.00364 
0.00149 
0.000859 

















Theoretical Radiation 
Temperature 









Maximum Average 
on ok. 





F(eq)0 


0033670034 
,005215 

0649 
1053 
0714 
0254 















(1) 


(2) 
(3) 


(4) 


(5) 


Table B-1 (cont.) 
Superscript numerical notes of Table B-l. 


These values are based on 1950 observations of Kuiper. This 
diameter of Neptune is 3,300 miles less than usually given. 
The diameter and volume of Pluto by these observations is 
considerably less than earlier values. 


Has no solid Surface. 


Location of solid surface (below thousands of miles of dense 
atmospheric gases covering these planets) is not known; 
hence gravitational figures at surface are not reliable for 
the four giant planets. 


The “astronomical unit", a fundamental unit of astronomy, is 
the mean distance of the Sun from the Earth. “Solar parallax" 
is defined as an arc whose tangent is the radius of the Earth 
divided by the mean distance of the Earth from the Sun. 

For over 2000 years, increasingly accurate determinations 
of the solar parallax have been made. From about 1900-1940, 
astronomers were confident that the value of solar parallax 
was 8.80", with the error in its determination not greater than 
0.01", The asteroid Eros, roughly a cylindrical body 5 miles 
in diameter and 20 miles hong, provided a means for recent 
readjustment in the value of solar parallax. The nearest 
approach of Eros to the Earth is approximately 10 million 
miles; its last opposition (1930) was about 16 million miles. 
Approximately 2500 photographs of Eros, taken in 1930-31, were 
analyzed in great detail by Sir Harold Spenser Jones over a 
10-year period; his determination of solar parallax gave the 
unexpectedly low value of 8.790" with an estimated accuracy of 
1 part in 10,000. Considering the scope of Jones' work, his 
figure hardly seemed open to question, 

Recently, however, a new determination of solar parallax, 
along with other fundamental data, was published by Eugene 
Rabe. He made a full rediscussion of all the observations of 
Eros from 1926 to 1945 and secured a solar parallax value of 
8.79835" + 0.00039", which is much closer to the old value of 
8.80" than that of Jones. As a byproduct, the ratio of the 
masses of the Earth to Moon is 81.375 and of Sun to Earth 
plus Moon is 328,452. He also secured improved values for the 
masses of Mercury, Venus, and Mars. There are indeed several 
very small asteroids which approach the Earth closer than does 
Eros, and hence might be expected to give better results. 

Thus far, however, they have proved unavailable because of 
their extreme faintness and the uncertainty of their orbits. 


Period of rotation on accuse 


MALS 


(6) 


(7) 


(8) 


(9) 


(10) 


(er) 


(12) 


Table B-1 (cont.) 





Mercury's rotation period equals its period of revolution about 


the sun; hence one side always faces the sun. Temperatures 
here are as high as 750° F. above zero; on the other side, 
perhaps hundreds of degrees below zero. 


Mars' moons, Deimos and ‘Phobos, are about 10 miles in dia- 
meter. The inner one (Phobos) revolves so fast that it rises 
and sets three times in a Martian day. 


Four of the moons have diameters of 2,300 to 3,200 miles and 
revolve about Jupiter in 2-17 days. The other 8 moons are 
less than 100 miles in diameter. One, very close to Jupiter, 
revolves at over 1,000 miles per minute. The 12th moon was 
discovered by S.B. Nicholson of Mt. Wilson and Palomar Obser- 
vatories in 1951. It is only 14 miles in circumference, 

2 x 10~ miles from the planet, and circles Jupiter in 700 
days retrograde motion (east-to-west). Four of Jupiter's 
moons are characterized by retrograde motion. 


One of Saturn's moons is larger than the Earth's moon and is 
believed to have an atmosphere. All nine moons are outside 
the rings characteristic of this planet. 


The fifth satellite of Uranus, Miranda, was discovered in 
February 1948 by Kuiper. 


The second satellite of Neptune (Nereid) was discovered on 

1 May 1949 by Kuiper. It has a period of 359 days and ranges 
from 830,000 to 6,100,000 miles with major orbit axis of 
0.074 a.u., eccentricity of 0.76 (greater than for any other 
satellite), and inclination of 5° to ecliptic plane. The 
first satellite of Neptune is called Triton. 


The closest approach of any planet to Earth. 


WG 





Table B-2: Summary of Physical Data on Principle 
Natural Satellites of Planets 


Parent Distance from Period of Orbital 
Planet Planet (KM) Revolution Eccentricity 





ot Pa. 
Besos teres 9.35 x 102 ; ee 39.20 0.021 
Deimos Mars 2.35 x 10° 1° 6" 17,9" 0.003 
» 29 


JI Jupiter N22 ° X 102 

; Geri: tee ie a3. 

ie7t 108 
1766+*x 16 

1.815 x 10? 
10 

Lo, 

104, 
10 








Cy Cy Cy Cy Cy Cy Oy Cy 
a 


FMinW Fn FrHMinW 
Otn OFWEH OW nN 
2° 


Mimas 
Enceladus 
Tethys 
Dione 
Rhea 
Titan 
Hyperion 
Japetus 
Phoebe 


NO 
do 
> 


OO 


or 


st 


Nr FN 
NON 
py a 

© ® 


te 
MES SS Ss a 


ar 


® 


~~) 
jo 


® 


® 


Ariel* 1939 same 
Umbriel* 2.673 x 102 
i tania* 4.387 x 10% 
Oberon* 5.866 x 10 

Miranda 


10? 
gore 


Triton* 3 Be 
Nereid 1. 3eto 
x 10 


* Satellites indicated by asterisk are characterized by retrograde 
motion. Motion of Miranda and Nereid is not known. 


**k Most of this data was taken from reference (62). 
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Table B-2 (cont.) 


Theoretical 


Velocity Radiation 
Temp. °Kelvin 
Max. Ave. 













4.03 


Phobos es 
Deimos 8 



















Cy Cy Cy Cy Cy Cy Cy Cy Gy Gy Gy 





Mimas 
Enceladus 
Tethys 
Dione 
Rhea 
Titan 
Hyperion 
Japetus 
Phoebe 





Ariel 
Umbriel 
Tania a 
Oberon 
Miranda 


Triton | 4500 NOE? 2: pee: 9177 
Nereid 


2 
*** Mass of Earth = 5.966 + .001 x 10“? ems = 4.09 x 10 7 slugs 
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APPENDIX C 


GRAVITATIONAL MASS ATTRACTION AND THE 


ACCELERATION OF GRAVITY 


Geil: Gravitation 

According to Newton's law of gravitation, every particle in the 
universe attracts every other particle with a force varying directly 
as the product of the two masses and inversely as the square of the 


distance between them. 


MM 
Force = em C=!) 


% is a universal constant whose numerical value is approximately: 


x = 6.658 x 10714 ae 2 (c-2) 


KG = sec. 


A single fixed mass point of mass M will exert a force on any other 
particle according to equation (C-1). All space is affected by the 
mass M; a field of force is set up around the mass. The same general 
statement is true for any system of bodies; at each point in space a 
force vector can be drawn which is the vector sum of the forces acting 
on the mass M by the system of bodies M,, Mp, ..-M,. The tangents to 
the force vectors at each point in space give the direction of the 


force exerted on the mass at every point. The curves map out the field 


4449 


and are called lines of force. 

For the case of a single particle, the lines of force consist of 
straight lines radiating in all directions from the particle as a center, 
The lines of force give information as to the direction but not as to 
the magnitude of the gravitational forces. The latter is specified by 
introducing the intensity of the force field: the force per unit mass 
exerted on a particle at any given point. The intensity of the gravi- 


tational field is denoted by G; for a single particle of mass M: 





G= Force =-Fs= = (C=3) 


ma R 


G is a vector whose magnitude is given by (C-3) and whose direction 

is inward along the line connecting the masses. For a system of particles 
the intensity of the field is calculated by vector addition, the resul- 
tant intensity being the vector sum of the intensities due to the indi- 
vidual particles. 


(19) 


Gravitational forces are conservative » and it is normal to call 
the fields conservative*. It is convenient to introduce the gravita- 
tional potential @ , defined as the potential energy per unit mass. 

For the case of the field of a single particle, the gravitational 


potential is 


_ eM 
¢d — (C=) 


and the gravitational intensity vector G is the negative gradient of 


* A conservative field is defined as one in which the total work done 
is zero in moving a particle around any closed contour in the field. 
In moving a particle between any two points in a conservative field, 
the work done is independent of the path chosen between the points. 
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the potential ¢g : 


a (C5) 


The gravitational field intensity G at a point outside a spherical 
planet is determined by considering the planet to be composed of thin 
spherical shells, each of constant density. Each spherical shell acts 
as if its mass were concentrated at the center. Therefore, at any point 
outside the planet, the gravitational field intensity vector is directed 
toward the center of the planet with a magnitude directly proportional 
to the product of the gravitational constant times the mass of the planet 
and inversely proportional to the square of the distance from the center 


of the planet to the point. 


~ Igo 
G = = 


(C-6) 





I aI 


where My is total mass of the planet. This result is identical with 


that for a single particle given in equation (C-3). 


C.2 Gravitational Field of Oblate Spheroidal Planets 


If the planet were a nonrotating fluid mass, it would form a 
sphere under the influence of its own gravitational field. Because of 
rotation about an axis through the planet, the associated centrifugal 
force field causes the planet to bulge at the equator. In the case of 
Mercury, the axial rate of rotation is extremely small”, hence the planet 
must be very nearly a perfect sphere (except for surface anomalies -- 


* Mercury's day and year are both equal to approximately 88 Earth days. 
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mountains and valleys). Little is known about the rotational rate of 
Venus since the surface is not visible through the opaque atmosphere. 
The best guess is that Venus, too, has equal days and years (224.7 
Earth days). Consequently, Venus should be almost a perfect sphere 
(except for mountains, etc.) 

Earth, Mars, and the giant planets have relatively high rotational 
velocities . Mars rotates at almost Earth's rate, hence should have 
nearly the same equatorial bulge as Earth. Jupiter, Saturn, and Uranus 
have days that range from 9.8 to 10.7 hours, and Neptunes day is slightly 
under 16 hours. Consequently, the equatorial bulge of these giants is 
much more pronounced than Earth -- a fact verified by astronomic 


photographs. 


North 





Axis of rotation 


Angular velocity vector 
due to rotation about 
polar axis 


Fig. C.l: Gravitational Field of Oblate Spheroidal Planet 


* See Table B.1l. 
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Due to the interaction of gravitational mass attraction and the 
centrifugal force due to axial rotation, the actual shape of rotating 
planets is nearly ellipsoidal, or oblate spheroidal. Because of the 
non-spherical mass distribution of the planet, its gravitational field 
is not spherically symmetric. A first order approximation of the non- 
spherical component of gravitational field intensity may be made by 
considering it as a quadrupole moment. Fig. C.1 shows two sources and 
two sinks of equal strength located on the polar axis of the planet 
near its center to represent the non-spherical mass component. The two 
sources are located at the center and the two sinks are symmetrically 
spaced on the polar axis above and below the center. The positions of 
the sinks approach the center, while the product of strength and spacing 
remains constant. The quadrupole thus is similar to a pair of dipoles 


that have mirror symmetry. The gravitational potential including this 


oblateness quadrupole may be considered to be (63). 
Pp 
6 =- ty 1 eee ati ~ 3c0s 2h)| (C-7) 
G=-v7¢@ (C-8) 


In equation (C-7), 7 is a dimensionless constant for the planet, 
‘is geocentric latitude, and Ri eq)o is equatorial radius of the 


planet. For the planet Earth, reference (63) lists the value of Vv as 


follows: 


6 vw.= (1,637 = ces) x 107° 


E (C-9) 
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Ref. (56) gives a value of VY for the international ellipsoid’ of 


Barth as follows: 


6 V, = 1.638 X 107? . (C-10) 


and the value of Ye Mo for Earth as follows: 


(Geet 8a x 107?) (3.986329 X aioe om.?/sec.*) 


Be 
tl 


0.14077500 x lot” £t.3/sec.* (C-11) 


C.3 Gravity 


The surfaces for which G is constant are not oblate spheroids, 


nor does one surface coincide with the assumed geometric figure of the 


planet. The surface g equal to a constant are everywhere normal to 


the gravitational field intensity vector of the planet. The vector 


sum of the planet's gravitational specific force**and centrifugal force 


* 


The Hayford Spheroid of 1909 is the international reference 
ellipsoid of the Earth. At the Internationale Geodetic and Geophy- 
sical Union of 1924, the dimensions of the Earth summarized in 
Table A.3 were adopted (see Table 716, page 570 of ref. 55). 


Gravitational specific force G is identical with gravitational field 
intensity. Centrifugal specific force, a reaction force, is equal 

in magnitude but opposite in direction to the associated centripetal 
acceleration. Specific force is expressed in force-per-mass 
dimensions. 


+ 


per unit mass is defined as the planet's gravity field intensity.” 


g =- VO - Wy X (Wy X R) (C-12) 


where 


@ is given by equation (C-7) 


R is the vector from the center of the planet to a point 
on or external to the surface of the planet 
Wro is the angular velocity of the planet about its polar 


axis, 


The tangential component of the first term in Equation (C-12) 
prevents the oceans of the Earth from flowing toward the equator. 
The direction of gravity is normal to the surface of the planet; if 
gravity anomalies are neglected, this direction is unique at each 
point on the surface, The time variation of the direction of gravity 


(64) 


on the Earth is less than 0.05 microradians and is mainly due to 
tidal effects. Fig. C.2 shows the relation of Earth's gravitational 
specific force and gravity specific force. The angles are exaggerated 


in this figure for clarity. 


* The specific force of gravity (the force of gravity acting on a 
unit mass) is identical with the gravity field intensity. 


oo 


G tangential component 


W5sin A “ — Wry X (Weg X R) 
4’ssin A Z] 
R </ ye 
vA Balances G 
O Tangential 
a g component 


Note: 

Angles are exaggerated 

for clarity in presentation. Numerical values 
are approximate values for Earth. 


Fig. ©.2: Relation of Gravitation and Gravity. 


(19) The angular velocity of the Earth about its polar 
axis is 


Wro = 7-2921159 X 107 rad./sec. (C-13) 
The international standard gravitational conversion factor is the 


apparent acceleration of gravity at latitude 45° 32' 40" : 


E(m)E 7 32.17405 ft./sec.- ~ (C-14) 
Since 


E(m)E = °(m)E - We Rin)g cos A (C-15) 


the following is determined for Gm)E° 


G(nyp = 32-17405 + 0.05UN3 = 32.22648 ft./sec.° (C-16) 
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Appendix D 


FIGURE OF THE PLANET AND DEFINITION OF 


NAVIGATIONAL PARAMETERS 


D.1l Figure of the Planet 

The figure that would be taken by a fluid body with the mass 
distribution and axial rotation of the planet is defined as the figure 
of the planet. The figure of the planet is an equipotential surface 
of the planet's gravity field. 

In connection with the actual configuration of the Earth, study of 
the motion of Vanguard I has shown the Earth to be slightly pear-shaped, 
with the more narrow half above the equator and the sag below. There is 
an approximate 50 ft. rise in mean sea level at the South peas Fn, 

This deformation is interpreted to mean that the crystal layers of the 
Earth have more strength and less elasticity than previously assumed. 

The figure of other planets in the solar system is not nearly as 
well known as that of the Earth. Unmanned planetary probes, which are 
planned for launch in the 1960's should provide the first quantitative 
information obtained from data sources near the planet to augment astro- 
nomical observations and calculations accumulated over many years. The 
oblateness of the planets is generally estimated from astronomical 
photographs and from calculations based on observed rotational rates. 


In general, the terrestrial planets are believed to be more spherical 
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than the giant planets because of the higher rotational rates and smaller 
densities of the latter. 

The figure of the planet may be approximated by a reference ellip- 
soid, an analytical figure closely approximating the planet. The 
Clark Spheroid (1866) is the reference ellipsoid of the Earth used for 
North American triangulations. The Hayford Spheroid (1909) is the 
international reference ellipsoid of the Earth adopted by the Inter- 
nationale Geodetic and Geophysical Union of 1924. This is called the 
geoid, a "square shouldered" ellipsoid with a slightly undulated 


sarexee) | 


The surface of the geoid is represented by mean sea level. 
Variations in the elevation of the geoid relative to the closest 
reference ellipsoid are approximately one percent of the topographic 


variations in sie tion 


D.2 Navigation Parameters 

A manned vehicle entering the atmosphere of a planet must have 
some means of identifying its present position and the location of its 
destination. Three coordinates are required in order to uniquely specify 
position in three-dimensional space, Of the infinite number of possible 
sets of coordinates that may be used to identify position near a planet, 
three are most common: altitude, latitude, and longitude. 

Three navigational reference frames are discussed in this section 
and a comparison made between navigational parameters expressed in each 
of these frames: 

(1) <Astronomic Frame 
Recognizes the actual planet configuration. The 


coordinate frame is centered at the guidance point and 
has one axis along the true vertical (gravity vector). 
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(2) 


(3) 


The other two axes depend on the guidance grid desired 
(e.g., astronomic latitude and longitude). 


Geographic Frame 


Recognizes the ellipticity of the planet. The 
coordinate frame is centered at the guidance point and 
has one axis along the geographic vertical (normal to 
the surface of the reference ellipsoid). The other two 
axes depend on the choice of guidance grid (e.g., 
geographic latitude and longitude). 


Geocentric Frame 


Utilizes spherical polar coordinates. The coordinate 
frame is centered at the guidance point and has one axis 
along the geocentric radius. The other two axes depend on 
the choice of guidance grid (e.g., 1;, la , In triad or 
ip» lg , ly triad discussed in Appendix A). 


It was pointed out in Appendix C that, neglecting anomalies, the 


direction of the gravity vector Seeds unique at each point on the 


surface of the planet. The direction of the gravity vector is defined 


as the true vertical. The uniqueness of the true vertical is the basis 


for astronomical position. Two angles are sufficient to identify the 


position of any point on the surface of the planet: 


(1) 


(2) 


Astronomical latitude: 


The complement of the angle between a line parallel to 
the polar axis of the planet and the local gravity vector. 


Astronomical longitude: 


The angle about the polar axis of the planet between 
an arbitrarily chosen reference vertical and the local 
vertical. Due to the arbitrary choice of reference 
vertical, longitude is en inertially unique, although 
change of longitude is. ) 


The geographic vertical” is defined as the normal to the surface 


* 


In the strictest sense, the term "geographic" refers to the Earth 
(geoid) only. In this thesis, "geographic" is used in a more 
general sense to refer to quantities associated with the reference 
ellipsoid of any planet. 


si) 


of the reference ellipsoid. Because the figure of the planet does 
not have a smooth surface, the true vertical is not, in general, parallel 
to the normal to the reference ellipsoid at the same position. The angle 
measured from the normal of the reference ellipsoid to the true vertical 
is called the deflection of the vertical, or station error. For the Earth, 
station error is generally less than 0.3 milliradians; over continental 
land masses, it rarely exceeds 0.1 911) ae 

Identification of position on the surface of the reference ellip- 
soid is accomplished with the following two angles: 

(1) Geographic Latitude ( Ag) : 


Angle between the normal to the reference ellipsoid and 
its projection on the equatorial plane. 


(2) Geographic Longitude: 
Angle between the projection of the normal to the refer- 
ence ellipsoid on the equatorial plane of the planet and 
an arbitrarily chosen reference meridian. 

The geocentric vertical is defined as a unit vector directed from 
the point at which the guidance is taking place toward the center of the 
planet. The angle measured from the geocentric vertical to the normal 
of the reference ellipsoid is called the deviation of the normal. The 
maximum deviation of the normal for Earth is approximately 11.5 minutes 
of arc. 

identification of position by means of spherical polar coordinates 
may be accomplished fe ouen specifying radius from planet center (or 
altitude above the surface of the planet) and by two independent angles. 

One set of independent angles that may be used to specify position 
are great circle parameters, ie., angular position along and across the 


great circle track. Even though the term great circle has no meaning 
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on any surface other than that of a sphere, it is a useful concept in 
navigation over non-spherical bodies. An astronomic great circle course 
is defined by the property that the normals to the planet's surface 
along such a course are parallel to a single, fixed plane. A geographic 
great circle course is defined as a course for which the normals to the 
reference ellipsoid of the planet are parallel to a single, fixed plane. 
The astronomic and geographic great circle courses are not geodesics, 
therefore a point on the planet's surface following such a course would 
experience a geodesic acceleration.” 

Another set of spherical angles used to specify position are: 

(1) Geocentric Latitude (/\): 


Angle measured between the geocentric radius and its 
projection on the equatorial plane of the planet. 


(2) Geocentric Longitude (A): 
Angle in the equatorial plane of the planet between the 


projection of the geocentric radius on the equatorial 
plane and the longitudinal reference meridian. 


When the equatorial plane is taken as the true equatorial plane, 
geocentric and geographic longitudes are identical and astronomic 
longitude differs by a component of station error. A comparison of 


geocentric and geographic latitude is made in the next section. 


* A horizontal component of acceleration arising from the curvature 
of the course in the horizontal plane. 
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D.3 Comparison of Geocentric Latitude and Geographic Latitude 


Most of the analysis of this thesis is based on geocentric 
angles as prime dependent variables in the statement of the guidance 
problem of vehicles entering planetary atmospheres. Since position 
on maps and charts is conventionally specified in terms of geographic 


latitude, it is instructive to compare geocentric and geographic 


latitudes. 
/\ = Geocentric Latitude Normal to Reference Ellipsoid 
Ag= Geographic Latitude 
i 
Reale : R (5) = Equatorial and Polar o 


radii, respectively. 





Geographic 
Horizon 
Xo 
Fig. D.1l: The Reference Ellipsoid 
It is shown in reference (66) that: 
oO ae Al\ 
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dR (D.2) 


Therefore: tan (Ag - A) ee +5 dA. 


The equation of the ellipse in polar coordinates is: 


Z Z 








Ry ez 
ob = 
Rae ( R(pyo ! (Dp 3) 
where: 
R, =R cos /\ 
R, = Resin /\ 
Therefore: 
( cos /A | sin AY 
: \ Rreqyo “ (an) =a! al 
selying for R: 
R (ego 
R | cos AA + (ise) sin A\z (Ds) 


The eccentricity of the planet, E , is defined” by (see Fig. A.6): 


Ripro_ 
<_. e 


Equation D.5 thus becomes: 
Rieq)o 


R = 2aeinzA 
‘gi a — Ae (D.6) 


* The ellipticity of the Hayford ellipsoid is given as: 
Ellipticity = Frege Ripye = 1/297 


Riege 
Thus: C2 = 7593 x (Ellipticity) 
ee = 0,0067226700 
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Using Eq. (D.6) to eliminate R from Eq. (D.2) gives: 


2 
7 e~ sin AcosA _ _e sinAcosA 
tan(A_ -A) (Fe) + ESineA ) = (1 = €2c0s@A ) CDi) 
Equation (D.7) may be written: . 
tan( A, -/\) = €~ sinNcosA + Sagas A at (D.8) 


This expression is an exact* relation between geocentric latitude /\ and 
geographic latitude lees Examination of this equation shows that 

Ch -/\ ) © 12 minutes of arc over the Earth's surface. It is probably 
less than 12' for Mercury, Venus, and Mars; but (A, =e Nes) is probably 
somewhat greater than 12' for the giant planets (high rotational speeds -- 
small density). Using the small angle approximation of Ap =/\ )ofeg 


terrestrial planets give : 


Mercury 


iN me /\ 4+ e~ sinAcos/\ — only (D.9) 


Mars 


D.4 Radius of Spheroidal Planets as Function of Latitude 





Equation (D.6) gives an expression for the radius vector to a 
point on the surface of the ellipsoidal planet in terms of geocentric 
latitude, eccentricity, and equatorial radius. This may be expanded 


by the binomial theorem to give: 


2 L (esinAy’ acl (esin 
= R —sfny +e i 
"surface eqpl é (1-€?) i= 6 i] fe ca a 


« "Exact" to the extent that the reference ellipsoid approximates the 
figure of the planet. 
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For the terrestrial planets where RU)o and R do not differ 


(eq)O 
greatly, Equation (D.10) may be reduced to a first order in €*: 


Mercury 
~ ee Venus | 
Hee Pe ate - “(BGO (1 - eee Earth ony (D.11) 
Mars 


D.5 Position Reference for First-Time Entry into the Atmospheres of 
trange Planets 


First-time entry into the atmospheres of strange planets presents 
special problems in specifying position when compared to navigation over 
a well-mapped planet such as the Earth. The choice of a suitable lan- 
ding site must necessarily be based on reconnaissance of the planetary 
surface while in orbit around the planet. The orbital altitude for the 
reconnaissance phase must be high enough such that a prolonged orbit may 
persist, yet low enough that fairly accurate mapping of the terrain is 
feasible. 

The navigation of a vehicle flying from one point to another point 
on the surface of a planet is generally based on navigational parameters 
measured with respect to the planet, i.e., latitude, longitude, and 
altitude. It is less common to use parameters identified with the 
particular mission, such as angular displacements measured with respect 
to the great circle course along-track and across-track. 

The entry mission to the surface of a strange planet, on the other 
hand, does not originate from a point on the planet's surface. The entry 
mission originates from the reconnaissance orbit, which, if the perigeal 
altitude is sufficiently great, is very slowly changing with time. The 


* "Strange planets", as used in this thesis, signifies any planet on 
which a vehicle with human occupants has not landed. 
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non-spherical component of the planet's gravitational field causes the 
line of nodes to rotate slowly with respect to the inertial framework. 
Drag forces result in energy transfer from the vehicle to the planetary 
atmosphere, but for sufficiently high orbits, this transfer causes 
negligible change in the satellite orbit over periods of time comparable 
to that required for entry once retro-rocket thrust is generated. 

A basic position reference available during the course of entry is 
the original reconnaissance orbit. If the entry vehicle is launched 
from a mother satellite, then the mother satellite, which remains in the 
reconnaissance orbit, may track both the entry vehicle and the pre- 
selected landing site and transmit this tracking information to the 
navigational computer of the entry vehicle. In this way, the parent 
satellite replaces ground tracking stations which are used as an external 
source of tracking information for Earth satellites and entry vehicles. 

In the event that no parent satellite exists, then a navigational 
satellite to serve the same purpose may be @®posited in the reconnaissance 
orbit by the entry vehicle prior to initiating the entry phase. 

Since a navigational scheme such as outlined briefly above uses the 
recorinaissance orbit as the basic reference from which to measure 
positions, it may be found convenient to express position and to carry 
out the guidance computations in terms of elliptical parameters such as 
those presented in Chapter 4. The landing site may be considered to be 
a target moving in three dimensional space with respect to the near- 
stable reconnaissance trajectory represented by the mother or naviga- 
tional satellite. The entry vehicle is also moving with respect to the 
reconnaissance trajectory. The entry problem is therefore similar to the 


Fire Control problem with the entry vehicle (projectile) fired from the 
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parent satellite (gun)to hit the moving landing site (target). The 
problem is much more severe than the conventional fire control problem, 
however, because the projectile must be constrained to paths for which it 
will not burn up or encounter accelerations beyond tolerable levels. 

The nominal or programmed path of the entry vehicle may be computed 
in advance as one which the vehicle would fly under standard atmospheric 
conditions starting from the particular initial point and ending at the 
landing site selected in advance. This trajectory must be consistent 
with tolerable accelerations and heating rates (Chapter 7). 

It was shown in Chapter 4 that there are six elliptical elements 
required to specify the position and path of the vehicle. One set of 
six such elements are: 

(1) W, =. (to specify the instantaneous orientation 
of the plane of the trajectory). 

(2) SEs Co c, (to specify the ellipse which matches 
instantaneously the dynamical state of the vehicle) 

(3) (to specify the position of the vehicle in this ellipse). 

P, W, and Aj;7 are constant or very slowly varying with time in the 
reconnaissance orbit. fc and Co are sinusoidal with a very slowly 
changing magnitude of oscillation. d increases monotonically. The 
instantaneous state of the entry vehicle, with respect to the navigational 
satellite in the reconnaissance orbit, may conveniently be specified in 
terms of the six elliptical quantities. The orbit of the navigation 
satellite should be predictable to a fairly high degree of accuracy. 

The predicted values of six elliptical elements for this orbit are part 


of the data stored for use during the entry mission. 
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APPENDIX E : 


THE ATMOSPHERE OF THE PLANETS AND THEIR NATURAL SATELLITES 


E.1 Composition of Atmosphere of the Planets and Their Natural Satellites 


The presence and stability of planetary atmospheres can be predicted 


from the kinetic theory of cies”, 


For a stable atmosphere to exist, 
the root-mean-square molecular velocity of the atmosphere should be less 
than 20% of escape velocity. The RMS molecular velocity is a function 
of temperature (which depends on distance from the Sun), atmospheric 
composition, planetary rotational rates, etc. The relative lilslihood of 
major bodies of the solar system possessing an atmosphere is given in the 


(67). 


following list The higher a planet or moon is on the list, the 


higher is the probability that there exists an atmosphere around it. 


(1) Jupiter* 
(2)  Saturn* 
(3)  Neptune* 
(4)  Uranus* 
(5) Earth* 
(6)  Venus* 
(7) Pluto 
(8) Triton 
(9) Mars* 


(10) Titan* 
1) eaJowiiene TLL 
(12) Jovian I 


* Proof of the existence of an atmosphere by spectroscopic or other 
means has been established for all bodies marked with an asterisk(*). 
Evidence for all other bodies is inconclusive. 
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(13) JovienlV¥ 
(14) Jovian II 
(15) Mercury 
(16) Moon 


A brief description of the atmosphere and climate of each of the 
planets and their natural satellites, adapted primarily from references 


(7) and (58), is summarized below: 


(1) Mercury 


Mercury is a small rocky sphere, about half again as 
large as Moon, that always has the same side exposed to 
the Sun. Maximum surface temperature on the "hot" side 
is estimated at 650-750° F. while that on the cold side 
approaches -400 to -415° F. Mercury is not known to 
have an atmosphere of any significance, nor would a 
permanent gaseous envelope be expected to occur under 
the condition existing on the planet. Its rocky surface 
is probably somewhat similar to that of Moon. 


(2) Venus 


A dense, dusty, turbulent atmosphere, containing 
much carbon dioxide and some nitrogen (but negligible 
free oxygen and water) conceals the planet's surface. 
The atmosphere of Venus contains white particles in 
suspension and is opaque to light of all wave lengths. On 
the basis of all available evidence, it may be presumed 
that the surface of Venus is probably hot, dry, dusty, 
Windy, and dark beneath a continuous dust storm; that the 
atmospheric pressure is probably several times the 
normal barometric pressure at the surface of the Earth; and 
that carbon dioxide is probably the mjor atmospheric 885+ (69) 
with nitrogen and argon also present as minor constituents : 


(3) Mars 


Many questions about surface conditions on Mars are 
still unanswered. It has an appreciable atmosphere and 
its surface markings exhibit seasonal changes in colora- 
tion. Its white polar caps* are apparently thin layers 
of frost of the order of inches or fractions of inches 
thick, The atmosphere is composed primarily of nitrogen, 
carbon dioxide, traces of water vapor, and ey, believed 
to be almost totally absent of free oxygen‘ ” 


* The polar caps of Mars appear during its winter and disappear in summer. 
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Topographically, its surface is flat, with no abrupt 
changes in elevation and no prominent mountains. The 
climate is imagined to be like that which would be encoun- 
tered on the Earth in a desert at 11 miles altitude. 
Noon summer temperatures in the tropics may reach 80° - 
90° F, while predawn temperatures may drop to -100° F, 

h 


There is evidence that some indigenous life forms 
may exist on Mars. The seasonal color changes, from 
green in summer to brown in autum, suggest vegetation. 
Recent spectroscopic studies give evidence of organic 
molecules being responsible for the dark areas 7) The 
objections raised concerning differences between the color 
and infra-red reflectivities of terrestrial organic matter 
and those,of.the dark areas on Mars have been explained 
by Tikhov‘ 0) | He has shown that arctic plants differ 
in infrared reflection from temperate and tropical plants, 
and an extrapolation to Martian conditions leads to the 
conclusion that the dark areas are vegetable lift. Human 
life could not survive without vast environmental modifi- 
cations, but a self-sustaining local animal colony is 
possible. 


(4) The Giant Planets 


The four giant planets are massive bodies of low 
density and large diameter. They all rotate rapidly, 
hence have considerable "flattening" at the poles 
("bulges" at the equator). It is shown in Table B-1 
that their densities vary from about 0.72 to 1.6 times 
the density of water. On the basis of this and spectral 
data, they are considered to be composed of a dense rocky 
core surrounded by a thick shell of ice and covered by 
thousands of miles of compressed hydrogen and helium.* 
Methane and ammonia are known to be present as minor 
wo nnemars io4) of the atmosphere, but no water has been 
observed 471) , 

Jupiter and Saturn show light and dark belts in the 
atmosphere parallel to the equator. These banded clouds 
are slowly changing in Jupiter; they are not as clear nor 
are they changing as rapidly on Saturn. Jupiter's great 
red spot, 20,000 miles long, seems more permanent than 
its cloud belts, but it is believed to be fading. Bright 
spots occasionally appear in the cloud banks of Saturn. 

Saturn's rings are its most prominent physical 
feature. They are probably composed of millions of tiny 
solid particles. They may be material which never formed 
into a satellite, or fragments of a close satellite torn 
asunder by the tidal pull of Saturn, or ice particles. 


ao ease ese se se se 2 =F FF FSF SF SOS ll 


* A "rock-in-a-snowball" structure -- reference (7). 
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(5) Pluto 
Little is known about Pluto beyond the physical 
characteristics listed in Table B-1 and the fact that it 
is extremely cold. 
(6) Moon 
The Moon has no appreciable atmosphere, and its 
surface is probably dry, dust-covered rock that is not 
homogeneous either in chemical composition or topography. 
The origin of the Moon's craters is still a matter of 


debate. Lunar mountains are higher than Earth mountains, 
presumably because of the absence of weathering effects. 


(7) Planetary Satellites 


A number of the satellites of Jupiter, Saturn, and 

Neptune are larger than the Earth's Moon; some may be as 
large as Mercury. At least one of Saturn's moons, Titan, 
is believed to have an atmosphere. Although reliable 
physical data on the satellites are lacking, it is possible 
that some of them may be more hospitable than their 
parent planets. 

There are two principal experimental methods for obtaining data 

concerning the composition of the planetary atmospheres: 
(1) Optical observations of the surface and atmosphere. 


(2) Spectroscopic measurements of radiation emitting from the 
surface, 


Both methods suffer severely because of interference with the Earth's 
atmosphere. 

In addition to the above experimental methods, other constituents 
of the atmospheres of the planets can be inferred by correlating esti- 
mates of the composition of the atmosphere at the time of formation of 
the planet with present measurements of physical characteristics of the 
planet and its atmosphere. Table E.1 gives estimates of composition of 
the atmospheres of the planets and their natural satellites from obser- 
vations and calculations recorded in reference (67). If the amount of 


the constituent is based on spectroscopic measurements, the amount is 
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recorded in Column 4 as the height in centimeters of an equivalent 
column of the gas at 0° C. and Earth standard sea level pressure 

(760 mm Hg). If the amount of the constituent is based on calculations 
and estimates, it is recorded in Colum 3 as a percent by weight of the 


total planetary atmosphere. 






Table E.1: Composition of the Atmospheres of the Planets 
and Their Major Natural Satellites. 












Notes: (1) If the amount of the constituent is based on spectroscopic 
measurements, it is recorded in Column 4 as the height in 

centimeters of an equivalent column of the gas at 0°C. and 

at 2 pressure equal to a 760 mm colum of Hg. 








If the amount of the constituent is based on estimates and 
calculations, it is recorded in Colum 3 as a percent by 
weight of the total planetary atmosphere. 


Column 1 Column 2 Column 3 
Body Gas Amount (4%) 


Venus 








Column 4 


Amount (see Note (1) above) 
cm 









100,000 





£100 
<100 








negligible 
negligible 






625,000 
168,000 
220 

, 


0 
0 


WEN 
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Table E.1 (Cont.): Composition of the Atmosphere of the 
Planets and Their Major Natural Satellites. 





Moon ~<0.0003 
<0.005 

Mars 

Jupiter 


Saturn 

Titan CHy, 33% 20 ,000 
NH} < 300 
A vk 
Ne 


* The total amount of atmosphere for Titan is estimated to be 60,000 cm. 
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Table E.1 (Cont.): Composition of the Atmosphere of the 


Planets and Their Major Natural Satellites. 





Column 1 Column 2 Column 3 Column 4 











Uranus 






220,000 
< 0.1 
< Oger 


370,000 


E.2 Models of the Atmosphere 
In much of the work of this thesis, the following two assumptions 


are made: 

(1) The planet and its atmosphere are spherically symmetric. 

(2) The atmospheric density varies exponentially with altitude. 
The first assumption is reasonable for the terrestrial planets” because 
of their slow rotational speeds. The giant planets have high rotational 
speeds, hence this assumption is not nearly so valid for these planets. 
The second assumption is based on the simple kinetic theory of an iso- 
thermal gas in a uniform gravitational field. This theory yields the 
well-known exponential approximation for the atmosphere: 

-KH 
nes) (E-1) 


* Mercury has no appreciable atmosphere. 


Woy 


where: 


free stream atmospheric density (slug/ft.7) 


“Oo 
I 


Qst)= "assumed" mean sea level atmospheric density” 
K -1 
H 


atmospheric density @cay parameter, ft. 


altitude of vehicle above surface of planet (ft.) 
K is equal to the mean molecular weight of the planet's atmosphere 
times the local acceleration due to gravity divided by the product of 


mean atmospheric temperature and the universal gas constant: 


K = g (Mean molecular weight of atmospheric gases) 
Rh, Timo (E-2) 


>. 


where: 


ee = universal gas constant = 8.31 X 10° ergs/°Kelvin 


g = acceleration due to gravity 
Timo = the mean temperature of the planetary atmosphere (°K) 
Two additional models of the planetary atmosphere are sometimes used: 


(1) An emperical matching of measured atmospheric data: 


free tld (H) (B-3) 


This model is useful for the Earth's atmosphere only, 
Since measurements of the density variation with altitude 
for other planetary atmospheres have not been obtained to 
date. The model represented by Eq. (E-3), with n = -5 to 
-8, generally matches more closely at high altitudes the 


data accumulated from Vanguard, Sputnik, and Explorer 


* This is not the true sea level atmospheric density (cst) . The 
"assumed" quantity (sc) is the intercept of the straight line which 
best fits the curve of log e vs altitude. 
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Earth satellites than does the exponential atmosphere with 
K constant. Since Eq. (E-3) is based on an emperical 
curve fit and not on any known physical law, the value of 
n for other planets cannot be predicted at this time. 

(2) The variation of density with altitude for an atmosphere 


in which temperature decreases linearly with altitude is: 


t 4 
$ oS ee temp. ) }P= Gene temp. )-(Lapse pate) 
( 


SL) Surface temp.)) ~ Surface Temperature 
(E-4) 
where: 
Y= ratio of specific heat ee 


If it is assumed that Y% and gravity are constant, the 


adiabatic lapse rate” is given by: 


(Mean Molecular Weight of atm. gases) 9im)o (%- 1) 


(Adiabatic Lapse Rate) = 
Ry % 


(E-5) 


E.3 The Venusian Model Atmosphere 
It is generally assumed that the atmosphere of Venus originated 


from volcanic gases. On Venus, water was converted to molecular oxygen, 
atomic oxygen, ozone, and hydrogen via photodissociation at high alti- 
tudes. This action was speeded up because of the nearness of Venus to 
the Sun. The ozone layer on Venus was probably formed at greater alti- 
tudes than on Earth. The oxygen was probably removed from the atmosphere 
by chemical combination with surface material. This chemical process 


* Lapse rate is the rate at which temperature decreases with increasing 
altitude. 
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should react about 30 times as fast on Venus as it does on Earth 
because of the higher surface temperatures on Venus. 

Loss of water via photodissociation to oxygen probably stopped 
the chemical reaction of carbon dioxide with surface material. Once 
this situation prevailed, it may be expected that carbon dioxide started 
to appear in the atmosphere of Venus in ever-increasing quantities. 

Venus should have about the same amount of nitrogen as Earth, and 
much greater quantities of carbon dioxide in the atmosphere. The much 
larger quantities of carbon dioxide are expected in the atmosphere of 
Venus when compred to Earth because of the high quantities of carbon 
dioxide combined with surface material on Earth.* 

Because of the above reasoning, Dore 668) proposed on atmospheric 
model containing 90% carbon dioxide and 10% nitrogen. The values of 


the constants in Equation (E-1) for Venus proposed by Dole are: 


Qa= 0.0326 slug/ft.- 


L) 


mot, — 
K = 4.88 X10 ~ ft. 
Romer‘ 72) , using Equation (E-4), determined Pcs.) = 0.0348 slug/ft.-. 
He used an atmospheric model of 10% nitrogen and 90% carbon dioxide 
with a mean molecular weight of 42.4 and a specific heat ratio of 
% = 1.31 at 400 °K. The lapse rate was computed to be 10.45 °K/KM. 


(8) suggests the lapse rate is 13° K/KMand % =1.4, while 


Kuiper 
aeeneculeums OC]? dedermines/s lapse rate of 10.4 °K/KM. 


Values given in (E-6) are used in this thesis for Venus. 


* Carbon dioxide was able to combine with surface material on Earth 
because of the presence of water. 


er? 


E.4 Martian Model Atmosphere 


Evidence for the existence of a Martian atmosphere is conclusive. 
Clouds are frequently observed. Stars passing behind the disc of Mars 
are seen longer and appear earlier than possible without atmospheric 
refraction’. A twilight zone on Mars caused by atmospheric refraction 
has been observed, Mars appears to be larger when photographed by 
reflected ultraviolet light than by reflected infrared light. (The 
atmosphere of Mars is included in ultraviolet photographs). A lower 
limit of 60 miles can be placed on the altitude of the Martian atmos- 
phere from these photographs. This is probably the altitude to which 
small ice crystals rise, thus creating a "blue haze". The blue haze 
often clears when Mars passes close to the Sun. On these occasions, 
damage is seen to occur to the primitive vegetation which is believed 
to exist on Mars. When present, the blue haze protects the vegetation 
from destructive ultraviolet radiation originating from the Sun. 

Although nitrogen is probably the main constituent of the Martian 
atmosphere, it has not yet been detected, This is due to the fact that 
absorption bands for molecular nitrogen are in the low ultraviolet 
region, These lines are effectively obscured from observation on Earth 
by ozone in the upper atmosphere of Earth. Atomic nitrogen in the 
Martian atmosphere could theoretically be detected through the atmosphere 


* It is interesting to note that the only time in recorded history that 
occultation of a first magnitude star (Regulus) by Venus occurred in 
mid-July, 1959. The total time that Regulus was obscured was ll 
minutes 4.8 seconds. This event was observed from a Spanish obser- 
vatory by a team led by Dr. Allen Hynek from the Smithsonian Astro- 
physical Observatory of Cambridge, Mass. Revised information on the 
Venusian atmosphere (its density variation with altitude, temperature, 
chemical composition, etc.) based on the results of this observation 
has not been published at this time. 
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of the Barth. However, there is a larger quantity of atomic nitrogen 
in the atmosphere of the Earth than in the atmosphere of Mars. Mars 
probably had less molecular nitrogen to begin with than Earth, and 
definitely has less solar energy to produce atomic nitrogen from molec- 
ular nitrogen. These two factors have combined to provide a weak spec- 
trum against a background that even sensitive Doppler techniques have 
proven futile for the detection of atomic nitrogen. 

It is generally thought that the atmosphere of Mars, like that of 
Venus and Earth, originated from volcanic gases. These gases consisted 
of water vapor, carbon dioxide, nitrogen, and other minor constituents 
in that order. Nitrogen was probably the only major constituent of 
volcanic gas which was not readily removed from the atmosphere of Mars 
by chemical reaction with rock formations or by escape since nitrogen 
is inert and highly resistant to photodissociation. 

Water was probably removed first from the Martian atmosphere. 
Atomic oxygen, ozone, molecular oxygen, and hydrogen are some of the 
gases produced as a result of photodissociation of water. Most of the 
hydrogen probably escaped from the Martian atmosphere. Since ozone would 
be produced at low altitudes on Mars, it could combine rapidly with 
surface material. Water probably remained close to the surface making 
the rapid combination of carbon dioxide with surface material possible. 
With the loss of oxygen via ozone formation and chemical reaction, there 
would result a shift in the chemical equilibrium to more ozone; this gas, 
in turn, could continue to combine with surface material until most of 
the water disappeared. 

After the disappearance of most of the water from the Martian 


atmosphere, carbon dioxide could no longer combine with surface material 
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at an appreciable rate. Consequently, it may be reasoned that the 
atmosphere of Mars has more carbon dioxide than the atmosphere of 
Earth. This has been verified by spectroscopic measurements (see Table 
E.1). In addition, water and oxygen are known to be absent from the 
Martian atmosphere in quantities which are significant in gas-dynamic 
problems (less than a few hundredths of a per cent.) The reasoning 
outlined in the foregoing paragraphs may account qualitatively for the 
composition listed in Table E.1. 


The Gazley‘t4) 


model for the atmosphere of Mars assumes that it 
consists of 95% nitrogen and 5% carbon dioxide. This model ignores 
minor constituents, especially argon, and is subject to a large error in 
the estimation of carbon dioxide content. Nevertheless, it is doubtful 
that a more accurate model would allow a significant improvement in the 
accuracy of gas-dynamic calculations. 
The average surface temperature of Mars was estimated by Gazley 
from the Stefan-Boltzman relationship to be minus 40° F.* He made 
the following assumptions in determining the average surface temper- 
ature: 
(1) The ratio of the heat flux reaching Mars and Earth is 
determined by their respective solar constants C7 Ter 
Mars and 1.0 for Earth). 


(2) The emissivities are the same for both planets. 


(3) The temperature of the surroundings to which the planet 
radiates is sufficiently low to be neglected. 


It is noted that a value of ~-40° F. was listed by Wanders (7) as an 


- Kuiper ‘67? and de Mabcodleurs\”-> use a surface temperature of 
273° K. (0° C.) The maximum diurnal temperature change is about 
50° ©. This is also approximately the maximum seasonal average 
temperature variation (occurring in summer and fall). 
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average experimental value. 

The ratio of the absolute mean-sea-level temperature of Earth to 
the absolute ground-level temperature of Mars is 1.23. The temperature 
lapse rate of the Martian atmosphere is lower than for the Earth” ; 
therefore, this temperature ratio decreases with increasing altitude. 

Gazley ‘14? uses a surface density value of 6.2 X 1077 Lbs. /ft.? 
for the Martian atmosphere which is in close agreement with the value 
of 1.94 X ies” slugs/ft.? given oar tte ene Gazley uses a value for 
K (see equations E-l and E-2) of 1.15 X ar ft.7t, Chapman‘15? uses a 
value for K of 1.665 X 107? ft.-+, In the numerical calculations of 
this thesis, the following are used for the values of the constants in 


Eq. (E-1) for the Martian atmosphere: 


~4 


des OSX 1,0 slugs/ft.7 


II 


PSL) 
K 


(E-7) 
1.665 X 107 ft.7+ 


II 


The maximum error in decelerations computed in this thesis should be no 
greater than about 20% for Mars, and the maximum error in heating rates 


and stagnation point temperatures should not exceed 10%. 


E.5 Model Atmosphere for the Giant Planets and Titan 


Information on the density characteristics of the atmospheres of 
Jupiter, Saturn, Titan, Uranus, and Neptune is fragmentary at best. 
Because of the cloud layer which obscures the surface of the giant 
planets, extrapolation of density characteristics to the surface is not 


* The lapse rate for Earth at sea level is 9.8 °K/KM. The average 
measured lapse rate at all altitudes is approximately 5 °K/KM, which 
indicates that the atmosphere ot) eee in adiabatic 


equilibrium, Most authorities have computed a lapse rate of 
3.7 - 3.9 °K/KM for Mars. 
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possible. Table E.2 summarizes briefly Kuiper's model atmospheres for 


the giant planets and Titan. 


Table E.2: Model Atmospheres far Giant Planets and Titan 


1. Jupiter 
Height of cloud tops from main cloud deck: 20-30 KM. 


Exponential decay parameter: K =1.65 X 107? ft.~1, 


Pressure at main cloud deck: 6-9 atmospheres. 


Temperature at top of atmosphere: 80-90°K. 


Atmospheric Model I: % = 1.46 
Mean molecular weight 2.47. 
Lapse rate 2.64 °K/KM 


Atmospheric Model II: % = 1.256 
Mean molecular weight 3.26 
Lapse rate 3.96 °K/KM 


Temperature at top of clouds: Model I: 165° K 
Medel If: 170° K 


Temperature at main cloud deck: Model I: 215-245°K 
Model II: 250-290°K 





2. Saturn 
Height of cloud tops from main cloud deck: 20-30 KM 


Temperature at top of atmosphere: 60-70° K. 


Atmospheric Model I: 7 A fs | 
Mean molecular weight 2.47 
Lapse rate 1.17 °K/KM 


Atmospheric Model II: Fo = 1.56 
Mean molecular weight 3.26 
Lapse rate 1.75 °K/KM 


Temperature at top of clouds: 155° K. 


Temperature at main cloud deck: Model I: 175-190° K. 
Model II: 190-210° K. 
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Table E.2: (Cont.) Model Atmospheres for Giant Planets and Titan. 


3. Uranus and Neptune: 
Mean molecular weight 3.55. 


Mean temperature of atmosphere: 78 °K. 
Pressure at main cloud deck: 9 atmospheres 


Exponential decay parameter: K = 1.87 X ie ft.7+ 


4, Titan 
% = 1.5 
Mean molecular weight: 20 
Mean temperature of atmosphere: 70-83° K. 
Lapse rate 1.4 °K/KM 


Pressure at surface: 8.9 X 107" atmospheres 


Osye 0.006 slugs/ft.- 





E.6 Model Atmosphere for Earth 

The 1956 ARDC model atmosphere is generally accepted as a reason- 
ably good representation of average atmospheric characteristics of the 
Earth below 100 miles. This profile was constructed with only one 
point above 100 miles. Data obtained from various sounding rockets 
and Earth satellites during 1957-1959, however, suggests that substan- 
tial revision of this model of the atmosphere is needed. Revised data 
obtained from satellites indicates that density of the ARDC model may 
be in error by an order of magnitude at low altitudes, and that den- 


sities at altitudes in excess of 100 miles are much higher than predicted 
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by the ARDC model. Furthermore, satellite data shows that the density 
at high altitudes is approximately twice as high in summer as in winter, 
and twice as high in the daytime as at night. Therefore, any model of 
the atmosphere is at best, a crude representation of conditions that may 
actually be encountered by a vehicle entering the atmosphere. 

Fig. E.l compares atmospheric density at high altitudes predicted 
by the ARDC model atmosphere with actual measurements made by recent 
satellites and sounding poseete A mean curve is shown on this 
figure which more closely matches measured densities at high altitudes 
than does the ARDC model; this mean curve matches the ARDC model at 
altitudes below 90 KM, and does not depart much from the ARDC model 
below 120 KM. 

The exponential approximation to the Earth's atmosphere was used 
in most of the numerical computations of this thesis. Values used in 
Equation (E-1) for Earth are: 


(5) 


K 


0.0027 slug/ft.? 


1/23,500 ft.7_ 


it 


(E-8) 


il 


The measured sea level atmospheric density of the Earth is 2.38 X 107? 


slug/ft.? 

The exponential approximation for Earth using the values given in 
Eq. (E-8) are compared in Fig. E.2 with the 1956 ARDC Model atmosphere. 
It is clear that a single value of K is a reasonable approximation 


below 400,060 ft. (approximately 80 miles)*. Numerical comparison of 


* It is shown in ref. (15) that peak decelerations and maximum aero- 
dynamic heating of an entry vehicle occurs well below this altitude. 
The region of most important heating and deceleration for a given 
vehicle occurs over a strip of altitude approximately 70,000 ft. 
thick across with the density changes by a factor of 20. Once the 
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the exponential model and the ARDC model is given in Table E.3 to four 
or five significant figures. 

Most of the derivations and numerical calculations of this thesis 
use dimensionless parameters in order that the results may be applied 
to any planet. Mean planetary radius is used as the reference length 
for non-dimensionalizing quantities having the dimensions of length or 


distance. Thus: 


—&h 
O = = - 
+. e (E-9) 


where: . 


k = dimensionless atmospheric decay parameter = KROn)O (B-10) 


h 


dimensionless altitude = H/R mo (E-11) 


G6 = atmospheric density ratio (E-12) 

Fig. E.3 shows a plot of the dimensionless exponential decay 
parameter Jk taken from ref. (15). k was determined for the ARDC 
model atmosphere at each altitude point by determining the mean slope 
eeeaeplot of log p vs altitude for the 70,000 ft. strip of altitude 
immediately above the particular altitude in question. It is seen 
from Fig. E.3 that below 400,000 ft. the variation in /k is no 
more than about 10% from a mean value of 30. These variations are due 
primarily to temperature changes with altitude (a constant exponential 
decay rate assumes an isothermal atmosphere). Since temperature changes 
as much as 15% with season and latitude(?”), than Jk (which is pro- 


altitude of this critical strip is determined for a given vehicle 
entering a given atmosphere, the exponential decay parameter K can 
be adjusted to correspond more nearly to that of the critical strip 
of altitude. 
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Dimensionless Parameter 
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portrorial to Blica\0 , see equation (E-2)) may fluctuate by 7% 
due to seasonal changes and latitude changes alone. 

For calculations of this thesis where the Earth is considered, 
k = 889 is used (corresponding to a mean atmospheric temperature of 


240° K (432° R)) together with © Pg, = 0.0027 slug/ft.’. 


E.7 Comparison of the Exponential Model Atmospheres of Venus, Earth, and Mars 

Fig. E.4 shows the isothermal atmospheric models used in this thesis 
for the terrestrial planets. It is apparent that the atmospheres of 
Venus and Earth are distributed in a similar way. Even with large errors 
in the assumed atmospheric eee of Venus, gas-dynamic heating and 
deceleration loads would, to a first approximation, be the same as they 
are on the Earth. A parachute on Venus, however, should be about three 
times as effective as on Earth. 

Maximum deceleration loads will be severe on Earth and Venus for 
some trajectories as pointed out in this thesis. These problems may be 
solved by resorting to lifting vehicles, by entering at lower entry angles 
(requires precision guidance) or coming in at slower velocities 
( requires large amounts of propellants). 

The atmosphere of Mars is characterized by a gradual decrease in 
density with increasing altitude and low surface density as compared to 
that of the Earth. The low ground-level density of this atmosphere will 
reduce the effectiveness of a parachute and will require the use of a 
high drag (low density) entry vehicle. The gradual change in atmospheric 
density of the Martian atmosphere will reduce the severity of peak gas- 
dynamic accelerations and heating loads. For identical size, weight, 


shape, and velocity, a body entering the Martian atmosphere would be 
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subjected to about 25% of the peak deceleration load and about 80% 
of the peak heating load that it would encounter during entry into the 


Earth's atmosphere. 


Table E.3: Comparison of ARDC Model Atmosphere (1956) 


with Exponential Approximation ry 
( Ast 22.7072) (UK = 2)395eh fo”) 


H (feet) ARDC (slug/ft.2) 


2.3769 x 1072 
1.7556 x 102 
132673 x 1072 
8.9068 x 1074 
5.8727 x 107+ 
3.6391 x a 
292561 ae 107 
1.3993 x 107% 
8.6831 x 107 
5.2531 x 107 
3.2114 .% 107) 
10% 2.0014 x 10-5 
Lo; 1.2697 x 1079 
10," 8.1894 x 10-6 
10 5.3640 x 1076 
10+ 3.5642 x 1076 


10;¢ 2.4329 x 10-6 
10 1.6929 x 10-© 
104 1.1995 x 10-6 


10} 8.5890 x 107? 
10} 6.0583 x 107 
10) 2.8710 x 107? 
10; 1.2580 107% 
=. me ‘ 
Zo x I 2702 
30 x 10° 6.065 
32.5 x 10} 1.610 
Boeolx ie 4.957 
37.5 x 10 1.718 
ho.0 x 107 6.565 
42.5 x 10), 26 
45.0 x 10 Lela 
W765 x 10} 5.348 
50.0 x 10 2 862 
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APPENDIX F 


EXTERNALLY-AIDED SELF ADAPTIVE CONTROL OF THE ENTRY VEHICLE 


The operational extremes experienced by a vehicle during entry into 
a planetary atmosphere suggest a control system that adapts itself to 
the changing environment. The Conservation Parameter, defined in 
section 1.7, is a measurable function that has sharp behavior at the 
boundaries between operational phases; behavior of this function is 
discussed in some detail in Chapters 8 and 9. The conservation parameter 
may be useful as a switching function and as a variable sensitivity factor 
to improve control system operation. 

The generalized guidance and control system represented functionally 
in Fig. 1.4 is a three-dimensional system. The command signal to the 
control system from the guidance computer may be in the form of accel- 
eration components along three axes fixed to the entry vehicle. Fig. 

F.1 suggests a relatively simple functional instrumentation of the control 
system for one of these axes, the normal axis, which utilizes the 
conservation parameter for adapting to the operational environment. 
Similar adaptive techniques may be applied to other control channels. 

The function of the acceleration control system of Fig. F.1 is to 
produce normal accelerations of the vehicle in response to command 


inputs, while minimizing changes in flight direction that result from 
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gust disturbances. The operation of this system is not unlike that used 
to control guided missiles and high performance aircraft except for the 
unique feature of incorporating the conservation parameter as an input. 
The environmental logic system compensates for environmental influences 
on the effectiveness of the elevators in controlling the vehicle. 

The normal specific force measuring system in the feedback loop 
detects the sum of all specific force components in the normal direction. 
With measurements of bank and pitch angles and altitude, it is possible 
to introduce compensations for gravity. Such compensation is not shown 
ay Pie, Fl. 

After compensating for gravity, the measured normal acceleration 
is compared to the input command voltage. When a difference exists, a 
correction is introduced. The correction signal goes to the environ- 
mental logic system and the integrating gyro unit. The integrating gyro 
unit produces a signal containing two terms, one proportional to the 
integral of the correction voltage and one proportional to the integral 
of the pitch rate. In order to produce a large elevator angle for 
given pitch deviation such as may be required in the subsonic phase of 
flight near the termination of the entry trajectory, it is necessary to 
provide a signal to the elevator servo proportional to pitch rate 78) | 

The environmental logic system takes the time derivative of the 
Signal from the integrating gyro unit and combines it with the correction 
voltage. The output voltage of the environmental logic system is sent 
either to the reaction controls, to the elevator servo, or to both, 
depending on the magnitude of the conservation parameter and its time 
behavior. The conservation parameter, which comes from the guidance 


computer, indicates the phase of flight (see Chapters 8 and 9) and 
Ou, 


reflects the degree of effectiveness that can be expected from the 
elevators due to flight conditions. 

The system of Fig. F.1 is not a self-adaptive control system, 
i1.e., it does not monitor its own performance and automatically adjust 
itself for optimum performance in a changing environment. The system is 
"adaptive", however, in the sense that it uses environmental data from 
the guidance computer to adjust its performance characteristics. 
Augmenting the adaptive features of a system similar to Fig. F.1 with 
self-adaptive inner loops to optimize control system operation is 


suggested for future investigations. 
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